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package edu . st anf o rd . f acs . 1 ogl i ke ; 

import java.awt.*; 
import [java.awt.geom. 
import lavax. swing.*; 
import javax. swing. border.*; 

/it it 

* <p>Title: Log Like Data Transform</p> 

* <p>Description: </P> ^ . . 

* <S>Copyright: copyright. (c) 2002</p> 

* <p>company: Stanford umversity</p> 

* ©author wayne A. Moore 

* ©version 1.0 
*/ 

public class AxisDialog extends JDialog 

' ^ Zf£r^i''£^^orT^^^^ BorderLayoutO; 

double[][] nomoRelatiye; _ „ , , ,._ „ .._ qs", "-.02", "-.01". "O 

Stringn nornoLabel,= { -1 , -J... .'.2.. / "5.'. , "lO". "20", "50". "100 . 

.01 i -vZ J -vJ , ■ J- » 

"200", "500" }; 

private : Panel q Panel 1 = new ^ Panel O; 
private DLabel jLabell = new ^LabelQ, 
nr-iv3te 3Label iLabel2 = new DLabel Oi 
?ri^l« GridiiyJut gridLayoutl = new GndLayoutQ ; 
nrHvatp ipanel i panel 2 = new DPanelUi ^ 
griJI"! DWonint axi scomponent = new ^Component 0 

^ public void pai ntcomponent ( 
Graphics graphics ) 

^ MTl JSlSidtSo'-'inlPtileft - inset.right; 
in^ h = gSHeightO - inset. top - inset.bottom; 
int xO = inset. left + w/3; 
int xl = inset. left + 2*w/3; 

iSgSni-.m^Rl??^ ?rsS?^lf??" e '.top, w, h 

graphics. drawLineC xl, inset. top, xl, inset. top + 
if (nomoRelative == null) 

int^^yo"clip = inset. top, yl_clip = i^set-top; 

int ascent = graphics. getFontMetncsQ .getAscentQ , 

for (int 1 = nomoRelative. length - 1; i >= 0 . - t) 

^ -ini- vn - inset too + h - (int) Math, round ( notnoRelative[i] [0] * h ); 
in? yl = Ynsfx.'Xol X h - (int) Math.round( nomoRelative[i] [1] h ), 

'RK?Sngl42rri??^4ra^hics!g^^^^^^ .getStringBounds( nomoLabel [i] , 

^'■^'''''"if'(yO_clip<yO- ascent) 

^ graphics.drawstring( nomoLabel [i ] , xO - em/4 - (int) rect.getwidthQ . yO 

^' graph! cs.drawLine( xO, yO, ^^^^^^l^' ^' 



h - 1 ); 
h - 1 ); 



parks_et_al_l.txt 

yO_clip = yO; 
if (yl_clip < yl - ascent) 

yl_clip = yl; 

} 

} 

J^late BorderLayout borderuay^^^^^^^^^ BorderLayoutQ ; 
private 3 Label 3 Label 3 = new DLabelU, 

?Slfc%5?s'D?alogCF?ai4 fra.e. String title, boolean modal) 

^ super (frame, title, modal); 
try 

^ jbinitO; 
packO ; 

catch (Exception ex) 
^ ex.printstackxraceO ; 
} 

} 

public AxisDialogO 

^ this (null, "", false); 

private void jbinitO throws Exception 

',tSfl2:lSrriSiaufg^™ent(swin,constants.cENTEK); 

1 Label 2 . setText (" Logi cl e ) ; 

iPanel2 setLayout(borderLayout2), . 
^plSIll"sltMinimumSize(new Dimension(60, 17)) , _ 
?S!!lli:StPreferredsize(new Dimension (180, 17)), 
gri dLayoutl . setcol umns (3) ; 
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package edu.stanford.facs.loglike; 

import java.awt.*; 
import java.awt, geom.*, 
import javax. swing. , ^. 
import 3avax. swing. border. , 

/r<p>Title-. Log Like Data Transform</p> 

* <P>Description: </P> 2002</p> 

:Mrnfst»^ 

* ©author Wayne A. Moore 

* ©version 1.0 

public class Ax-isDialog oxtends lioialog 

"200"', "500"'}; 

^ public void paintcomponent ( 
Graphics graphics ) 

S R5^cl:Str„l^ fi; HIS?; xi, mset.top . h - 1 ), 

IcnomoRelative == null; 

return; . ^-i (-lin = inset. top; 

int yO_clip = ^"set.top, yi-cnp ^^^^^ ^^^^^^^ 

id'^urti-cF^rajSic^Vt^^^^ 

»"^^'"jV(y0.cl1p < yO - ascent, _ ,,,,,etwidtKO , VO 

' graphics.drawStrinsC noBOLabeUi], xO em/i 

>• graphics.drawLineC xO, yO, xO^- em/8, yO); 
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yO_clip = yO; 
if (yl_clip < yl - ascent) 

^ graph! cs.drawStringC nomoLabel[i], xl + em/4, yl ); 
graphics. drawLineC xl, yl, xl + em/8, yl ), 
yl_c1ip = yl; 

} 

} 

lUvate BorderLayout borderLayoutZ = new BorderLayoutQ ; 
private J Label i Label 3 = new JLabelQ; 
';i]^?^^'^U^ogS"^ fran,e, string title, boolean modal) 

^ superCframe, title, modal); 
try 

jbinitO; 
packO; 

catch (Exception ex) 

^ ex.printStackTraceO; 
} 

} 

public AxisDialogO 

^ this (null, false); 

private void jbinitO throws Exception 

BoLe?S?;Jr;.createCompoundBorderCBorderFac^^^ 

ED, col or. white, col or. white, new color(103, 101, 98), new ColorU4B, 

140)) BorderFactory . createEmptyBorder C6 , 6 , b , b;; , 

panell.setLayout(borderLayoutl); . ^ ^ ^cm-tcd^. 

1 Label 1. setHori zontal Al ignmentCSwi ngConstants . CENTER) , 

i Label l.setTextC'FlowDo"); , . ^ ^ ^,.K,-rcD^. 

1 Label 2 . setHori zontal Al i qnment (Swi ngConstants . CENTER) , 

jLabel2.setText("Logicle ); 

1 panel l.setLayout(gridLayoutl); 

axi SComponent. setBackgroundCCol or. white); 

axi SComponent. setMi nimumsize (new Dimension(60, 120)) , 

axi SComponent. setPreferredSize(new Dimension(180, 360)), 

jPanel2.setLayout(b9rderLayout2); _ 

1 Panel 1 . setMi ni mumSi ze (new Di mensi o" CoO^ on 1 4^ ^ . 

iPanell.setPreferredSize(new Dimension (180, 17)), 

gridLayoutl.setColumns(3) ; 

panel 1. setBorder(borderl) ; . 

j Panel 2 . setBackground (Col or . whi te) ; 

getContentPaneO .add(panell) ; ,.«„™>. 

panel 1 . add ( j Panel 1 , BorderLayout . NORTH) , 

i Panel 1. add (i Label 1, null); 

iPanell.add(3Label3, null); 

j panel 1. add (j Label 2, null); 

panel 1 . add ( j Panel 2 , Borde r Layout . center) , _ 
1 panel 2. add (axi SComponent, BorderLayout. center) , 

} 

} 
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package edu . Stanford . f acs . 1 ogl i ke ; 

import java.awt.*; 
imoort i ava . awt . event , 
impSrt java.io.BufferedReader; 
import qava.io.Fi1e; 
import lava.io.FileReader, 
import iava.io.iOException; 
import lava.util .Arrays; 
import iava.util.ArrayList; 
import lava. util. Iterator; 
imoort lava.util .Map; 
ilport jiva.util .StringTokemzer; 
import 3avax. swing.*; 
import javax. swing. table. , 



import org. i sac.*; 



^**<p>Tit1e: Log Like Data Transform</p> 

: SSop5?iShi:"^o;^?ight (c) 2002</p> 
* :?;Span?: Stanford university</p> 



* ©author Wayne A. Moore 

* ©version 1.0 

*/ 



public class converterPrame extends JFrame { 

Jlra^e r«tSix?rSS; S.Tl2l^lfr §;.«enux«™0 = 

nr-ivate JToolBar jToolBar = new DTooiBaru, 

nAyltl JButton jButtonl = new JButtonQ 

lAllle jiStton iButtonZ = new ^ButtonQ 

lAwlll llutton ^Buttons = new DButtonQ; 

private imageicon imagel; 

private imageicon imageZ; 

private imageicon imagei; -..^u^ir^- 

//construct the frame 
public converterFrameU 

^ enableEventsCAWTEvent .WINDOW^EVENT^MASK) ; 



try 

jbinitO ; 



catch (Exception e) 

^ e.printStackTraceO; 
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//component initialization ^-^^ r 

private void jbinitO throws Exception { 

imagSU""anford . facs J ogl i ke . converter Frame . cl ^ 

LglSn(edu"t?anford.facs.log1ikeXonverterFrame.class.get^^^^^ 

W%ldu?l?anford.facs.1oglikeXonverterFrame.class.getResou^ 
//seticonimageCToolkit.getDefaultToolkitO.createimageCconverterFrame.class.getResou 

'""^^Kt'panS^Po^anel) this.getContentPaneO; 
content Pane . setLay9ut (borde r Layoutl) ; 
this.setSizeCnew Dimension (400, 300)), 
this.setTitleC'Log Like Data Converter ), 
statusBar.setTextC; "); 
iMenuFile.setTextC File ); _ 

TMenuFileExit.setTextC Exit ); . . ^ /-^ / 

jSeSuFiltixit.addActionListenerCnew ActionustenerQ { 
public void actionPerformed(ActionEvent ej I 
jMenuFileExit_actionPerformed(e; ; 

} 

}) ■ 

iMenuHel p . setText ("Hel p") ; 

iMenuHelpAbout.setTextC About ); ci-<>nprr^ i 

] MenuHel ?About . addActi onLi stener (new Acti onLi stener () i 
public void actionPerformed(ActionEvent e) I 
j MenuHel pAbout_actionPerformed(e) ; 

} 

pSi:lddic?S^ 

public void acti on Performed (Acti onEvent e) 

^ data_action(e) ; 

} 

iButtonl.setToolTipTextC'Open File") ; 

i ButtonZ . setlcon (i mageZ) ; ^ . , ... . 

iButtonZ.setToolTipText ("Close File ), 

1 Buttons. seticon(image3); , 
j Button 3. setToolTipText( Help ); 

:iSe«™;ld5Ic"^nSffe™5rC^4w java.awt event.Acti i 
public void actionPerformed(ActionEvent e) { 
spectrum_action(e) ; 

} 

jSreSirddlcJiffiSPni^ 

public void acti onPerformed(Acti onEvent e) 
{ 

data_action(e); 

} 

jMenuitem2 . setText("Save") ; 
iTool Bar . add (q Buttonl) ; 
iTool Bar . add (j Button2) ; 
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iToolBar.addCjButtonS); 
TMenuFile.addCmenuSpectrum;; 
iMenuFile.add(jMenuiteml) ; 
]MenuFile.add(DMenuitem2); 
iMenuFi 1 e . add (iMenuFi 1 eExi t) ; 
1 MenuHel p . add (3 MenuHel pAbout) ; 
]MenuBarl.add(iMenuFile) ; 
jMenuBarl.add(]MenuHelp; ; 

this.setDMenuBar(jMenuBarl); M«B-ru\. 
JontentPane.addCjToolBar, Border-Layout. NORTH) 
contentPane.addCstatusBar, fo rderLayout. SOUTH) , _ 

lilWl iof5lMlS^F??eixu!a«?onPerfo™ed(ActionEve« e) ^ 
System. exit (0) ; 

i„tS%S55;iSAK«!™?ionPerfo™edCActionEvent e) { 

//overridden so we can exit when window is closed 
protected void processWindowEvent(WindowEvent e) I 
super. processwi ndowEvent (e) ; . ^. j 

if (e getlDO == wi ndowEvent. window closing) { 
jwenuFi 1 eExi t_acti onPerf ormed(nu 1 1 ) ; 

} ^ 

vate static final String[] stringArray = new string[0]; 
vate FCSFile fcs; 

] sensor_names; 
int_data; 
raw_data ; 
compensated_data ; 
compensati on_mat ri x ; 



pri 
pri 
pri 
pri 
pri 



vate String 
vate int[][ 



^. .vate double __ 
private double [! 
private double""" 

void compensate () 
throws lOException 

^ int n = fcs.getTotalQ; 
int m = sensor_names. length; 



raw^data = new double [ ") ][ "^^^ . . 
compensated_data = new double [ m ]L n J, 
for (int i = 0; i < m; 

^ FCSParameter p = fcs .getParameterC sensor names [i] ); 
int[] input_data = int„data[p.getindexO-l] , 
doubled output_data = raw_data[ij; 
if (p-isLogO; 

^ doubled antilog = new double [ p.getRangeQ ]; 
antiioaFOl = p.qetMinimumO ; 

iSter?5lateC antilog, 0, p.getRangeQ, p.getMaximumQ 
for Cint i = 0; j < n; ++t) 
output_data[j] = antilog[ input_data[]] J; 

} 

else 

^ double g = p.getGainO; 

for (int j = 0; j < n; ++i) 
output_data[j] = g * input_data[3] ; 
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} 

} 

for (int i = 0; i < m; 

^ doubleG mult = compensation_matrix[i] ; 
doublet] output_data = compensated_data[i] , 
for (int j = 0; j < n; ++J) 

double X = 0; 

for (int k = 0; k < tn; ++k)., . 

X += raw_data[k][]] * mult[k]; 
output_data[j] = x; 

} 

} 

for (int i = 0; i < m; 

^ FCSParameter p = fcs.getParameteK sensor_natnes[i] ); 
double[] input_data = compensated_dataLi J , 
int[] transformed_data = new intLnj; 
if (p.isLogO) 

^ doubled lookup = new double[ p.getRangeO ]; 

{S?^"r?ilL^C^Srp?'o';™^^^etRang.O. p.,etMaxi.u.O ); 

for (int j = 0; j < n; 

^ int k = Arrays.binarysearch( lookup, input_data[j] ); 

if (k < 0) ^ ^ ^.^ . 
transformed_dataLDJ = -k; 

^ transformed_data[j] = k; 

} 

} 

else 

^ double g = p.getGainO; 
'°Jr^isforjed.dltalj] =' (i^t) Math.floor( input_data[3] / g ); 

} 

} 

} 

void spectrum_action ( 
ActionEvent e ) 

^ int returnval = spectrumchooser.showOpenDialog( this ); 
if(returnval == jFileChooser.APPROVE_OPTlON) 

{ 

try 

^ BufferedReader br = new BufferedReader( new FileReader( 
spectrumchooser .getSelectedFileO ) ); 
br.readLineO ; 

sJn-ng?oken^^^^^ st = new stringTokenizer( br.readLineO. "\t" ): 
ArrayList al = new ArravustO; 
while (st.hasMoreTokensO) 

fSr (int r = 0; r < sensor_nanies. length; ++r) 
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^ St = new stringTokenizerC br.readLineO , "\t" ); 
for (int c = 0; c < columns. length; ++c; 

^ rows[r][c] = new DoubleC Double. parseDoubleC st.nextTokenQ ) ); 
} 

} 

for fint c = 0; c < columns, length; ++c) .... 
fows[sensor_names. length ] [c] = new DoubleC 0 ); 

spectrum_table.setModel( new oefaultTableModel ( rows, columns ) ); 

} 



/ 



} 



catch (Exception ex) 

^ ex.printStackTraceO; 
} 



compensation_matrix = new double[ sensor names length ][ sensor.names. length ] 
for (int i = 0; i < sensor.names. length; ++^) 
compensation_matrix[i][i] =1.0; 



} 

private void interpolate ( 
double [] func, 
int i , 
int n, 
double x ) 

^ . 

int m=n/2; r-T*«^. 
double y = Math.sqrt( func[i] * x ), 

func[ i + m ] = y; 

if (m > 1 ) 

^ interpolate( func, i, m, y ); 
interpolate^ func, i + m, m, x ;; 

} 

} 

void data_action(ActionEvent e) 

^ int returnval = spectrumChooser.showOpenDialog( this ); 
if(returnval == jFileChooser.APPROVE_OPTlON) 

{ 

try 

^ FCSFile fcs = new FCSFile( spectrumchooser .getselectedFileO ); 
int n = fcs.getTotalO; 
int m = fcs.getParametersO; 
int [in int_value = new int[m][n]; 
FCSHandler Imi = fcs .getinputiteratorQ ; 
for Cint i = 0; Imi .RasMoreEventsO ; 

for Cint i = 0; Imi .hasMorevaluesO ; ++:) 
int_value[j][i] = Imi . readvalueO ; 
this. fcs = fcs; 

catch (Exception ex) 

^ ex.printStackTraceO; 
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package edu . Stanford . f acs . 1 ogl i ke ; 

import :java.awt.*; 
import lava.awt. event.*; 
import ^ava.io.BufferedReader; 
import lava.io.File; 
import lava.io.FileReader; 
import ;java.io.lOException; 
import lava.util -Arrays; 
import ]ava.util .ArrayList; 
import java.util .Iterator; 
import java.util -Map; 
import java.util .St ringTokemzer; 
import lavax. swing.*; 
import javax. swing. table.*; 

import org.isac.*; 



parks«et«al_4.txt 



* 



<p>Title: Log Like Data Transform</p> 

* <p>Description: </p> 

* <p>Copyright: copyright Cc) 2002</p> 

* <p>Company: Stanford university</p> 

* ©author Wayne A. Moore 

* ©version 1.0 
V 

public class converterFrame extends 3Frame { 
private D Panel content Pane; 
private JMenuBar jMenuBarl = new JMenuBarQ; 
orivate JMenu jMenuFile = new DMenuQ; _ 
brivate DMenuitem jMenuFileExit = new DMenultemQ; 
orivate DMenu jMenuHelp = new JMenuQ; 
private DMenultem jMenuHelpAbout = new DMenultemO , 
private JToolBar jToolBar = new JToolBarCj; 
private DButton jButtonl = new JButtonO; 
private JButton iButton2 = new DButtonU; 
private JButton ]Button3 = new 3ButtonC); 
private imageicon imagel; 
private imageicon image2; 
private imaqeicon image3; u ir^. 

nn'vatp iLaBel statusBar = new JLabeKJ, 
IrVvlTe Bo?dlrL5oJlt borderLayoutl = new BorderLayout () ; 
pHvate DScrollPane jscrollPanel = new 3ScrollPaneO ; 
private DTable spectrutn_table = new ^TableQ; 

Pi^ 5??r.JE^;sri?SfiSo5Sr'^'SSi'!?nSkhoos=r(5ys«».9etProperty( 

"user. home", "" ) ); , 
private JMenultem jMenulteml = new JMenuitemQ, 
private jMenuitem jMenuitemZ = new DMenuitemQ; 

//construct the frame 
public ConverterFrameQ 

* enableEvents(AWTEvent.wiNDOW_EVENT_mSK) ; 
cry 

jbinitC); 

catch (Exception e) 

e.printstackTraceO ; 
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} ...... 

//Component imtianzation 

private void jbinitO throws Exception { 

imagel = new ^ ^-i 

imagelcon (edu . Stanford . f acs . 1 ogl i ke . ConverterFrame . cl ass . getResourceC openFi I e . gi t ) 

); 

image2 = new , t ^-t -^n 

imagelcon (edu . Stanford . f acs . 1 ogl i ke . ConverterFrame . cl ass . getResou rce ( cl oseFi 1 e . gi f 

)); 

images = new ^ _ .xnx> 

imagelcon (edu . Stanford . f acs . 1 ogl i ke . ConverterFrame . cl ass . getResou rce ( he! p . gi f ) ) ; 

//seticonimage(Tool kit .getDefaultToolkitO .createimage(ConverterFrame. class. getResou 

rce("[Your icon]"))); 

contentPane = (DPanel) this.getContentPaneO ; 
contentPane.setLayout(borderLayoutl) ; 
this.setsize(new Dimension(400, 300)); 
this.setTitle("Log Like Data Converter"); 
statusBar.setTextC' "); 
jMenuFi 1 e . setText ("Fi 1 e") ; 
iMenuFi 1 eExi t . setText ("Exi t") ; 

]MenuFileExit .addActionListener(new ActionLi stener() { 
public void actionPerformed(ActionEvent e) { 
j MenuFi 1 eExi t_acti onPerf ormed (e) ; 

})1 

iMenuHel p . setText ("Hel p") ; 
iMenuHel pAbout . setText ("About ") ; 

jMenuHel pAbout .addActionLi stener(new ActionLi stenerO { 
public void actionPerf ormed (ActionEvent e) { 
jMenuHelpAbout_actionPerformed(e) ; 

i Buttonl . seticon (i magel) ; 

3 Buttonl . addActi onLi stener (new j ava . awt . event . Acti onLi stener () 
public void actionPerf ormed (ActionEvent e) 
data„action(e) ; 

} 

}) ; 

jButtonl.setToolTipText("Open File") ; 
1 ButtonZ . seticon (i mage2) ; 

3 ButtonZ . addActi onLi stener (new j ava . awt . event .Acti onLi stener () 
public void actionPerf ormed (ActionEvent e) 
save_action(e) ; 

} 

j ButtonZ . setTool Ti pText ("cl ose Fi 1 e") ; 
1 Buttons . seticon (i mageS) ; 
3 Buttons . setTool Ti pText ("Hel p") ; 
menuSpectrum. setText ("Spectrum") ; 

menuSpect rum. addActi onLi stener(new Java . awt . event .ActionLi stener () { 
public void acti on Performed (ActionEvent e) { 
^ spectrum_action(e) ; 

iMenulteml. setText("Data") ; 

3Menulteml. addActionLi stener(new java.awt .event .ActionLi stener () 

public void actionPerf ormed (ActionEvent e) 
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{ . . ^ 

data_actionCe) ; 

})• 

iSiSSSllSz : IddiSi onLi s?en4r(new java . awt . event . Acti onLi stenerO 
public void acti onPerformed (Acti onEvent e) 
save_action(e); 

} 

iToolBar.add(jButtonl) ; 

iTool Bar.addQButtonZ) ; 

iTool Bar. add(j Buttons) ; 

iMenuFi 1 e . add(menuSpectrum) ; 

iMenuFi 1 e . add ( j Menuiteml) ; 

iMenuFi 1 e . add ( jMenultem2) ; 

iMenuFi 1 e . add (iMenuFi 1 eExi t) ; 

iMenuHelp.add(]MenuHelpAbout) ; 

]MenuBarl . addqMenuFi 1 e) ; 

jMenuBarl.add(]MenuHelp) ; 

thi s . setDMenuBarCjMenuBarl) ; 

contentPane . add( jTool Bar , BorderLayout . NORTH) ; 

contentPane . add(statusBar , BorderLayout . SOUTH); 

contentPane . add( j Scrol 1 Panel, BorderLayout .CENTER) ; 

jScrollPanel.getviewportO .add(spectrum_table, nui U ; 

//File I Exit action performed j,^- ^ ^«^J 

public void jMenuFileExit_actionPerformed(ActionEvent e) i 

System. exit(O); 

//Help I About action performed ^- ^ ^ j 

public void jMenuHelpAbout_acti onPerformed (Acti onEvent e) i 

//Overridden so we can exit when window is closed 
protected void processwindowEvent(WindowEvent e) { 

super. processwi ndowEvent(e) ; 

if (e.getlDQ == windowEvent .WIND0W_CL0SING) { 
jMenuFi 1 eExi t_acti onPerformed (nul 1 ) ; 

} ^ 

private static final String[] stringArray = new string[0]; 
private FCSFile fcs; 
private String [] sensor_names; 
private int[][] int_data; 

' " [] new„data; 
e[][] raw.data; 
e [] [] compensated_data ; 
e[] [] compensation_matrix; 



private int_ 
private doub 
private doub" 
private doub"" 



void compensate () 
throws lOException 

int n = fcs.getTotalC) ; 

int m = sensor_names. length; 

raw_data = new double [ m ][ n ]; 
for (int i = 0; i < m; ++i) 



^ FCSParameter p = fcs .getParameter( sensor„names[i] ); 
int[] input^data = int_data[p.getlndex()-lj ; 
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double [] output_data = raw_data[i]; 
if (p.isLogO) 

^ double[] antilog = new double[ p.getRangeQ ]; 
?5p?'o|aJean?nog^^ p.getMaximu^O ); 

^°output_data[j] ="'antiiog[ input_data[j] ]; 

} 

else 

^ double g = p.getGainO; 
for (int j = 0; j < n; 
output_data[j] = g * input_data[]] ; 

} 

} 

compensated_data = new double [ m ][ n ]; 
for (int i = 0; i < m; 

double [] mult = compensation_matnxli J ; 
double [] output_data = compensated_dataLi J , 
for (int j = 0; j < n; 

^ double X = 0; 
for (int k = 0; k < m; ++k)^, ^ 
X += raw_data[k][j] * mult[k]; 
output_data[j] = x; 

} 

} 

new_data = new int[ m ][ n ]; 
for (int i = 0; i < m; 

^ FCSParameter p = fcs.getParanieter( sensor_names[i] ); 
double[] input_data = compensated data [i] ; 
int[] transformed_data = new_data[i]; 
if (p.isLogO) 

^ double[] lookup = new double[ p.getRangeO ] ; 
lookup [0] = p.getMinimumO ; , , ,n n . 

interpolate( lookup, 0, p.getRangeQ, p.getMaximumQ ); 
for (int j = 0; j < n; ++j) 

^ int k = Arrays. binarysearch( lookup, input_data[j] ); 

if (k < 0) 
transformed_data[3] = -k; 

else ^.^ , 

transformed_data[]] = k; 

} 

} 

else 

^ double g = p.getGainO; 
^°?ransforjed_ditab] 2= (iJt)^ Math. floor ( input_data[j] / g ); 

} 

} 

} 

void spectrum_action ( 
ActionEvent e ) 
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int returnval = spectrumChooser.showOpenDialogC this ); 
ifCreturnval == DFileChooser.APPROVE_OPTiON) 
{ 

try 

^ BufferedReader br = new BufferedReaderC new FileReader( 
spectrumChooser.getSelectedFileQ ) ); 
br. readLineO ; 

br.readLineO; , . ^ u uA^^r^ ^- 

StringTokenizer st = new stnngTokemzerC br. readuneC) , \t J, 
ArrayList al = new ArravListQ ; 
while (st.hasMoreTokensQ) 

al.addC st .nextTokenO ); .... ^ • * ^ 

StringH columns = sensor^names = (Stnng[]) al . toArrayC stnngArray ) ; 
Double[][] rows = new Double[sensor_names. length + 1] [columns, length j , 
for (int r = 0; r < sensor_names. length; ++r) 

^ st = new stringTokenizerC br.readLineO, "\t" ); 
for (int c = 0; c < columns. length; ++c) 

^ rows[r][c] = new Double ( Double. par seDouble( st. nextTokenO ) ); 
} 

spectrunutable.setModeK new DefaultTableModel ( rows, columns ) ); 
catch (Exception ex) 
ex . pri ntStackTrace O ; 

} 

} 

compensation_matrix = new double[ sensor.names. length ][ sensor^names. length ]; 
for (int i = 0; i < sensor^names. length; ++i) 
compensation_matrix[i] [i] = 1.0; 

} 

private void interpolate ( 
doubled func, 
i nt i , 
int n, 
double X ) 

int m = n / 2; 

double y = Math.sqrt( func[i] * x ); 

func[ i + m ] = y; 

if (m > 1 ) 

interpolate( func, i, m, y ); 
interpolateC func, i + m, m, x ); 

} 

} 

void data_action(ActionEvent e) 

^ int returnval = spectrumchooser .showopenDialog( this ); 
if (returnval == JFilechooser.APPROVE_OPTlON) 
{ 

try 

^ FCSFile fcs = new FCSFile( spectrumchooser. getSelectedFileO ); 
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int n = fcs.getTotalO; 

1"^ ' = fcs.getParametersQ; 

]nly^n int value = new int[m][n]; 

intLJLJ _ ff-c aetinputlteratorO; 

'?r"(in't i = 0;' Imi '.SasSSfeEventsQ ; 

for fint j = 0; Imi .hasMorevaluesQ , ++D) 
^°int 5alue[j][i] = Imi .readvalueO ; 



this.int_data = int.value; 
this.fcs = fcs; 

catch (Exception ex) 

^ ex.printStackTraceO; 

} 

} 

} 

void save_action(ActionEvent e) 

{ 

try 



^ f?fl2Si 1 e( ' spectrumchooser . getsel ectedFi ^ e Q ) ; 
for'CiSt i = 0; i < sensor_names. length, 

^ PCSPara.eter oldp = fc|.getParajeterC -nsor^^ 
FCSParameter newp = f cs . addParameter w , 
newp . setBi ts( ol dp . getBi ts() ) . _ 
S.setRangeC oldp.getRangeQ ). 
if (oldp.isLogO) 

^ newD setoecadesC oldp.getDecadesQ ); 
SS'-lSSinimumC oldp.getMimmumQ ). 



else 



^ newp . setDecades ( Doubl e . NaN ) ; 
newp.setcainC oldp.getGannQ ), 

} 

int n = fcs.getTotalQ; 
iSHaSdle^^iSf^^fcrpJSSputiteratorO; 
for (int i = 0; Imi .RasMoreEventsQ ; ++0 

^ rini- i = 0* 1 < 1 nt_data . 1 ength ; 
%i.;rit^ValS4( lnt_data[j][i] ); 



Irm.closeO; 



} 
} 

} 

} 



catch (Exception ex) 
^ ex.printStackTraceO; 
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package edu . Stanford . f acs . 1 ogl i ke ; 

mport java.awt.*; 

mport ^ava.awt.event. ; 

mport Java.io.BufferedReader, 

mport lava.io.File; 

mport lava.io.FileReader; 

mport lava.io.iOException; 

mport lava.util .Arrays; 

mport lava.util .Arrayust; 

mport lava.util. Iterator; 

import 3 ava.util .Map; 

i mport 1 ava . uti 1 . Stn ngTokem zer ; 

import qavax. swing,*; 

import javax. swing. table.*, 



import org. i sac.*; 



/ 



* 



<p>TitTe: Log Like Data Transform</p> 
<p>Description: </p> -inn-y^/n-^ 
<D>copyright: copyright (c) 2002</p> 
SJSKpany: Stanford university</p> 



©author Wayne A. Moore 
* ©version 1.0 

*/ 



public class converterFrame extends DFrarae { 

=x&5Se EH^-pSoif i^J«e„uxte,0 . 

IrV^Vxl DButton 3Button2 = new DButtonQ 

SriJSl JButton 3Button3 = new JButtonQ; 

private imageicon imagel; 

private imageicon imageZ; 

orivate imageicon image3; -..^u^ir^. 

//construct the frame 
public converterFrameO 

^ enableEventsCAWTEvent.WlNDOW_EVENT_MASK) ; 
try 

^ jbinitO; 

catch (Except! on e) 

^ e.printStackxraceO; 
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//component initialization .^^ ^ 

nrivate void jbimti,; T:nruw:> t.yvv.^K «^^="^ 
Jai&^^^^^^ 
IgSU-anford.f^^^ 

//seticonl.age(ToolMt.get0efaultToolMtO.createI.age(ConverterFra.e.c 

'"^cirentPre^=^ ol>anel) this.getcontentPaneO ; 
coS^Spfne . setLayoutCb9rderLayoutl) , 
this setSizeCnew Dimension (400, i^^Jj',,^, 
^his slCTitleC'Log Like Data Converter ). 

statusBar . setText C , ; , , 
jMenuFile.setTextC Fije )- 

jMenuFi 1 eExi t_actionPerformed(e) . 



} 

}); 



jM4nuHe1p.setText("Help';); , 

jiJienuHel pAbout_acti onPerf ormed(e) , 

} 

feni:lfdiS 

public void actionPerformedCActionEvent e) 

^ data_action(e); 
} 

jBittonl.setTooTripTextropen File"); 

rBSliddicJiffilteJl^ 

^ public void actionPerformedCActionEvent e) 

^ save_action(e); 

} 

jBittonZ.setToolTipTextC'close File") ; 

3 Buttons. seticon (images); 

1 Buttons. setToolTipText( Help ), 

ienuspectrum. setTextrspectrum ) ;^ java. awt . event .ActionLi stener Q { 

"^Sfb??f Sdlctf^^^^^^^^^^ ' 
spectrum_action(e) ; 

} 

ISSlildd^^^^^ 

^ public void actionPerformed(ActionEvent e) 
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^ data_action(e); 
} 

^ public void actionPerfomedCActionEvent e) 
{ 

save_action(e); 

} 

f incti onchoi ce . setModeT (f uncti onModel ) ; 

iTool Bar . add(iButtonl) ; 

nToolBar.addQButtonZ); 

iToolBar.addClButtonS); . 

ItooI Bar . add (f uncti onchoi ce , null ) , 

nMenuFi 1 e . add(menuspectrum) ; 

iMenuFi 1 e . add( jMenuiteml) ; 

iMenuFi 1 e . add(:}Menuitem2) ; 

{SenuFi 1 e . add qMenuFi 1 eExi t) . 

iMinuHelp.addqMenuHelpAbout) ; 

i MenuBarl. add (]MenuFi 1 e) ; 
3 MenuBarl . add (: MenuHel p) ; 

thi s . set JMenuBar (DMenuBarl) ^ . . north) ; 

contentPane.add DToolBar |°^^3frLlyout . south) ; ^ 
contentPane . addCstatusBar , Boroer y CENTER) ; 

rsS^^rpfnii^ge^if^^^^^^^^^ -^^^ ' 

kVi^c iofd^5Mln;i??e?^;S?onPerfor.edCA« e) { 

system.exit(O) ; 

i„?5°f«e;rHMKrt?onPerfo™edCACtionEv.nt e) { 

jMln2FileEiit_actionPerformed(tiu11); 

private static final StringU strirgArray = new StringCO]; 

Dfivate FCSFile fcs; 
brivate stringG sensor^names; 
private int[][] int.data; 
^ . .-^nn new_aata; 



private int[ a^^^- 
nrivate double[][] raw^data, ^ . ^ 



raw^data; 
compensated^data; 
compensati on_mat n x ; 



); 

private double getModel Parameter ( 
int row, 
int col ) 
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Object obj = model_parameters.getvalueAt( row. col ). 

if (obj == null) 

return Double. NaN; 
if fobi instanceof Number; 

retuVn ((Number) obj) .doubleValueQ ; 
if (obj instanceof String) 
{ 

try 

^ return Double. par seDouble( (String) obj ); 

catch (NumberFormatException ex) 

^ return Double. NaN; 
} 

else 

^ model_parameters.setvalueAt( null, row, col ); 
return Double. NaN; 

} 

} 

orivate void invert ( 
doublenn matrix) 

^ int row = 0, col =.0. n = matrix. length; 

int pivotC] = new intLn]; 

int row_index[] = new intM, 

int col_index[] = new int[nj; 

for (int i = 0; i < n; ++i) 
' ?r(?nt'1 Si j<n; ..j) 

''%^'rV^-=l] k<n; ..k) 
^ if (pivot[k] == 0) 

^ double abs = Math.abs(matrix[j] [k]) ; 
if (abs >= big) 
{ 

big = abs; 
row 
col 

} 



row = ^; 



^SiritllgaliJgumentException( "matrix is singular" ); 

} 



++pivot[col. 
row_index[i. 
col_index[i" 



= row; 
= col ; 



if (row != col) 

^ system. out. println( "pivoting!" ); 
for (int k = 0; k < n; ++k) 

^ double t = matrix[row][k]; ^^^^ ^ 
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matrix [row] [k] = matrix [col] [k] ; 
matrix[col][k] = t; 

} 

} 

'^ff '£°9tegllrguSLtExcepn "matrix is singular" ); 

double inverse = l/matrix[col] [col] ; 

matrix[col] [col] =1; 

for (int j = 0; ] < n; ++3) 

matrix[col][j] *= inverse; 
for (int j = 0; j < n; ++3) 

if (j != col) 

^ double t = matrix[j] [col] ; 
matrix[j][col] =0; 
for (int k = 0; k < n; ++k) 

matrix[j][k] -= matnx[col] [k] * t, 

} 

} 

for (int i = n-1; i >= 0; .—)) . 
if (row_index[i] != col_index[i]) 
for (int j = 0; j < n; ++3) 

^double t = matrix[j][row_index[i]]; ^ 
matrix[i][row_indexti]] = matrix[3] [col_index[i]] , 
matrix[3][col_index[i]] = t; 

} 

} 

void compensate () 
throws lOException 

^ int n = fcs.getTotal () ; 
int m = sensor_names. length; 

double[][] compensation_matrix = new double[m] [m] ; 
for (int i.= 0; i .< m; 

compensation_matnx[i] [1] = 1.0. 
for (int i = 0; i < m; 

for (int j = 0; 3 < 

^ double coefficient = getModelParameter( i, j ); 
if (Double. isNaN( coefficient )) 

coefficient = 0; . 
compensation_matrix[3] [i] = coefficient, 

invert( compensation_matrix ); 

raw_data = new double [ m ][ n ]; 
for (int i = 0; i < m; ++i) 

^ FCSParameter p = fcs.getParameter( sensor names [i] ); 
int[] input_data = int_data[p.getindex()-l] ; 
double [] output_data = raw_data[i]; 
if (p.isLogO) 

^ double [] antilog = new double [ p.getRangeO ]; 
antiloq[0] = p.getMinimumO; ■ \. 

interp5late( antilog, 0, p.getRangeQ, p.getMaximumQ ); 
for (int j = 0; j < n; ^ . , 

output_data[j] = antilog[ input_data[3 J ], 
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} 

else 

^ double g = P-g«^^^^"^J,i 
for (int 3 = 0; 3.< 



} 

} 

} 

new_data =.new int[ m ][ n ]; 
for (int T = 0; i < m. 

rt^'"transf"o^SeS"a^a"Tl:Jata[i]; 

if (p.isLogO) 

^ doubleG lookup =.new doubleC p.getRangeQ ] ; 

lookupCO] = p.getMi"imumO . ^ p.getMaximumQ ) ; 

rublf'ra?Llret"^'getMSdllParameter( Ol 

if (! Double. isNaNC parameter )) 

^ switch (functionChoice.getselectedindexO) 

' "^'"^°k5Vh?!ogC^raJi:efe^ p.getMini.u.Q )) * (doubl 

.getRangeO -^f-cidesO^ fp??cl°?S!geUngeO - D 

^ double delta = Clookup[splice.l] - lookup [splice-l])/2-. 
Sbll blse = lookupLsp ice]; 
%rpC3r=2ise^-^s?Vi^cl'-l)^^* delta; 

} . 



break; 
case 1 



oarimeter *= p.getDecadesQ /,P-g«^'}^"?eO ; 
for (in? j =0: j < lookup.length. 
lookupCj] += parameter D, 



break; 

} 



for (int j = 0; j < n; 

Snt k = Arrays. binarysearch( lookup. inputJataUl ). 
^^transformed_data[j] = -k; ^^^^ g 
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^^transformed_data[j] = k; 

} 

else 



} 

} 



^ double g = p.getGainQ; 

'°?r^isforjed_dUj] 2' (irt)^ Math, round ( input_data[3] * 9 ); 

} 



void spectrum_action ( 
Act ion Event e ) 



i f (retJ rn5a! == 3 Fi 1 eChooser . APPR0VE_0PTI0N) 



{ 

try 



^ BufferedReader br = new BufferedReaderC new FileReaderC 
fileChooser.getSelectedFileO ) ); 
br.readLineO; 

S?ring?ok2S^^^^^ st = new stringTokenizer( br.readLineO. "Xt" ); 
ArrayList a1 = new ArravListQ; 
while (st.hasMoreTokensO) 

al.addC st.nextTokenQ ); fstrinqH) al.toArrayC stringArray 

fSrCint r = 0; r < sensor_names. length; ++r) 

^ St = new stringTokenizerC br.readLineO, "\t" ); 
for (int c = 0; c < columns. length, ++cj 

^ rows[r][c] = new Double( Double. par seDouble( st.nextTokenO ) ); 
} 

model parameters = new DefaultTableModel ( rows, columns ); 
spect7SSle.setModel( model .parameters ); 

catch (Exception ex) 

^ ex . pri ntStackTrace O ; 
} 

} 

} 

private void interpolate ( 
double □ func, 
int i , 
int n, 
double X ) 

{ _ / ?• 

dSSble V= Math.sqrtC func[i] * x ); 

func[ i + m ] = y; 
if (m > 1 ) 

^ Page 7 



parks_et_al_5.txt 

interpoTateC func, 1- yj'^ 
interpolateC func, i + m, m, x 

} 

} 

void data_action(ActionEvent e) 

if(retu?nval == jFilechooser .approvloption) 
{ 

try 

^ File f = fileChooser.getselectedFileO; 
FCSFile fcs = new FCSnieC r j; 
int n = fcs.getTotaU); 
i nt m = f cs . getParameters 0 , 
intnn int value = new intLmJLnJ. 
pStiidler Tmi = fcs.getinputlteratorQ ; 
for Cint i = 0; Imi .RasMoreEventsQ ; 
for Cint i = 0; Imi .hasMorevaluesQ ; ++:) 
in^l!)IluJ[j][i] = Imi.readvalueO; 
this.int_data = int_value; 

statiisBarTsetTextC f.getNameO ); 
catch (Exception ex) 

^ ex.printstackTraceO; 
} 

} 

} 

void save_action(ActionEvent e) 

^ fileChooser.setDialogTitleC "save New fcs Data File" ); 
Inr rpturnval = fileChooser.showsaveDialogC this 
i f (re?; rilSal == JFi 1 echooser . approve_OPTION) 
{ 

try 

fr'inf" Pi 0; i < sensor_names. length; ++i) 

^ Fcsparameter oldp = fcs.getParameter( sensor_naines[i] ); 

prIplSter newb = fcs.iddParameterQ ; r--, v 

newD sl?At?ribute° "$P", "N". "["+sensor_narnes[i]+ ] . ). . 

!!e!;p:lSA«ribute( "$P", "s", oldp.getAttribute( $P , S ) ). 

newp . setBi ts ( ol dp . getBi ts () ) ; 

newp.setRange( oldp.getRangeO ), 

newp.setLog( oldp.isLog() ); 

newp.setDecades( oldp.getDecadesQ ), 

newp. setMini mum ( oldp.getMirnmumO 

newp.setGain( oldp.getGainO 

} 

fcs setFile( fileChooser .getSelectedFileQ ); 
FCSHandler Imi = fcs.getOutputiteratorQ ; 
for (int i = 0; Imi .fiasMoreEventsQ ; 

^ for (int. j =0; j < i^^^?^^' ^^^^ 

Imi .writevalue( int_data[3J Li J )< 
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for (int 3=0; j < new.data. length; ++]) 
Imi .writevalueC new_data[]J [t J )\ 

Imi .closeQ ; 

catch (Exception ex) 

^ ex.printStackTraceO; 
} 

} 
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package edu . Stanford . f acs . 1 ogl i ke ; 

mport java.awt.*; 
mport java.awt .event.*; 
mport qava.io.BufferedReader; 
mport :java.io.File; 
mport lava.io.FileReader; 
mport :java.io.iOException; 
mport lava.util .Arrays; 
mport :java.util .ArrayList; 
mport qava.util .Iterator; 
mport java.util .Map; 
mport java.util .stringTokemzer; 
mport lavax. swing.*; 
mport ]avax.swing.table.*; 



import org. i sac.*; ,^ . t ■ ^ 

i mport javax . swi ng . event .Tabl eModel Li stener , 

/** 

* <p>Title: Log Like Data Transform</p> 

* <p>Description; </p> 

* <p>copyright: copyright (c) 2002</p> 

* <p>Company: Stanford Urn versity</p> 

* ©author Wayne A. Moore 

* ©version 1.0 
V 

public class converterFrame 
extends JFrame 

^ private JPanel contentPane; 

private JMenuBar jMenuBarl = new JMenuBaru; 

private JMenu jMenuFile = new 3Menu(); 

private DMenultem jMenuFileExi t = new jMenultemu; 

private JMenu jMenuHelp = new JMenuQ; 

private jMenultem jMenuHelpAbout = new DMenuItemc; , 

private JToolBar jToolBar = new DToolBarC); 

private 3Button iButtonl = new JButtonQ; 

private JButton ]Button2 = new JButtonQ; 

private DButton jButtonB = new JButtonQ; 

private imagelcon imagel; 

private imagelcon imageZ; 

private imagelcon image3; 

private DLabel statusBar = new DLabelQ; 

private BorderLayout borderLayoutl = new BorderLayoutC) , 

private 3ScrollPane jScrollPanel = new DScroMPaneO ; 

private 3Table spectrumTable = new JTableQ; 

private JMenultem menuSpectrum = new ^'^enuitemQ; . _^ nPtProoertvC 

private JFileChooser fileChooser = new JFileChooserC System. getProperxyi 

"user, home", ) );. ^ -.Monnri-omn - 

private 3Menuitem jMenulteml = new 3MenultemU, 
private 3Menultem jMenultemZ = new 3MenultemU; 

//construct the frame 
public ConverterFrameO 

^ enableEvents(AWTEvent.wlNDOW_EVENT_MASK) ; 
try 

spect rumTabl e . setModel ( model Parameters ); 
^ Page 1 
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catch (Exception e) 
e . pri ntStackTraceO ; 

} 

//component initialization 

private void jblnitO throws Exception { 

image;?oS(edu?s?anf ord . f acs J ogl i ke XonverterFrame . cl ass . getResource("ope^ 
); . 

imagS(edu"sLnford.facs.loglike.ConverterFrame.class.getResource("closeFile.gif" 

)); . 

imagS(edu"l?anford.facsJoglike.ConverterFrame. class. getResourceC-help.gif")); 

//seticonimageCToolkit.getDefaultToolkitQ.createimageCconverterFratne. class. getResou 

rce("[Your icon]"))); , . ^a. 

contentPane = (JPanel) thi s.getContentPaneQ , 
contentPane . setLayout(borderLayoutl) ; 
this.setsize(new Dimension (400, 300)); 
this.setTitleC'Logicle Data converter ); 
statusBar.setTextt" ") ; 
jMenuFi 1 e . setText ("Fi 1 e") ; 

jMenuFileExit.setText("Exit ); , -.^.^^^^r.^ t 

iMenuFileExit.addActionListener(new ActionustenerO I 
public void actionPerformed(ActionEvent e) I 
jMenuFileExit_actionPerformed(e); 

})i 

iMenuHel p . setText("Hel p") ; 

]MenuHelpAbout.setText( About ); , -^^^^^^r^ x 

lMenuHelpAbout.addActionListener(new ActionustenerC) i 
public void actionPerformed(ActionEvent e) { 
jMenuHelpAbout_actionPerformed(e); 

}); 

jButtonl.setlcon(imagel); «. a^*^««, ,-ei-onorn 

] Buttonl . addActi onLi stener (new ] ava . awt . event .Acti onu stener () 

public void actionPerformed(ActionEvent e) 

data_action(e) ; 

} 

1 Buttonl . setTool Ti pText ("open Fi 1 e") ; 

nButton2.setlcon(image2); . ^e^-^narn 

] Button2 . addActi onLi stener (new java . awt . event . Acti onLi stener (; 

public void acti onPerformed (Acti onEvent e) 
save_action(e) ; 

} 

3Button2.setToolTipText("Close File"); 

] Buttons. seticon (i mage3) ; 

] Buttons . setTool Ti pText ("Hel p ) ; 

raenuSpectrum.setlcon(new h^nn. 
imageicon(ConverterFrame. class. getResource( spectrum. jpg ))), 

menuSpectrum.setText("Spectrum") ; ^-^^i / 

menuspect rum . addActi onLi stener (new ]ava . awt . event .Acti onLi stener (J l 
public void acti on Performed (Acti onEvent e) { 
spect rum_acti on (e) ; 
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} 

jMenulteml.seticon(imagel) ; 

irnS"S:SIc"onS"i4rCne»:ava.a»t..vent.ActionLlst.rerO 
public void actionPerformed(ActionEvent e) 

^ data_action(e); 

} 

jMenuitem2.setlcon(image2) ; 

pui?r™2:lddIc"^nu?steUr(n.wjava.awt.eve«.ACtionti«enerO 

public void actionPerformed(ActionEvent e) 

^ save_action(e); 
} 

finctionChoice.setModel(functionModel); 
iButton4.setMaximumSize(new DimensionQbl, ^/jj, 
nButton4.setMinimumSize(new DimensionQbi, ^/)}, 

imaai?"c"n(2onvfrS^^ 

^i5S?!on4 addActionListenerCnew java.awt.event.ActionListenerO 

public void action Performed (Action Event e) 

^ spectrum_action(e) ; 

} 

]ToolBar.add(jButton4, null); 

]Tool Bar . add(] Buttonl) ; 

jTool Bar . add( j ButtonZ) ; 

iTool Bar. add(j Buttons); _ 

iTool Bar . addCf uncti onChoi ce , null ) ; 

iMenuFile.add(menuSpectrum) ; 

iMenuFile.addCjMenuiteml) ; 

iMenuFile.add(iMenultem2) ; 

i MenuFi 1 e . add (1 MenuFi 1 eExi t) ; 

iMenuHelp.addQMenuHelpAbout) ; 

]MenuBarl.add(jMenuFile); 

]MenuBarl.add(jMenuHelp); 

this,set3MenuBar(jMenuBarl) ; moo-tu^. 
contentPane . add ( jTool Bar . Bo rderLayout . NORTH) , 
contentPane.addCstatusBar, Bo rderLayout. SOUTH) , 
contentPane . add ( j Scrol 1 Panel , BorderLayout . center) , 
j?croll panel . getvi ewport () . add (spect rumTabl e , null); 

iof3lMl?ui??e^E^l!a???onPerfor.ed(ActionEvent e) { 
System. exit (0) ; 

bi? i<,id°j5enSM«!S???onPerfor..dCAct1onEve„t e) { 

//overridden so we can exit when window is closed 
protected void processWindowEvent(WindowEvent e) { 
suDer.DrocesswindowEvent(e) ; „,..^n r 

if (e getiDO == windowEvent.wiNDOW_CLOSlNG) { 
jwenuFi 1 eExi t_acti onPerf ormed (nul 1 ) ; 
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} 

vate static final String[] stringArray = new String[0]; 



pn 
pri 
pri 
pri 
pri 
pri 
pri 



vate FCSFile fcs; 

vate Strin^[] sensorNames = new StnngLOJ; 



vate int[. 
vate int[ 
vate doub 
^. .vate doub 
private doub 



originalData; 
□ new^data; 
e[][] raw_data; 
e [] [] compensated^data ; 
e[] [] compensation^matrix; 



private class Logi cleFunction 

^ public void transform ( 
int column, 

FCSParameter parameter ; 
^ if (parameter. isLogO) 

^ doubled lookup = new double[ parameter .getRangeQ ]; 

]S?^??^?iti(«?'6?Tara:S^^^^ parameter. getMaxi.u.O ); 

} 

} 

} 

private Logi cleFunction logLinearSplice = new Logi cleFunction Q 

public void transform ( 
int column, 

FCSParameter parameter ; 
^ if (Double. isNaNCmodelParameters.getParameterC sensorNames. length, column )) 
1 1 ! paramete r . i sLog Q ) 

^ super. transformC column, parameter ); 
return; 

} 

double[] lookup = new double[ parameter. getRangeQ ]; 

' if (Double. isNaNC parameterl )) 
• 1 1 break * 

(S5;;;?lo;c'p"Le?e?i^r-'«ath.loaC para.e«r.g«Mini«0 » ' double) 
para™eter.getRa„g.O/^para»e«r.g«De^ 

^ double delta = (lookup [splice+1] - lookup [splice-l])/2; 
double base = lookup[splice] ; 
for (int j = 0; j < splice; ++3) 

lookupLjJ = base - (splice - :) * delta, 

} 

} 

}; 

private class Model Parameters 
extends AbstractTableModel 
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private double[][] parameters = new double[OJ [J ; 

public int getRowCountQ 

^ return parameters. length; 
} 

public int getColumnCountQ 

^ return sensorNames. length; 
} 

public String getcol umnName ( 
int columnlndex ) 

^ return sensorNames [ columnlndex ]; 
} 

public Class getColumnClass ( 
int columnlndex ) 

^ return super .getColumnClass ( columnlndex ); 
} 

public boolean iscell Editable ( 
int rowindex, 
int columnlndex ) 

{ 

return true; 

} 

public Object getvalueAt ( 
int rowindex, 
int columnlndex ) 

^ double value = parameters[ rowindex ][ columnlndex ]; 

if (Double. isNaNC value )) 

return null; 

else n -1 ^ 

return stnng.valueOf ( value ;; 

} 

public void setvalueAt ( 
Object avalue, 
int rowindex, 
int columnlndex) 

^ double dvalue; 

if (avalue instanceof Number) 

dvalue = ((Number) avalue) .doublevalueO ; 
else if (avalue instanceof string) 
{ 

try 

^ dvalue = Double. valueof( (String) avalue ).doubleval 

catch (NumberFormatException ex) 

^ dvalue = Double. NaN; 
} 

} 

else 

dvalue = Double. NaN; 

Page 5 



parks_et_al_6.txt 
if (Double. isNaNC dvalue) && rowlndex < sensorNames. length) 
if ( rowlndex == colutnnindex) 

dvalue = 1; 
else 
dvalue = 0; 

parameters [ rowlndex ][ columnindex ] = dvalue; 
super. setvalueAtC avalue, rowlndex, rowlndex ); 

} 

public void setParameters ( 
double[][] parameters ) 

this. parameters = parameters; 
thi s . f i reTabl eStructureChangedO ; 

} 

public double getParameter ( 
int rowlndex, 
int columnindex ) 

^ return parameters [rowlndex] [columnindex] ; 
} 

private Model Parameters model Parameters = new Model Parameters Q ; 
private JComboBox functionchoice = new JComboBoxQ ; i„„ ^ ^ h" 

private. St ring[] functions = { "log linear splice", "x + a log x + b , 

^p?i late' Def ail tComboBoxModel functionModel = new Def aul tcomboBoxModel ( functi 

private double getModel Parameter ( 
int row, 
int col ) 

^ Object obj = model Parameter s.getvalueAtC row, col ); 
if (obj == null) 

return Double.NaN; 
if (obj instanceof Number) 

return ((Number) obj) .doubleValueO ; 
if (obj instanceof String) 
{ 

try 

^ return Double. parseDouble( (String) obj ); 
catch (NumberFormatException ex) 
return Double.NaN; 

} 

} 

else 

^ model Parameters. setvalueAt( null, row, col ); 
return Double.NaN; 

} 

} 

*/ 

private void invert ( 
double[][] matrix ) 

^ int row = 0, col = 0, n = matrix. length; 
int pivot[] = new int[n]; 
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int row_index[] = new int[n]; 
int co1_index[] = new int[n]; 

for (int i = 0; i < n; ++i) 

double big = 0; . 

for (int 3=0; ] < n; ++3) 

^ if (pivot[j] != 1) , 

for (int k = 0; k < n; ++k) 

^ if (pivot [k] == 0) 

^ double abs = Math.abs(matrix[j] [k]) ; 
if (abs >= big) 
{ 

big = abs; 
row = i ; 
col = k; 

^ } 

'"■throJ ^R'eHllegalirgumentException( "matrix is singular" 

} 

> ■ r n 

++pivotLcol] ; 
row_index[i] = row; 
col_index[i] = col ; 

if (row != col) 

^ system. out. println( "pivoting!" ); 
for (int k = 0; k < n; ++k) 

^ double t = matrix[row][k]; 
matrix[row][k] = matnx[col] [k] ; 
matrix[col][k] = t; 

} 

} 

'^^S'';;ew°l]legllA;guS^entException( "matrix is singular" ); 

double inverse = l/matrix[col] [col] ; 

matrix [col] [col] =1; 

for (int j = 0; 3 < ++3) 

matrix[col][j] *= inverse; 
for (int j = 0; j < n; ++3) 

if (j != col) 

^ double t = matrix[j] [col] ; 
matrix[j][col] =0; 
for (int k = 0; k < n; ++k) 
matrix[j][k] -= matnx[col] [k] * t; 

} 

} 

for (int i = n-1; i >= 0; —i) 
if (row_index[i] != col_index[i]) 
for (int 3=0; j < n; ++3) 

^ double t = matrix[i][row_index[i:i]; 
matrix [j][row_index[i]] = matnx[3] [col_index[i]] , 
matrix[3][col_index[i]] = t; 
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^ } 

nrivate static final double ln_10 = Math.logC 10 ); 
lAllll DsStton jButton4 = ne« JButtonQ; 

void compensate Q 
throws lOException 

Snt n = fcs-getTotalQ; 
int m = sensorNames. length, 

raw_data = new double[m] [n] ; 
for (int i = 0; i < m; 

^ i.r«;p;,rameter p = fcs.getParaiTieter( sensorNames [i] ) ; 
int[] inSSLdSta = or?ginalData p.getindexQ-l] . 
doubleC] output_data = raw_data[i], 
if (p.isLogO) 

^ doubleC] antilog = new doubleC p.getRangeQ ]; 

?S§fSeC ^an?frg?^r^^^etRangeC). p.getMaximu.Q ); 

'°Su^i"uLdita[jj Ja^t;iog[liput_data[3] ]; 
else 

^ double g = p.getGainQ; 
'%u^iSLdit=a[Si =^Sp"uLdrtatj] / g; 

} 

} 

doubl e [] [] compensati on_mat ri x = new doubl e [m] [m] ; 
for (int i = 0; i < m; 

'°[oSeJsajionVrix[ji[if=^"odelParameters.getPa i. J ); 

invert( compensati on_matrix ); 

compensated_data =.new doubl e[m] [n] ; 
for (int i = 0; i < m; 

^ doubleC] output_data = comperisated_data[i] ; 
for (int j = 0; j < m; ++3) 

' ffbll" J?'=VofplnsaS^^^^^^^ 

nu^ruLd'ata[Si i ?n"pit_ra^'a[k] * mult; 

} 

} 

new_data = new int[m][n]; 
for (int i = 0; i < m; 

^ Fcsparameter p = fcs.getParameter( sensorNames [i] ); 
doublecf inpu?-data =\ompensated data[i] ; 
int[] transformed_data = new_data[ij, 
if (p.isLogO) 

^ doubled lookup = new doubleC p.getRangeQ ]; 
lookup CO] = p.getMinimumO; ^^^^ ^ 
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inter 
switc 
{ 

case 0 



•polateC lookup, 0. P-get'^^"9^P„L P>?f ""^^"'"'"'"^^ ^' 
■W (f uncti onchoi ce . getsel ectedindexO) 



]ookup[0] = p.getMinimumO; p.getMaximumQ ); ^ 

if (Double.isNaNC parameter! 
break " 

C^^^^o;,%l'l^^^^^ p.getMini.u.O )) * (double) 

p.getRangeO ./^ P^.getD|cadesO /p^^-f ' .getRangeO - D 

^ double delta = (lookup [splice+1] - lookup[splice-l])/2; 
double base = lookup [splice]; 
for (int 1=0; j < splice, 
lookup[j] = base - (splice - ]) delta. 

} . 
break; 

case 1: ■ ■ 

]^t7c^ltl(^\ll^^^^^ p.getMaxi:n,u.O ); 

iSrIKe?! = mSdelParameters.aetParameter( m, i ). 
if (Double. isNaN( parameterl )) 

doSblf'parameterZ = model Parameters. getParameter( m + 1. i ); 
if (Double.! sNaN( parameters 

paraSIS^f t:^:r°a;eterl *,r/^^Sl^R^anreO ' ' '"''^^ 

lookSpCj] += parimeterl * j + parameterZ; 
break; 



case 2; 



?S?e??S^^^^^^ 0, p.getRangeO, Math.expC p . getoecades 0 * ln.10 ) 

Darameterl = model Parameters. getParameter( m, i ); 
if (Double.isNaNC parameterl }) 

parameterl = 0; . 
double X = p.getMimmumO; 



look 
break; 



lor (int j = 0; j < n; 

Snt k = Arrays. binarysearch( lookup. input_data[j] ); 

if (k < 0) . 
transformed_data[3] = -k; 



el se 

transformed_data[j] = k; 

} 

} 

else 



^ double g = p.getGainO; 

'^tr^nSforiedJUj] 2' (iJt)^ Math. round ( input_data[3] * g ); 
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} 

} 

void spectrum_action ( 
ActionEvent e ) 

^ fileChooser.setDialogTitleC "Load spectrum File" ); 
int returnval = fileChooser.showOpenDialogC this J, 
ifCreturnval == jFileChooser.APPROVE^OPTlON) 
{ 

try 

^ BufferedReader br = new BufferedReaderC new FileReaderC 
fileChooser.getSelectedFileC) ) ); 
br. readLineO ; 

String?o^^^^ st = new stringTokenizerC br.readLineQ , "Xf 

ArrayList al = new ArravListQ; 
while (st.hasMoreTokensQ) 

al.addC st .nextTokenQ ); . 
sensorNames = (StringH) al.toArrayC stnngArray ), 
double [][] parameters = new double [sensorNames. length + 
2] [sensorNames . 1 engtn] ; 

int r; n u ^ 

for (r = 0; r < sensorNames. length; ++r; 

^ st = new StringTokenizerC br. readLineO , "\t" ); 
for Cint c = 0; c < sensorNames. length; ++c; 

parameters [r][c] = Double. parseDoubleC st. nextTokenQ ); 

for (: r < parameters. length; ++r) 

for (int c = 0; c < sensorNames. length; ++cj 
parameters [r] [c] = Double. NaN; 
model Parameters. setParametersC parameters ;; 

catch (Exception ex) 
ex . pri ntStackTrace () ; 

} 

} 

} 

private void interpolate ( 
double [] func, 
int i , 
int n, 
double X ) 



{ 



intm=n/2; 

double y = Math.sqrt( func[i] * x j; 
func[ i + m ] = y; 
if (m > 1 ) 

interpolate( func, i, m, y ); 
interpolated func, i + m, m, x ); 

} 



} 

void data^acti on (ActionEvent e) 
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■ fileChooser.setDialogTitleC "Load Data File" ); 
int returnval = fileChooser.showOpenDialogC this ), 
if(re?umval == 3FileChooser.APPR0VE_0PTlON) 
{ 

try 

^ File f = fileChooser.getselectedFileO; 
FCSFile fcs = new FCSFileC f j; 
int n = fcs.getTotalQ; 
int m = fcs.getParametersQ; 
intnr] int_value = new int[ni]Lnj; 
FCSHandler Imi = fcs.getinputlteratorQ ; 
fS cint i = 0; Imi .RasMoreEventsQ ; ++0 
for Cint i = 0; Imi .hasMo revalues Q ; ++]) 
int_valui[j][i] = Imi .readValueQ ; 
this. original Data = int_value; 

statiisBarTsetTextC f.getNameO ); 

catch (Exception ex) 

^ ex.printstackTraceO ; 
} 

} 

} 

void save_action(ActionEvent e) 

^ filechooser setDialogTitleC "save New. fcs Data File" ); 
int re?SrnJal = fileChooser.showSaveDialogC this ); 
ifCreturnval == 3FileChooser.APPROVE_OPTlON) 
{ 

try 

compensate O ; 

• "iO ); 




^ FCSParameter oldp = fcs.getParameterC sensorNames[i] ); 
Fr<;paraffleter newp = nf cs .addParameterC) ; , 

newp.setBitsC oldp.getBitsQ ); 
newp.setRangeC oldp.getRangeQ 
newp.setLogC oldp.isLogO j; 
newp.setoecadesC oldp.getDecadesQ 
newp.setMinimumC oldp.getMinimumO 
newp.setGainC oldp.getGainQ j; 

} 

FCSHandler Imi = nfcs.getOutputlteratorO ; 
for (int i = 0; Imi .hasMoreEventsO ; 

^ for (int j = 0; j < original Data length; 
lmi.writevalue( originalDatah] [i] ) , 
for (int j = 0; ] < new_data. length, ++D) 
Imi .writevalue( new_dataLjJ Lt J -> . 



Imi .closeO ; 
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catch (Exception ex) 

^ ex . pri ntStackTraceO ; 
} 
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package edu . Stanford . f acs . 1 ogl i ke ; 

import java.awt.*; 

imoort 1 ava.awt. event. 

iKrt java.io.BufferedReader; 

import lava. 10. File; 

import lava.io.FileReader, 

import lava.io.iOException; 

import 5ava.util. Arrays; 

import lava.utiT .ArrayList, 

import iava.util. Iterator; 

imoort lava.util .Map; 

impoft java.util .stringTokemzer; 

import javax. swing.*: ^_ 

import javax.swing.taoie. , 

5 ja?ai!Ki*na • event .Tabl eModel Li stener ; 

^**<p>Title: Log Like Data Transform</p> 

* <p>Description: </p> onn^^/n^ 

* <D>coDvriqht: Copyright (c) 2002</p> 

* ^E>cSan?: Stanford university</p> 

* ©author Wayne A. Moore 

* ©version 1.0 
*/ 

public class converterFrame 
extends : Frame 

* ^ 3SL«Bfr!^. new :«.nusar0; 

nr-ivate JButton jButtonZ = new JButtonu. 

lAATe jButton3 = new JButtonQ; 

private imageicon imagel; 

private imageicon imageZ; 

orivate imageicon image3; -,,^u.^ir^- 

private J Label j Label 1 = new JLabelQ, 

//construct the frame 
public converterFrameO 

^ enableEvents(AWTEvent.WlNDOW_EVENT_MASK) ; 
try 
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jpS™lable.setModeK model Parameters ); 

catch (Except! on e) 

^ e.printStackTraceO ; 
} 

/«?Sd1irnftlf rh?Sws Exception ( 
imag»da?rtanford.facs.log1ikeXonverterPrame.class.^ openP.le.g.f ) 

lgS(edu"f?anford.facs.loglike.Conver« 
rmlgSU"s?anford.facs.loglikeXo^ 

//setxconxmageCToo1kit.getoefau1tToolkit0.createXmageCconverterPrame.c1ass.get.esou 

'""^cireStPa°n"e^=^ O^-anel) this.getContentPaneQ ; 
cSSlStptne . setLayout (b9rderLayoutl) 
this setSizeCnew Dimension (400, 300)), 
^hil.litTitleC'Logicle Data converter ), 
statusBar . setText ( ) ; 
jMenuFi1e.setText( Fije ) . 

?SMlll^U1ddISonLrs?en4r(new ActionListenerQ { 
^"•pSbl c vJid1ctionPerformed(ActionEvent e) { 
jMenuPileExit_actionPerformed(e) , 

jMenuHelp.setText("Help;;); 

rnreltfbrt:lfdIS^nU?fenlr(new ActionListenerQ { 
^ oSblic'^Coid actionPerformed(ActionEvent e) { 
jMeSuHelpAbout_actionPerformed(e) ; 

ps«rni:iddic?iJirf^^^^^ 

public void actionPerformed(ActionEvent e) 

^ data_action(e); 
} 

jBittonl.setToolTipTextC'open Pile") ; 

TsKzllfdiSffi^^^^^^^ 

^ public void actionPerformed(ActionEvent e) 

^ save_action(e) ; 
} 

jBittonZ.setToolTipText ("Close File") ; 

1 Buttons . setlcon(image3) ; 

] Buttons. setToolTipText( Help ), 

menuSpectrum.setText( Spectrum ), ^^^^ ^ 



spectrum_action(e) ; 



}); 



iM^nuiteml.seticonCirnagel); 

^ public void actionPerformedCActionEvent e) 
{ 

data_action(e) ; 

} 

jM4nultem2.setlcon(itnage2); 

^ public void actionPerformedCActionEvent e) 
^ save_action(e); 

^ public void actionPerformedCActionEvent e) 

^ functionChoiceCe) ; 

} 

iBitton4.setMaximumSizeCnew DimensionCSl, 27)); 
?Blltton4.setMinimumSizeCnew DimensionCSl. 27)J, 
jButton4.setlconCnew netResourceC'spectrum. jpg"))) ; „ 

^ public void actionPerformedCActionEvent e) 

^ spectrum_actionCe); 
} 

iLibell.setPreferredsizeCnew DimensionC17, 17)); 
iToolBar.addC3Button4, null); 

ItooI Bar . add Cd Buttonl) ; 

iToolBar.addCiButton2) ; 

iToolBar.addC^ButtonB); 

ijool Bar . add ClLabel 1 . ""1 ' .... 

ItooI Bar . addCf uncti onchoi ce , null). 

iMenuFile.addCmenuSpectrum) ; 

iMenuFi 1 e . addCiMenulteml) ; 

1 MenuFi 1 e . add Cl Menultem2) . 

iMenuFile.addseparatorC); 

1 MenuFi 1 e . add C j MenuFi 1 eExi t) ; 

5SHelp.addC3MenuHelpAbout); 

iMenuBarl.addgMenuFile); 

BMenuBarl.addCDMenuHelp) , 
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contentPane.addCstatusBar Borderuyout^S^^^^ 

kWl iofd'S«S;i??erxjS«?onPerfo™ed(Actio.Evert e) { 

System.exit(O); 

iofd°fMeSSlrpX;;S?onPerfor.ed(ActionEvent e) i 

jM^nSpi 1 eExi t_acti onPerf ormed (null); 

} 

. ^- ^ c4-r.-innn strinaArray = new String[0]; 
private int[]t] onginalData, 



private int. 
private doub 
private doub 
private dou 



Die 



newData ; 
][] rawData; 
'in compensatedData; . 
'] □ compensati onMatn x ; 



private class LogicleFunction 

^ public double □ transform ( 
int column, 

FCSParameter parameter ; 
^ doubleC] lookup = new double[ parameter. getRangeQ ]; 
S^^^^ parameter.getMaxi 



} 

> 



»lte(^?^^^^^^^^ parameter. getMaximumO ); 

return lookup; 



n , -.^^^r-cniirp - new LoqicleFunction Q 
private LogicleFunction logunearSplice - new Log 

public doubleD transform ( 
int column, 

FCSParameter parameter ; 
, . ...e,.— ^^^^^^^^ 

doubleG lookup = new double[ parameter. getRangeQ ]; 

cS5li'!o^('j-L^?-i^r"'^atK lo,C para^eter^getMini^O )) ' Rouble) 
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double delta = (lookup [splice.l] - lookup [spnce-l])/2. 
double base = lookupCsplice], 

^"^o^pt:] : 55^^ delta; 

} 

return lookup; 

J^ivate LogicleFunctionD transformer = { logLinearSplice }; 

private class Model Parameters 
extends AbstractTableModel 

^ private doubleGG parameters = new double[0][]; 
public int getRowCountO 
^ return parameters. length; 
} 

public int getColumnCountO 
^ return sensorNames. length; 
} 

public string getCol umnName ( 
int columnindex ) 

^ return sensorNamesC columnindex ]; 
} 

public Class getcolumnclass ( 
int columnindex ) 

^ return super. getColumnClassC columnindex ); 

} 

public boolean iscell Editable ( 
int rowlndex, 
int columnindex ) 

^ return true; 
} 

public Object getvalueAt ( 
int rowlndex, 
int columnindex ) 

^ double value = parameters [ rowlndex ][ columnindex ]; 
if (Double. isNaN( value )) 
return null; 

return string.valueof ( value ); 

} 

public void setvalueAt ( 
Object avalue, 
int rowlndex, 
int columnindex) 

^ double dvalue; 

if (avalue instanceof Number) ^^^^ ^ 
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dvalue = ((Number) avalue) .doublevalue() ; 
else if (avalue instanceof String) 
{ 

"Llue = Double.valueOf( (String) avalue ) .doubleValueO ; 

catch (NumberFormatException ex) 

^ dvalue = Double. NaN; 
} 

} 

TfelirNalt-dfaiue) «^ rowindex < sensorNames. length) 
if ( rowindex « columnindex) 

dvalue = 1; 
else 

dvalue = 0; 



} 

' ?!;^ i:frrira£?sit7u??«^4d() •. 
} 

public double getParameter ( 

int rowindex, 

int columnindex ) 
^ return parameters [rowindex] [columnindex] ; 
} 

^' -.^^^iDaraitiPters = new Model Parameters () ; , 

private Model Parameters model Parameter^ ^ ^^^..^ "x + a log x + b , 

^■private static final double ln_10 = Math.log( 10 ); 

nrivate void invert ( 

double [][] matrix ) 
^ int row = 0. col = 0 n = matrix.length; 
int Divot[] = new intLnJ. 
S ?ow_indexn = new int[n]; 
int col_index[] = new int[n], 

for (int i = 0; i < n; ++i) 
'fo;M = o: i<n; ..3) 
' ^%ir(fnt^^k''= S; k<n; -k) 

^ if (pivot[l<] == 0) p^g^ g 
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^ double abs = Math.abs(matrix[j] [k]) ; 
if (abs >= big) 

^ big = abs; 
row = i; 
col = k; 

} 



hy^r^ ^SiritlegalirgumentExceptionC "matrix is singular" ); 

} 



= row; 
= col; 



++pivot[col 
row_inclex[i 
col_index[i 

if (row != col) 

^ system. out. printlnC "pi voting!" ); 
for (int k = 0; k < n; ++k) 

^ double t = matrix[row][k]; 
matrix [row] [k] = matnx[col] [k] , 
raatrix[col][k] = t; 

} 

} 



^^^ror^K^^SluAE^^^^^^^ --trix is singular- ); 
double inverse = l/matnx[col] [col J , 
matrix [col] [col] =.1; . 
for (int D = 0; ] < +11^ 



or iinr ] = j -^ 

matrix[col][j] *= inverse; 

for (int j = 0; ] < n; ++]) 
if (j != col) 



} 
} 



^ double t = matrix[3][col]; 
matrix[j][col] = 0; 
for (int k = 0; k < n; ++k) ^ 
matrix[j][k] -= matrix[col] [k] t, 



for (int i = n-1; i >= Oj ."'t^T^ 
if (row index[i] != col_Tndex[i]) 
for (int j = 0; 3 < n; ++]) 



matrix[j][col_index[i]] = t, 
} 

} 

void compensate () 
throws lOException 



Snt n = fcs.getTotaK); 
int m = sensorNames, length, 

rawData = new double [m] [n] ; 
for (int i = 0; i < m; ++1) 
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f/-c nPtParameterC sensorNatnesLi J J. 
?S?[]"?n5u^ d'ata ^^rfJnll'oataCp^getindexO-l] ; 
double [] output_data = rawData[i]; 
if (p.isLogOj 

^ doubled antilog = new double[ p.getRangeQ ]; 

fnrelS^l^aleC ^an?!rog?^S^"^^'^etRangeO . p.getMaxi^u.Q ); 

'%u^j"uUtI[Si =\^t;iogrrnput_data[j] ]; 
else 

^ double g = p.getGainO; 
^u^pSUtlui ^^SpCt-drtlb] / g; 

} 

} 

double[][] compensationMatrix = new doubleW [m] ; 
for (int i = 0; i < m; . 

invertC compensationMatrix ); 

cotnpensatedData = new double[tn] [n] ; 
for (int i = 0; i < m; 

^ double □ output_data = cotnpensatedData [i ] ; 
for (int j = 0; j < m; ++]) 

' «l".u^rz'ofplnsaKf?i^ ; 
'%uii"uUta^i 4 ?n;it_raS[k] * .ult; 

} 

} 

newData = new int[m][n]; 

for (int i = 0; i < m; ++iJ 

^ FCSParameter p = fcs.getParameter( sensorNames[i] ); 
double!] inpu?-data =^compensatedData[i] ; 
int[] transformed_data = newData[i]; 
if (p.isLogO) 

^ doublen lookup = new double[ p.getRangeQ ]; 

swi tch (f uncti onchoi ce . getsel ectedindexO) 
{ 

case 0: • - r\. 

lookup[0] '.PiaetmniMO, p.getMaximumO ); ^ 

;;J??f'ra?L^°t?x''i L?i?rrS?s:g?tlarame«rC ™, 

if (Double.isNaNC parameterl 
b peak " 

^"^ (Sdi"og('iM?i^r-'^ath.log( p.getMini.u.O )) * (double) 
P.getRangeQ P(s??i?r>^"E {pli^e°<^^.getRangeQ - D 

^ double delta = (lookup [splice^l] - lookup[splice-l])/2; 
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); 



double base . lookup[|Pli«J.; 



} , 

break; 
case 1 



ra tSl?r= -2 "?ar°a;e?e?s oetParameterC m, . ). 
if (Double.! sNaNC parameterl 

do2Kf>ra.eter2 = .odel Parameters. getParameterC . . 1, ^ ); 
if (Double. isNaN(parameter2 )) 

P^:^:f !!:^oar?i.eterl * Math.log( p.getMinimutnO ) / ln_10, 
llrmTevl *= ? ge?DecadesO / P-g^^'^^"9eO ; 
parameter J. h y lookup, length ++V 

'"ioSSliJ paieterl '"j . paran,eter2; 



break; 
case 2 



1°?»(Mookup. 0. p.getRangeO. Math.exp( p.getoecadesQ * ln_10 ) 

narameterl = model Parameters. getParaineter( m. i ); 
?f (Souble.isNaN( parameterl )) 
parameterl = 0; . 

: ^x^' ilooSS5-^-T/^5ikurc3L - par-ten, 

break; 

for (int j = 0; j < n; 

Snt k = Arrays. binarysearch( lookup, input_data[3] ); 
^^transformed_data[j] = -k; 
*^transformed_data[j] = k; 

} 

} 

else 

^double g = p.getGain(); 

'°[r^isforied_dUj] ^ (i^t)' Math, round ( input.data[D] g ). 

} 



} 



void spectrum_action ( 
ActionEvent e ; 

' fileChoosensetoialogTitU 
?f'(refu"rral'==rFilShSoser.APPROVE_OPTiON) 

{ 

try , , , 

' BufferedBeader_, br , new sufferedReaderC ne« F.leKeaderC 
fileChooser.getselectedFileO ) ), 

br.readune(;; ^^^^ 9 
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br.readLineO; _ stringTokenizerC br.readLineO. "Xf 

al.addC st.nextToke^^^^^^ Il.toArrayC stringArray ); 
sensorNames = CStnngu.' T^„ni-h + 

..p^ - new double[sensorNames. length + 
doubleGG parameters - new aou 
21 [sensorNames. length J , 

fo^r (s''- 0- r < sensorNames. length; ++r) 

parameters [rjLcJ - 

} 

for (; .r < P^-'^^'c^^sSsSi-Names. length; ++c) 

'%a?amete'rs[r]tc] = Double.NaN; 
„odelParameters.setParameters( parameters ); 

catch (Exception ex) 
^ ex . pri ntstackTraceO ; 
} 



} 

} 



private void interpolate ( 
^ doubleD func. 



{ 



int 1 . 
i nt n , 
double X ) 

iSSble V=Vath.sqrtC funcUl • X V, 

funcL i + m ] = y. 
if (m > 1 ) 

} 



} 



void data_action(ActionEvent e) 

f • T /' "I r\if\ prs Data ft le j 9 

^ fileChooser.setDialogi^tleC Lo^^^^^^ ); 

rfCrSm^a? « 3F]fecrooler.APPROVE_OPT^^^^ 



{ 

try 



LI y 

< Pil^ f . fileChooser.getselectedFileOl 
^gFiil fcs . new fCSFileC f ). 

FCSHandler Imi = tcs.getJ.iw ^^^^ 
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this. original Data = int_value, 

SitutBlrTsetTktC f.getNameO ); 
catch (Exception ex) 
^ ex.printStackTraceO; 
} 

} 

} 

void save_action(ActionEvent e) 

' filechooser.setDialogTitleC ^S-e N|w F^^^^^^^^ 

{ 

try 

^ File f = fileChooser.getselectedFileO; 
compensate O; 

fo"onri'- 0;"l < sensopNames.length; 

{ ij -Pr-c nPt parameter ( sensorNames[i] ); 

FCSParameter oldp = J^=^-9|55J^S3erO ; 
FCSParameter newp - nrcs-au .,r.. nsorNames[i ]+„] .). , 
newp.setAttnbuteC $P , n„. o\dp ogtAttributeC "$P , S ) ). 
newD setAttnbuteC *p >. 
S : setii ts C ol dp . getBi t s 0 ) ; 

Slwp.setGainC oldp.getGainQ ). 

] .transformC i . newp ); 

for (int j =0. n = fcs.getTotalO; j < n; 
Snt k = Arrays.binarySearchC lookup, input_data[]] ); 
^^transformed_data[j] = -k-l". 
^^transformed_data[j] = k; 

} 



else 



^double g = "ewp-getcainO^ ^^i) 
^°ir^irforJeLSat"aH]'="?n^^^^^^^^^^^^ 



* g ); 
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} 

{ • n. -1 < oriainal Data. length; 

Imi .closeO ; 
latch (Exception ex) 
^ ex . pri ntstackTraceO ; 
} 

} 

} 

void functionChoiceCActionEvent e) 

i ^ r- ^ A n- i < sensorNames. length; 

for (int T = O'o.M < this.N PARAMETERS; ++]) 
for Ont 3 = 0, 3 < S.1N sensorNatne- 



} 



} 
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package edu . Stanford . f acs . 1 ogl i ke ; 

import java.awt.*; 

r/p^r? faralfo^BrffSedReader; 

import iava.io.iOException, 
import ]ava.util. Arrays; 
import nava.util.ArrayUst, 
import java.util. Iterator; 
-imoort lava.utT I .Map; . 
impSrt java.util .stringTokemzer; 
import :! avax. swing. ; 
import ] avax. swing. table. , 

import com. borl and. ibcl. layout. , 
import j avax. swing. border.*, 

^r<p>Title: Log Like Data Transform</p> 

I Sopy?ight:"wright CO 2002</p> 

* ^gJcoSpany: Stanford university</p> 

* ©author Wayne A. Moore 

* ©version 1.0 
V 

public class converterFrame 
extends J Frame 

' lZ^..?]^'ZTrV- new JMenuBarO; 

ShvI« JbUmor j Buttons - new JButtonO, 

private imageicon imagel; 
private imageicon imageZ; 

^ f.5;fJS"nl ie?J°Pr nditorPaneC) , 
??i«atc apanel aboutMessage = new ^PanelU. 
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Private GndLayout gn dLayout^^^^^^^ . 
I'AZTe com? onent co.ponent2 , . 
PpSrafe PI 6 = new 3LabeU); 

bi^areM^^ 

private Uabel 3^"=" _ ™ jpanelO; 

' enableEventsCAWVEvent.«KOO«-EVENTJ«SK) : 

try 

helpDialog.pacKO; 



catch (Except! on e) 
^ e.printstackTraceO; 
} 



fe"vS?d1i?i?U5"t'h?S«s exception < ,,„esourceC"open.il..g1f) 
.I,»u"rta„ford.facs.logme.converterPra™..cUss.,etK 

l,;frce3.-.?an.0...„^ 

images = new ^ , , ^^^e converterFraa,e.class.getResourceC"help.gif ))- 

rir;r.-r:-.«o...„^^^^ 

rce("[Your l^o"! -^j^^j.^nel) this.getContentPaneQ ; 
contentPane =„^^'^^"!'ieverticalstrut(8) ; 
componentl = |JJ-^'4|^|vfrticalStrut(8) ; 5^ 
component2 = createEmptyBorder(10,lb,iu,^3^ . 

borderl = border Factory. crea^^^M^ . 
components = ^ox-createve ^^^^ 



} 
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about_action(e); 



} 

}); 



^ public void actionPerformedCActionEvent e) 
{ 

data_action(e); 
|rtro"n:lft«^ - Pile") ; 

^ public void actionPerformedCActionEvent e) 

^ save_action(e) ; 
} 

ibSilMtS^^^^^^^ -ended PCS Pile") ; 

PtSiiraffi^^^^^ java.awt.event.ActionListenerO 
^ public void actionPerformedCActionEvent e) 

^ help_actionCe); 
} 

jBitton3.setToolTipTextC"Help"); 
imagSn?rore?ilfp'rre^c^^^ 

"""'mSS"sSe?trum . Text C spectrum ^ , event .Acti onLi stener Q { 

Tub?fc'-id -^"^^^^^^^^ ^ 
spectrum_actionCe) ; 

} 

ijinwenultem . setEnabl ed Cf al se) ; 
openMenultem.setlconOmagel). 

2re"nSS:l?dIS^ 

^ public void actionPerformedCActionEvent e) 

^ data_actionCe); 

} 

saieMenuitem. setEnabl edCfalse) ; 
saveMenuitem.setlconOmage2J, 

saveMenultem . setText C save as . . . ; , ^^^^^ ^^.^^ onLi stener Q 

saveMenultem . addActi onLi stener knew j a 

^ public void actionPerformedCActionEvent e) 

^ save_actionCe); 

JiLt1onCho1ce.setFontCre« java.awt.FontrDialo,". 1, 
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^ public void actionPerformedCACtionEvent e) 

^ function_action(e) ; 
} 

jBitton4.setMaximumSize(new Dimension (51, 27)); 

^ public void actionPerformedCACtionEvent e) 

^ spectrum_action(e); 
} 

iLibell setPreferredsizeCnew Dimension(17 , 17)) ; 
iplcfrlmSbl e . setRowsel ecti onAl 1 owed(f al se) ; 

KiigsiisiriddicSina^reni^ 

^ public void actionPerformedCACtionEvent e) 
^ help_dismissCe) ; 

|rn!;si:2:ifdis^^^^^^^ 

public void actionPerformedCACtionEvent e) 

^ help_actionCe); 
} 

heipPane. setPreferredsizeCnew DimensionC400. 400)); 

!!l'lS?rnrift?eiu"<Sffl^^^^^^ <head>\r\n\r\n </head>\r\n <body>\r\n <P>\r\n 
Logicle Data^^°"-ert;r^;^- </body>\r\n</html>\r\n") ; 

tfbl1l«"";?^o^??f°niuoresce„ce compensation", = 
ttSfll:t"Sir?;SS&'l9n£t?sSngconstants.cE^^^^ 

1 Label 3 . seticonTextGapCO) ; 
iiabel3 setTextC'wayne A. Moore ), 
iboS^iissIge . setLayout Cgri dLayoutl) ; 
gridLayoutl.setColumnsCl); p^^^ ^ 
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?[ibef4°setHoH zo:?il Ai i gn.ent (swi ngConstants . CENTER) ; 
3 Label 4 . setlconTextGap(O) ; „ 

1 Label 4. setTextC August 2002 .n-Q.. 2 10)); 

iLabelS.seticonTextGapCO); ^ f Trustees"); 

iLabelS-setTextC'Copynght 2002 by xne BOdru 

tlblr6;K^^eW?K^°-land Stanford .r. universnty"); 

] Label 7. setText ("David R. Parks ), 

statusBar.setTextC ); 

j panel 2 . setLayout (borderLayoutS) ; _ 

progressBar . setstri ng("compensate ) , 

lToolBar.add(jButton4, null); 

iTool Bar . add (openButton) ; 

iToolBar.add(saveButton) ; 

]ToolBar.addqButton3); 

ItooI Bar . add (1 Label 1 , nul 1 ) , 

iTool Bar . addCf uncti onChoi ce , nul 1 ) , 

1 Menu Fi 1 e . add (menuspect rum) ; 

iMenuFi 1 e . add (openMenultem) ; 

iMenuFi 1 e . add(saveMenuitem) ; 

iMenuFile.addseparatorC); 

1 MenuFi 1 e . addC j MenuFi 1 eExi t) , 

]MenuHelp.addC3Menuitem2) ; 

iMenuHel p . addseparatorO ; _ 
iMenuHelp.addCijMenuHelpAbout) , 

iMenuBarl.addqMenuFile) ; 

1 MenuBarl . add(3MenuHel p) ; 

ilJ?ollPa3e2.getviewportO .addChelpPane, null) , 

aboutMessage . add ( j Label 2 , nuT 1 ) 5 . . 
aboutMessage . add (cotnponentl 
aboutMessage. addQLabel 3, null) , 
aboutMessage. add (Uabel 7, nul i) , 

aboutMessage . add (3 Label 4 , nul 1 ) , 
aboutMessage.add(component2, nul 1) , 
aboutMessage. addCuabel 5, null) , 
aboutMessage . add (3 Label 6 , nul 1 ) , 
aboutMessage. add(component3, nul u . 

i„f3lMS:i?1e?xfflonPerfo™edCActionEvent O < 

system. exit (0) ; 
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private int 
private int._ 
private double 
private doub"*'" 
private doub 
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//Help I About action performed ^ . ^- ^ . 
//overridden so we can exit when window is closed 
protected void processWindowEventCwindowEvent e3 i 

super. processwindowEvent(e); ^. ^^^'K.^^ / 

if (e. get ID O == windowEvent. window closing) { 
j MenuFi 1 eExi t_acti onPerf ormed (nul 1 ) ; 

} ^ 

private static final String[] stringArray = new String[0]; 
private static final int N_PARAMETERS = 2; 
private FCSFile fcs; 

private String fcs_name; . 
private String[] sensorNames = new StringLOJ, 
^ • — - .--^nfi originalData; 

newData; 
rawData; 

compensatedData ; 
compensati onMat ri x ; 
private int progress; 

public class LogicleFunction 

^ public double [] transform ( 
int column, 

FCSParameter parameter ) 
^ doubled lookup = new double[ parameter .getRangeO ]; 
]S?^S?Lec''«'6r~TPnith, parameter.getMaximu™0 ) ; 

return lookup; 

} 

} 

public final LogicleFunction logLinearSplice = new LogicleFunction Q 

public double □ transform ( 
i nt col umn , 

FCSParameter parameter j 
^ double linear.range = model Parameters. getParameterC sensorNames. length. 

column ) ; . , n • ^^ 

if (Double. isNaNC linear.range )) 

return super .transform( column, parameter ;; , t- ^ 4. ?■ 

double In range = Math.logC parameter .getMaxi mum Q / tlnC 
Saramlte^TsSolcadesC Injknge / ln_10. parameter.getMaxnmumC) / Math.expC 

ln_range ) ); 

double[] lookup = new doubleC parameter .getRangeO ]; 
^S?Wa«c'S"'"6rffirT"nith, para.eter.get«axi-() ) , 

Mlfh!"gC l']SlL!?al!ge°rp^ran,«er getMinimumC) ) * Cdouble) 

Darameter.qetRangeQ / parameter. getoecadesQ / ln_iu J, 
parameter .ge^.^^y^w /^^^^^p^.^g < parameter .getRangeO - 1) 

^ double delta = (lookup [splice+1] - lookup [splice-l])/2; 
double base = lookup[splice] ; 
for (int i = 0; j < splice; ++d) 
lookupM] = base - (splice - j) * delta; 
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} 

return lookup; 



Luc final uo,Ule.unctio„ hyperbolic = new uo9lcU..nction O 

{ 



public doubleG transform ( 
int column, 



S°a"™l«rTS"«^^^^ para™eter.9etMaxi«0 / Math.expC 

ln_range ) ); 

^ M.n lookuD = new double[ parameter. getRangeO J, 

fSyrfet loo.p.Un,.K, pa™e»..e™a.«.0 / sc.e.fac«. 



}; 



} 



'°^„rpd : ?icoUc]r?/lorpb3ri'scaU.factor, 
return lookup; 



• T««Diii<: - new LoaicleFunction O 
public final LogicleFunction logPlus - new Logi 

^ public doubleG transform ( 
int column, . 

FCSParameter parameter ; T^„„i.h rolumn 

< ,o.ble para«eterl = ™odelPara^«rs.,e..ara.«erC 

double para«eter2 - «del.ara»eters.,et.ara.eterC «nsorsa.es.length X, 

"'"^'<Sr%ire?"^r^al!SSS^Jo|L, parameter ); 

if (Double. isNaNC parameterZ )) 

parameterZ = 0; 
doublen lookup =. new doublet parameter.getRanaeO 3; 
irS?slkc''?S5!^6rSr&^ para.eter.,et«axi-0 

return lookup; 



// 



};■ 



} 



private class Model Parameters ^^^^ ^ 
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extends AbstractTableModel 

^ private double[][] parameters = new double[0][]; 
public int getRowCountO 
^ return parameters. length; 
} 

public int getColumnCountQ 
^ return sensorNames. length; 
} 

public String getcol umnName ( 
int columnindex ) 

^ return sensorNames [ columnindex ]; 
} 

public class getcolumnclass ( 
int columnindex ; 

^ return super. getcolumnClassC columnindex ); 

} 

public boolean iscell Editable ( 
int rowlndex, 
int columnindex ) 

^ return true; 
} 

public Object getvalueAt ( 
int rowlndex, 
int columnindex ) 

^ double value = parameters: rowlndex ][ columnindex ]; 
if (Double. isNaN( value )) 
return null; 

^^Sturn string. valueofC value ); 

} 

public void setvalueAt ( 
Object avalue, 
int rowlndex, 
int columnindex) 

^ double dvalue; ^ . . 
if (avalue instanceof Number) 

dialue = ((Number) avalue) .doublevalueQ . 
else if (avalue instanceof String) 
{ 

try 

^ dvalue = Double.valueOf( (String) aValue ) .doublevalueQ ; 

catch (NumberFormatException ex) 

^ dvalue = Double. NaN; 
} 
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If'^CDJubllisS^^^^^ ^ rowindex < sensorNames. length) 

if ( rowindex == columnindex) 

dvalue = 1; 
else 

dvalue = 0; 

parametersL rowindex ][ columnindex ]= lvalue; 
fireTablecell updated ( rowindex, columnindex 

} 

public void setParameters ( 
double [][] parameters ) 

^ this. parameters = parameters; 
thi s . f i reTabl eStructurechangedO ; 

} 

public double getParameter ( 
int rowindex, 
int columnindex ) 

^ return parameters [rowindex] [columnindex] ; 
} 

}; 

private Model Parameters model Parameters = new Model ParametersQ ; 
private string[] functionNames = ^ ino x + b" >■ 

{ "log linear splice", "hyperbolic , x + a log x + D t, 
private LogicleFunction[] logicleFunctions = 

Sr]°vI^rSSUSmboBoxte'?;n^ctioni°o§r^ Ik. Def aul tcomboBoxModel ( 

^"Snire's^ItU final double In.lO = Math.logC 10 ); 
private BorderLayout borderLayout2 = new BorderLayoutQ , 

iT^rsTriir^l^^^^^ are not compatible with the 

spectrum." ; 

private void interpolate ( 
double [] func, 
int i , 
int n, 
double X ) 

intm = n/ 2; r-n* ^. 

double y = Math.sqrtC func[i] * x ), 

func[ i + m ] = y; 

if (m > 1 ) 

^ interpolateC func, i, m, y ); 
interpolateC func, i + m, m, x j; 

} 

} 

public void invert ( 
doubled [] matrix ) 

^ int row = 0, col = 0, n = matrix. length; 
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int pivot[] = new int[n]; 
int row_index[] = new intLnj; 
int col_index[] = new int[n]; 

for (int i = 0; i < n; ++i) 

double big = 0; 

for (int J = 0; 3 < n; ++]) 

* if (pivotCj] != 1) , 

for (int k = 0; k < n; ++k) 

^ if (pivot [k] == 0) 

^ double abs = Math.abs(matrix[j] [k]) ; 
if (abs >= big) 

big = abs; 
row = j ; 
col = k; 

^ } 

^^throw ^new°il legal Argument Except! on ( "matrix is singular" ); 

} 

> • r -n 
++pivot[col J ; 

row_index[i] = row; 

col_index[i] = col ; 

if (row != col) 

for (int k = 0; k < n; ++k) 

^ double t = matrixCrow] [k] ; 
matrix[row] [k] = matnx[col] [kj ; 
matrixCcol] [k] = t; 

} 

^^throw''^r5ew°illegalArgumentException( "matrix is singular" ); 

double inverse = l/matrix[col] [col] ; 

matrixCcol] [col] = 1; 

for (int j = 0; j < n; ++j) 

matrixEcol] [j] *= inverse; 
for (int j = 0; j < n; ++]) 

if (j != col) 

double t = matrixCj] [col] ; 
matrix[j][col] =0; 
for (int k = 0; k < n; ++k) ^ 
matrix[j][k] -= matnx[col] [k] - t; 

} 

} 

for (int i = n-1; i >= 0; ^--j) ^ 
if (row_index[i] != col_index[i]) 
for (int j = 0; j < n; ++3) 

^double t = matrix [i][row_index[i ]]; 
matrix[j][row_indexLi]] = matnx[3] [col_index[i]] , 
matrix [j][col_index[i]] = t; 

} 
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void compensate Q 
throws lOException 

int n = fcs.getTotalQ; 
int m = sensorNames. length; 

rawData = new double[m] [n] ; 
for (int i = 0; i < m; ++i) 

^ FCSParameter p = fcs. get Parameter ( sensorNames [i] ); 
int[] input_data = onginalData[p.getlndex()-lJ ; 
double [] output^data = rawData[i]; 
if (p.isLogO) 

^ double [] antilog = new double [ p.getRangeQ ]; 

antilogCO] = p.getMinimumQ ; . ^^ ^ 

interpolateC antilog, 0, p.getRangeO , p.getMaximumQ j; 
for (int j = 0; j < n; ++;)) 

output_data[j] = antilog[ input^dataLj] J; 

} 

else 

double g = p.getGainO; 
for (int j = 0; j < n; ^ , 

output_data[j] = input_data[]] / g; 

progressBar.setvalue( ++progress ); 

doubled [] compensationMatrix = new double[m] [m] ; 
for (int i = 0; i < m; ++i) 

for (int j = 0; i < m; ^ r • • ^ 

compensationMatrix[j] [i] = model Parameters. getParameter( i, ] j; 

invert ( compensationMatrix ); 

compensatedData = new double[m] [n] ; 
for (int i = 0; i < m; ++i) 

^ doubled output_data = compensatedData [i ] ; 
for (int j = 0; j < m; ++j) 

doubled input^data = rawDataCj] ; . r-nr-^ 
double coefficient = compensationMatnx[i] []] ; 
for (int k = 0; k < n; ++k) , 
output_data[k] += input_data[k] * coefficient; 

progressBar.setvalue( ++progress ); 

} * 

void spectrum_action ( 
ActionEvent e ) 

if (spectrumTable.isEditingO) . 
spect rumTabl e . getcel 1 Edi tor () . stopCel 1 Edi ti ng () ; 

fileChooser.setDialogTitle( "Load Spectrum File" ); 
int returnval = fileChooser.showOpenDialog( this ); 
if(returnval == JFileChooser . APPR0VE_0PTI0N) 
{ 

try 
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BufferedReader br = new BufferedReaderC new FileReader( 
fileChooser.getSelectedFileO ) ); 
br. readLineO ; 

sJring^^^^^ st = new stringTokenizerC br. readLineQ , "\t" 

ArrayList al = new ArravListQ; 
while (st.hasMoreTokensQ) 

al.addC st.nextTokenQ ); , • . ^ 

sensorNames = (StringG) al.toArrayC stnngArray ), 

double[][] parameters = new double [sensorNames. length + 
2] [sensorNames. length] ; 

int r; i. ^ 

for (r =0; r < sensorNames. length; ++r; 

^ st = new StringTokenizerC br . readLineQ , "\t" ); 
for fint c = 0; c < sensorNames. length; ++c; 

parameters [r][c] = Double. parseDoubleC st.nextTokenQ ); 

} 

for (; r < parameters. length; ++r) 

for (int c = 0; c < sensorNames. length; ++0 
parameters[r] [c] = Double. NaN; 

model Parameters. setParametersC parameters ); 
openButton.setEnabledC true ); 
openMenultem.setEnabledC true ); 

catch (Exception ex) 

ex.printStackTraceO ; 

} 

} 

} 

void data_actionCActionEvent e) 

^ if (spectrumTable.isEditingO) n^^-^-^^r^ . 

spect rumTabl e . getcel 1 Edi tor () . stopCel 1 Edi ti ng () , 

fileChooser.setDialogTitleC "Load PCS Data File" ); 
int returnval = fileChooser.showOpenDialogC this ;; 
ifCreturnval == JFileChooser.APPRQVE^OPTlON) 
{ 

try 

^ File f = fileChooser.getSelectedFileO; 
statusBar.setTextC f.getNameQ ); 
progress =0; _ 
progressBar.setMimmumC 0 ); 
progressBar.setvalueC 0 ); 

FCSFile fcs = new FCSFileC f ); ^ ^ 
for (int i = 0; i < sensorNames. length; 

if (fcs.getParameter( sensorNames [i] ) == null; 

^ DOptionPane.showMessageDialogC this, SPECTRUM_ERROR, 
"Error", JOptionPane. error^message ); 
return; 

} 

int n = fcs.getTotalO ; 

DroqressBar.setMaximum( n ); 
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progressBar.setEnabledC true ); 
int m = fcs.getParametersq; 
intnn int value = new intLmJLnJi 
FCSHandler Tmi = fcs.getinputiteratorO ; 
for (int i = 0; Imi .RasMoreEventsQ ; ++0 

^ for (int i = 0; Imi .hasMorevaluesO ; 
int_va1ue[j][i] = litii . readvalueQ ; 
progressBar.setvalueC ++progress 

this. original Data = int_value; 
this.fcs = fcs; 
fcs_naine = f .getNameO ; 

statusBar.5etText( fcs name ); 
progressBar.setvalue( 0 ); 
progressBar.setEnabled( false 
saveButton . setEnabl ed ( t rue ) ; 
saveMenuitem.setEnabledC true 

catch (Exception ex) 

^ ex.printStackTraceO ; 
} 

} 

} 

void save_action(ActionEvent e) 

fileChooser.setDialogTitle( "Save New. FCS Data File" ); 
int returnval = fileChooser.showSaveDialog( this ), 
if(returnval == 3FileChooser.APPROVE_OPTlON) 
{ 

try 

^ File f = fileChooser.qetSelectedFileO; 
Logi cl eFuncti on 1 ogi cT e = 1 ogi cl eFuncti ons [ 
functionChoice.getselectedindexO ] ; 

if (f.existsO) 

^ Object[] message = { f. "already exists. Do you want to replace it?" }; 
^^3rt^o"nPrnI:?KS??rm5ialog( this, message "overwrite" 

DO?tionPane.OK_CANCEL_OPTION. 30ptlOnPane.WARNING_MESSAGE )) 

return; 

} 

progress = 0; 

int m = sensorNames. length; 
progressBar.setMinimum( 0 J! ^ „ * „ ^. 
prog res sBar.setMaxi mum ( 2*m+m*m 
progressBar.setvalueC 0 ); 
progressBar.setEnabled( true ); 
statusBar.setTextC "working ); 

rawData = null ; 
compensatedData = null; 
newData = null ; 
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^ compensateQ ; 

catch cmegalArgumentException ex) 

^ DOptionPane.showMessageDialogC this. "Spectrum Matrix is singular 

' joptionpane.ERROR_MESSAGE ); 
return; 

rawData = null ; 

Fr5File newFcs = new FCSFile( fileChooser.getselectedFileQ ); 
newSatI = new int [sensorNames. length] [fcs.getTotal ()]; 
Ss get?extsegient0.readFrom( fcs.getTextSegmentO ), 
for (int i = 0; i < sensorNames. length; 

^ FCSParameter oldp = fcs.getParameter( sensorNames [i] ); 
if (oldp == null) 

^ joptionPane.showMessageDialogC this, spectrum_error. 
"Error", ] option Pane. error_MESSAGE ); 
return; 

FCSParameter newp = newFcs.addParameterO ; 
newp.setBits( oldp.getBitsO ); 

S;r«SL('o?d?"elSlc?deU), oldp.get«ini»mO ); 
newp.setGainC oldp.getGainQ ); 

double[] input^data = compensatedpata[i] ; 
int[] transformed_data = newData[ij; 
if (newp.isLogO) 

^ double [] lookup = logicle.transformC i, newp ); 
for (int j = 0, n = fcs .getTotal Q ; 3 < n; ++3) 

^ int k = Arrays. binarySearchC lookup, input_data[j] ); 

if (k < 0) ^ ^ rn 

transformed„data[]J = -k-i, 

else . . ^.^ , 

transformed_data[]] = k; 

} 

} 

else 

^ double g = newp.getGainO; . 

for fint i = 0, n = fcs.getTotal () ; 3 < ^^y^ r-i * « ^ 
tranSoriedJataCj] = (int) Math.round( input^data[]] g ) 

progressBar.setvalue( ++progress ); 

compensatedData = null; 

progress =0; ^ ^ ^ 

progressBar.setvalue( 0 ^ 
proqressBar.setMaximum( fcs.getTotal U 
statusBar.setText( f.getNameO ); 

FCSHandler Imi = newFcs.getoutputiteratorO; 
for (int i = 0; Imi .hasMoreEventsO ; 
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^ for (int j = 0; j < originalData. length; 

Imi.writevalueC o^iginalDatanKi] ). 
for (int j = 0; 3 < newData length, ++3^ 

Imi .writevalueC newData[3] Li J 
progressBar.setValueC ++progress 

Imi .closeQ ; 
newData = null; 

progressBar.setvalueC 0 ); 
progressBar.setEnabledC false ), 
statussar.setTextC fcs_name ;; 

catch (Exception ex) 

^ ex.printStackTraceO; 
} 

} 

} 

void function_action(ActionEvent e) 
' ^^^iSSfe^.gei^^^^^ 

^°modelParimiters.ie^valueAt( null. sensorNames. length ^3.1). 

} 

void help_action(ActionEvent e) 

^ helpDialog.setLocationRelativeTo( this ); 
helpDialog.setvisible( true j; 

} 

void help_dismiss(ActionEvent e) 

^ helpDialog.setvisible( false ); 
} 

void about_action(ActionEvent e) 

^ 30ptionPane.showMessageDialog( this, aboutMessage . "About Logicle", 
DOjtionPane.PLAiN_MESSAGE ); 

} 
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package edu . Stanford . f acs . 1 ogl i ke ; 

mport java.awt.*; 

mport lava.awt. event.*; 

mport qava.io.BufferedReader; 

mport iava.io.File; 

mport lava.io.FileReader; 

mport qava.io.lOException; 

mport java.util .Arrays; 

mport iava.util.ArrayUst; 

mport java.util .Iterator; 

mport java.util .Map; 

mport java.util .St ringTokemzer; 

mport lavax. swing.*; 

mport javax. swing. table.*; 

^Zorl i a3ai"wi ng . event .Tabl eModel Li stener ; 
i mport com . borl and . i bcl . 1 ayout . * ; 
import javax. swing. border.*; 

^* <p>Title: Log Like Data Transform</p> 

* <p>Description: </p> 

* <p>Copyright: copyright (c) 2002</p> 

* <p>Company: Stanford university</p> 

* ©author Wayne A. Moore 

* ©version 1.0 
*/ 

public class converterPrame 
extends 3 Frame 

^ private 3Panel contentPane; 

private JMenuBar jMenuBarl = new DMenuBarC;, 

private JMenu jMenuFile = new JMenuQ; 

brivate JMenuitem jMenuFileExit = new JMenultemQ, 

orivate JMenu jMenuHelp = new JMenuO; 

oAwlte JMenuitem jMenuHelpAbout = new JMenultemQ; 

brivate JToolBar jToolBar = new JToolBarQ; 

private JButton openButton = new JButtonU; 

private JButton saveButton = new JButtonC;; 

private JButton jButton3 = new JButtonQ; 

private Imagelcon imagel; 

private imagelcon imageZ; 

rri^a^e ^rdfrL^aVuHfr^erLayoutl = new BorderLayoutQ ; 
Srivate DScrollPane jscrollPanel = new J Scroll Pane 0, 
orivate DTable spectrutnTable = new JTableQ; 

"user. home", "" ) ); ..T^-^^r\^ 

orivate JMenuitem openMenultem = new DMenuitemQ , 
private JMenuitem saveMenuitem = new DMenultemQ , 
private JButton spectrumButton = "^w JButtonQ , 
private JComboBox functionchoice = new JComboBoxQ, 
private JLabel jLabell = new ^LabelC); 
private JDialog helpDialog = new JDialogO, 
private JPanel j Panel 1 = new jPaneK); 
orivate JButton helpDismiss = new JButtonQ;, 
pAwlll JScrollPane^ scroll Pane2 = new J Scroll PaneQ; 
pH^ati JEditorPane l^elpPane = new JEditorPaneQ; 
orivate JPanel aboutMessage = new JPaneKJ, 
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^ i-gL new uabeUV, 

private Border Dorderl, 

vate component component3 j^^^ . 

private 3Label ^Laoei/ "« ipaneU); 

//construct the frame 

public converterFrameO 

^ enableEvents(AWTEvent.WlNDOW_EVENT_MASK) ; 

try 

helpDialog.pacRQ; 



catch (Exception e) 
^ e.printStackTraceO; 



//component initialization^ Exception { 

private void ]bimtC) tnrows „PtResourceC"openFile.gif") 
x„agSU"s?anford.facs.loglike.conver«rFra™e.class.g«Msour.C 

;Urrce-."r«nfo...facs.1o.meXo„ve.e..a..cUss.,e«es„u.«C"c.o 

images = new i„nke.converterFra»e.class.getResourceC"help.gif ')) ; 

rceC" [Your icon] ^j): -,^ this.getcontentPaneO ; 

components = 5°J-^^||i|^|f^icalStrut(8) ; 
componentl = |°^-^^|g|virticalstrut(8) ; 
?o;^n^Pa^e:srt^a5^ut(borderLayoutl ; ^ 

r;i:iftri^!Mr^^ 

jMenuFile.setTextC File ) 

jwenuFileExit.setTextC Exit ActionListenerQ { 

'TuSril^vS^d-lSS^^^^^^^^^^^^ ^ 

P jMeSuMleExit_actionPerformed(e). 



} 
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'T„S^?»4ftMo™3c«^?onEvent .) ( 
about_action(e); 

} 

^ public void actionPerformedCActionEvent e) 
^ data_action(e); 

^ public void actionPerformedCActionEvent e) 

^ save_actionCe) ; 
} 

S^t^nllM^^^^ -tended PCS Pile"); 

PSSi:lfdAc?i^ju^^^^^^^^^ 

^ public void actionPerformedCActionEvent e) 

^ help_action(e); 
} 

5 Buttons . setToolTi pText C'Hel p") ; 
ipectrumMenultem.seticonCnew ^^^^^^ 

spectrurruaction Ce); 

} 

opinMenuitem. setEnabl edCfal se) ; 
openMenulteiti.setlconCirnagel). 

srerMfnKifdi^^^^^ 

^ public void actionPerformedCActionEvent e) 

^ data_actionCe); 
} 

saieMenultem. setEnabledCfalse) ; 
saveMenultem.setlcongmagez; , , 
saveWenultem.setTextC save A^ 
saveMenuItem . addActi onLi stener <.new ja 

^ public void actionPerformedCActionEvent e) 
^ save_actionCe) ; 
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Snctionchoice.setpreferredsizecn^ ^ Transformation ), 

^ public void actionPerformedCActionEvent e) 
^ function^actionCe); 

>itru™Button.set«axi«|i|eCn-D^ 
ir.«™Br«o"n;re™lTrpTe«^"Load Spectrum «atr,x ). 
SMCtrumButton.setlconCncm „|,gsource("spectrum. jpg"))) ; 

* public void actionPerformedCActionEvent e) 

^ spectrunuactionCe) ; 

^ public void actionPerformedCActionEvent e) 
^ help-disitiissCe); 

fei:ardlrt?an?enera4»iava.a»t.event.«tion.,^ 
^ public void actionPerformedCActionEvent e) 

{ 

help_action(e); 
Logicle Data converter^;^^. </body>\r\n</html>\r\n") ; 

Jrt£^12ll"e«rS1i?^?huorescence compensation") ; 
tfbfl'3:y™enW^^^^^ 

tfbllI-«tKnTS«lK^^^^^ . 
ttbelliemxtC-wayne A. Moore ), ^ 
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aboutMessage.setLayout(gridLayoutl); 
qridLayoutl.setColumnsUJ , 

^rafef4TetArMz°on?^li^l^gnment(S«in^ 
i Label 4 . setlcon-rextGapCO) ; _ 

1 Label 4. setText ("August 2002 ^' 2, 10)); 

jtlbelS.setFontCnew pva awt^F^^^^^^^ . 

tfbi^iferHrMforttnK 

iLabel5.setlconTextGap(0); ^ Trustees"); 

tfbfl|-|Stl°oMr«llTe;fpSirt.^lr«XingConstants.^ 



nill); 



1 Label 8 . setxexxi. vei =» • ""^x; 
3too1 Bar . addCspectrumButton , 
iTool Bar . add (openButton) ; 
iToolBar.add(saveButton) ; 
ijoolBar.addqButtonS); 
ixool Bar . add (iLabel 1 , nul 1 ) , . 
iTool Bar . addCtuncti onChoi ce , nul 1 ^ , 
lllmAVe . addCspectrumNienuitem) ; 
iMenuFi 1 e . add (openMenultem) ; 
3MenuFile.add(saveMenultem) ; 

iMenuFile.addseparatorC; ; 
^SuFile.addqMenuFileExit); 

iwenuHel p . add(]Menuitem2) , 

iMenuBarl. addC^MenuFi le) ; 
iMenuBarl.addCjMenuHelp) ; 

thi s . setDMenuBar (DMenuBarl) 5 ^ ^ . NORTH) ; ^ 

contentPane.addOToolBar, Boraer y ^^^^^^ 
contentPane . ^ddUScro 1 Panel . Border l y^^^^^ , 
contentPane • add ( jPanel 2 , Boraer y ^ . 

^PrnS^.td'dPp^rg^ SSaV Bor^^^^ ^ 

|^£°ora?rg^^efc5;!li:!rpS?6 S^cfprnS. BorderLayout.souTH); 

iPanell.addChelpDismiss null); 2. BorderLayout. CENTER) ; 

irifJaSergftSfeSSS^^^^^^ -^^^ ^ 

iboutMessage . idd C j Label 2. 
tbSStSfssage . add Ccomponent4 nu 1 ) ; 
aboutwessage . add p Label 8 , nul 1 ) , _ 
iboutMessage . add (componentl . nul 1 ) . 
aboutMessage . addQ Label 3 , nul 1) , 
aboutwessage . add Q Label 7 , nul 1 ) . 
aboutMessage . add (3 Label 4 nul 1 ) , _ 
aboutMessage . add (component2 null). 
aboutMessage. addOLabelS. null). 
aboutMessage. add (D Label 6, nuii;, 
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aboutMessage. add (components , null) ; 

km Uy5MSSi??eK;S?on«rfo™ed(ActionEvent e) I 

system.exit(O) ; 
//Help I About action performed ^^^^^^ 

jMln2MleEiit_actionPerformed(null). 

i iiii B ■ " 

^ sfr^ni!] slisorNames = new string[0]; 
private int[]t] origina Data, 
nr-ivate intf] [] newData, 
SnvSe doublet] □ rawData; ^ 
Kriva^ Sblin [] compensatedData; . 
Sn-vKl doSbllnn compensationwatrix; 
nrivate int progress;. ■^„„^r^a^ 
protected LogicleFunction logicle. 

public class LogicleFunction 

^ public double [] transform ( 
int column, , 
FCSParameter parameter ; 

^ doubled lookup = new double[ parameter. getRangeQ ] ; 

S»e«"6^«^^^^ parameter .getMaximumO ) ; 

return lookup; 

PU.11C final uogiclerunction log.ir.arsplic. = new ugUleFunction O 



public double □ transform ( 
int column. 



int column, - 

FCSParameter parameter ) 
t double linea..range - .odelPara.e«rs.,etPara™eterC sensorNa.es.length, 

Srrall«r:ScldSf iXfe'^r!"^^^ para..«r.g«Max«umO / Math.expC 
ln_range ) ) ; 

doubled lookup = new doubleC parameter .getRangeQ ] ; 
]rt»(«!^6?ToK^^ parameter.getMaximumO ); 

.ISl^.^SgC I^rn^e^ar:ifa5gr/"SU^^^ ^ * 
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narameter aetRangeO / parameter. getDecadesQ / in_10 ); 
parameter .getKanyeu /^"^^pi^^.g ^ parameter. getRangeO - D 

^ double delta = (lookup [splice+1] - lookup[splice-l])/2; 
double base = lookup[splice] ; 
for (int i = 0; j < splice; ++3) 
lookupCj] = base - (splice - 3) * delta. 

} 

return lookup; 

} 

}; 

public final LogicleFunction hyperbolic = new LogicleFunction () 
{ 

public double [] transform ( 
int column, 

FCSParameter parameter ) 

^ double linear_range = modelParameters.getParameter( sensorNames. length, 

column ) ; , ^ . 

if (Double.isNaNC linear_range )) ^ ^, 

return super .transform( column, P^/^meter ) , 
double scale factor = linear_range / (Math.E - l/Matn.E;, ^ . 

ln_range ) ); 



); 



double[] lookup = new double[ parameter. getRangeO ]; 

iS;e"r'pilLi(S^okfp'-o: lookup.length. parameter. getMaximumQ / scale_factor 



for Tint i = 0: i < lookup.length; 

lookSScj] = (ioolupCj] - l/loo^up[j]) * scale_factor; 

return lookup; 

} 

}; 

public final LogicleFunction logPlus = new LogicleFunction Q 

^ public double [] transform ( 
int column, 

FCSParameter parameter ) 
^ double parameterl = model Parameters. getParameter( sensorNames. length, column 
^' double parameter2 = model Parameters. getParameter( sensorNames. length + 1. 

column ); , ^ i n\ 

if (Double. isNaN( parameterl )) v 

return super .transform( column, parameter 
if (Double. isNaN( parameter2 )) 
parameter2 = 0; 

doubled lookup = new double[ parameter. getRangeO ]; 

]SS?lte(''Srp!^6n^^^^^ parameter.getMaximumO ); 

// oarameterZ += parameterl * Math.log( parameter. getMinimumO ) / ln_10; 

pa?ISfr! paftmeter.getDecadesO / parameter. getRange O ; 
for (int j = 0; j < lookup.length; ++D) 
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lookupCj] += parameterl * j + parameter2; 

return lookup; 

} 

}; 

private class Model Parameters 
extends AbstractTableModel 

private doubled [] parameters = new double[0][]; 

public int getRowCountQ 

return parameters . 1 ength ; 

public int getColumnCountO 
return sensorNames . 1 ength ; 

} 

public String getcol umnName C 
int columnlndex ) 

return sensorNames [ columnlndex ]; 

} 

public Class getcolumnclass C 
int columnlndex ) 

return super. getcol umncl ass ( columnlndex ); 

} 

public boolean iscell Editable ( 
int rowlndex, 
int columnlndex ) 

{ 

return true; 

} 

public Object getvalueAt ( 
int rowlndex, 
int columnlndex ) 

double value = parameters [ rowlndex ][ columnlndex ]; 
if (Double. isNaNC value )) 

return null; 
else 

return string.valueof ( value ); 

} 

public void setValueAt ( 
object avalue, 
int rowlndex, 
int columnlndex) 

double dvalue; 

if (avalue instanceof Number) 

dvalue = ((Number) avalue) .doublevalueO ; 
else if (avalue instanceof String) 
{ 

try 

Page 8 



parks_et_al_9.txt 
dvalue = Double.valueOf ( (String) avalue ) .doubleValueQ ; 

catch (NumberFormatException ex) 

^ dvalue = Double. NaN; 
} 

} 

else 

dvalue = Double. NaN; ^ n^,,«4-u^ 

if (Double. isNaNC dvalue) && rowlndex < sensorNames. length) 
if ( rowlndex == columnlndex) 

dvalue = 1; 
else 

dvalue = 0; 

parameters [ rowlndex ][ columnlndex ] = dvalue; 
fireTableCell updated ( rowlndex, columnlndex ); 

} 

public void setParameters ( 
doubled [] parameters ) 

this. parameters = parameters; 
thi s . f i reTabl eStructureChangedO ; 

} 

public double getParameter ( 
int rowlndex, 
int columnlndex ) 

^ return parameters [rowlndex] [columnlndex] ; 
} 

}; 

private Model Parameters model Parameters = new Model ParametersO ; 
private String[] functionNames = , 
{ "log linear splice", "hyperbolic", "x + a log x + b }; 
private LogicleFunction[] logicleFunctions = 

{ loqLinearSplice, hyperbolic, logPlus }; . ^ ^ i ^ 

private DefaultComboBoxModel functionModel = new DefaultComboBoxModeU 
functionNames ); . -, , > 

private static final double ln_10 = Math.loq( 10 ); 

private final String SPECTRUM_ERROR = "The data are not compatible with the 

spectrum."; 

private final int FCS^TIMER^TICK = 200; 

private void interpolate ( 
double [] func, 
int i , 
int n, 
double X ) 



{ 



intm=n/2; 

double y = Math.sqrt( func[i] * x ); 

func[ i + m ] = y; 

if (m > 1 ) 

interpolate( func, i, m, y ); 
interpolate( func, i + m, m, x ); 

} 
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} 

private void shuffle ( 
double[] sequence ) 

for (int i = 0; i < sequence. length; 

^ int i = i -f (int) Math.floor( Math.randomQ * (sequence. length - i) ); 
double t = sequence[i]; 
sequence[i] = sequence[j]; 
sequence[j] = t; 

} 

} 

public void invert ( 
doubled [] matrix ) 

int row =0, col =0, n = matrix. length; 

int pivot[] = new int[n]; 

int row_index[] = new int[n]; 

int col_index[] = new int[n]; 

for (int i = 0; i < n; ++i) 
{ 

double big = 0; 

for (int 3=0; j < n; 

if (pivotEj] != 1) 

for (int k = 0; k < n; ++k) 

if (pivot[k] == 0) 

double abs = Math.abs(matrix[j] [k]) ; 
if (abs >= big) 

big = abs; 
row = i ; 
col = k; 

^ } 

else if (pivot[k] > 1) 

throw new ll legal ArgumentException( "matrix is singular" ); 

} 

} 

++pivot[col] ; 
row_index[i] = row; 
col_index[i] = col ; 

if (row != col) 

for (int k = 0; k < n; ++k) 

double t = matrix[row] [k] ; 
matrix[row] [k] = matrix[co1] [k] ; 
matrix [col] [k] = t; 

} 

if (matrix [col] [col] == 0) 

throw new illegalArgumentException( "matrix is singular" ); 
double inverse = l/matrix[col] [col] ; 
matrix [col] [col] = 1; 
for (int j = 0; j < n; ++j) 

matrix[col] [j] *= inverse; 
for (int j = 0; j < n; ++j) 
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if (j != col) 

double t = matrix[j] [col] ; 
matrix[j][col] =0; 
for (int k = 0; k < n; ++k) 
matrix[j][k] -= matnx[col] [k] * t; 

} 



for (int i = n-1; i >= 0; --]}^^ 
if (row_index[i] != col_index[i]) 
for (int j = 0; j < n; ++]) 

^double t = niatrix[i][row_index[i]]; . 

matrix[i][row^indexti]] = matnx[j] [col^index[i J J ; 
matrixL]] [col_index[i]] = t; 

} 

} 

void make_busy () 

spectrumButton.setEnabled( false ); 
spectrumMenultem.setEnabled( false ); 

openButton.setEnabled( false ); 
openMenuItem.setEnabled( false ); 

saveButton.setEnabled( false ); 
saveButton.setEnabled( false ); 

progressBar.setEnabled( false ); 

} 

void make_ready () 

spectrumButton.setEnabled( true ); 
spectrumMenultem.setEnabled( true ); 

if (sensorNames != null) 

openButton.setEnabled( true ); 
openMenultem.setEnabled( true ); 

} 

if (fcs != null) 

saveButton.setEnabled( true ); 
saveButton.setEnabled( true ); 
statusBar.setText( fcs.getFileO .getNameO ); 

} 

else 

statusBar.setText( ); 

progressBar.setValue( 0 ); 
progressBar.setEnabled( false ); 

void spectrum^action ( 
ActionEvent e ) 

if (spectrumTable.isEditingO) , 
spect rumTabl e . getcel 1 Edi tor () . stopCel 1 Edi ti ng () ; 
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fileChooser.setDialogTitleC "Load Spectrum File )] 
int returnval = filechooser.showOpenDialogC this ); 
if (returnval == jFil eChooser. appR0VE_0PTI0N) 

^ make_busy(); 
try 

^ BufferedReader br = new BufferedReaderC new FileReaderC 
fileChooser.getselectedFileO ) ); 
br. readLineO ; 

's^ring?oi"en?z4r st = new stringTokenizerC br. readLineO . "\t" ); 
ArrayList al = new ArravListQ; 
while (st.hasMoreTokensO) 

seSo?Ser;"^^?rin^[]S ^i.toArrayC stringArray ); 
double[]n parameters = new double[sensorNames. length + 
2] [sensorNames . 1 ength] ; 

for (r'= 0; r < sensorNames. length; ++r) 

^ St = new StringTokenizerC br. readLineO , "\t" ); 
fnr rint c = 0; c < sensorNames . 1 ength ; ++c; 

paraSete^sCrKc] = Double. parseDouBleC st.nextTokenO ); 

} 

for C; r < parameters. length; ++r) 
for Cint c = 0; c < sensorNames .length; ++c) 
parameters [r][c] = Double. NaN; 
model Parameters. setParametersC parameters 

if Cfcs != null) ^ . . 

for Cr = 0; r < sensorNames. length; ++r; 
if Cfcs.getParameterC sensorNames [r] ) == null; 

fcs = null; 
break; 

} 

catch CException ex) 
^ ex.printStackTraceO; 
make_readyO ; 

} 

} 

FCS File fcs_data; 
FCSHandler fcs_handler; 
lOException fcs_exception; 

Timer fcs.timer = new TimerC fcs_timer_tick , new ActionListener O 

^ public void actionPerformed C 
ActionEvent action ) 

^ progressBar.setvalueC fcs_handler.getEventO ); 

} 

} ); 

void data_actionCActionEvent e) 
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if (spectrumTable.isEditingO) - 
spect rumTabl e . getcel 1 Edi tor () . stopCel 1 Edi ti ng () , 

fileChooser.setDialogTitleC "Load FCS Data File" ); 
int returnval = fileChooser.showOpenDialogC this ), 
ifCreturnval == 3FileChooser.APPROVE_OPTiON) 
{ 

try 

^ File f = fileChooser.getSelectedFileO; 

fcs^data = new FCSFileC f ); 

for (int i = 0; i < sensorNames.length; 

if (fcs^data.getParameterC sensorNames[i] ) == nulU 

^ jOptionPane.showMessageDialogC this, spectrum_error, 
"Error", 30ptionPane.ERROR.MESSAGE 
return; 

} 

make_busyO; . . 
progressBar.setmmmumC 0 j; 
progressBar.setvalueC 0 ); n 
progressBar.setMaximumC fcs^data.getTotal U J; 
progressBar.setEnabledC true ); .... 
statusBar.setTextC fcs^data.getFileQ .getNameQ J, 

new ThreadC data^thread, "FCS reader" ).startO; 
f cs_timer. start (); 

catch (Exception ex) 

ex . pri ntStackTraceO ; 
system. exit ( 1 ); 

} ^ 

} 

Runnable data^thread = new RunnableO 
public void run () 

^ Thread. currentThread().setPriority( Thread.MIN_PRIORITY ); 
fcs_exception = null; 
int n = fcs_data.getTotal(); 
int m = f cs_data. get Parameters () ; 
int[][] int_value = new int[m][n]; 
try 

^ fcs handler = fcs^data.getinputiteratorO ; 

for'(int i = 0; fcs_handler .hasNloreEventsO ; ++0 
for (int i = 0; fcs„handler.hasMorevalues() ; ++]) 
int_value[j] [i] = fcs_handler . readvalueO ; 
fcs_handler.close() ; 
original Data = int_value; 

catch (lOException ex) 

^ ex. pri ntStackTraceO; 
fcs.exception = ex; 

SwingUtilities.invokeLater( data_compl ete ); 
^ Page 13 
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}; 

Runnable data_compl ete = new Runnable 0 
^ public void run Q 

^ fcs_timer.stopO; 
if (fcs_exception == nullj 

^ fcs = fcs_data; 

make_readyO; 

} 

}; 

void save_action(ActionEvent e) 

^ if (spectrumTable.isEditingO) TTc^-i+n-nnr^ • 

spect rumTabl e . getcel 1 Edi tor () . stopCel 1 Edi ti ng () , 

fileChooser.setDialoQTitleC "save New FCS Data File" ); 
int returnval = filechooser.showSaveDialogC this ), 
ifCreturnVal == 3FileChooser.APPROVE_OPTlON) 
{ 

try 

^ File f = fileChooser.getselectedFileO ; ^^A-r^A^^r^ i- 

ligicle = logicleFunctTonsC functionChoice.getSelectedindexQ ]. 

if (f. exists O) 

^ Object[] message = { f. "already exists. Do you want to replace it?" 

^fogrio"nSSI:sfero;??rrDialog( this. StRNirMTSicE')) 
30ptionPane.0iC.CANCEL_0PTI0N, 30ptionPane.WARNING.MESSAGE )) 

return; 

} 

int n = f cs . getTotal O ; 
int m = sensorNames. length; 
rawData = null ; 
compensatedData = null; 
newData = null ; 

compensationMatrix = new double[m] [m] ; 
for (int i = 0; i < m; 

^'^o^p2Ssa'tionSitr?x[j]ti]^^^ ^ ); 

try 

^ invertC compensationMatrix ); 
catch (illegalArgumentException ex) 

^ 30ptionPane.showMessageDialogC this, "spectrum Matrix is singular", 



"Error' 



joptionPane.ERROiOiESSAGE ); 
return; 

} 

make^busyO; 

progress = 0; „ i>i 
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progressBar.setMinimumC 0 ); * „ 
progressBar.setMaximumC 2 * m + m m j, 
progressBar.setvalueC 0 ); 
progressBar.setEnabledC true ); 
statusBar.setxextC "working ); 

fr«; data = new FCSFileC fileChooser.getSelectedFileQ ); 
fcfcdSLgeJfex^igienio.readFromC fcs.getTextSegmentQ ); 

new ThreadC xf rm_thread, "Transform" ).startO; 

catch (Exception ex) 

^ ex.printStackTraceO; 
} 

} 

} 

Runnable xfrni_thread = new RunnableC) 
public void run () 
try 

int n = f cs , getTotal 0 ; 
int m = sensorNames. length; 

rawData = new double[m] [n] ; 
for (int i = 0; i < m; 

^ FCSParameter p = fcs.getParameter( sensorNames [i] ); 
int[] input_data = onginalData[p.getindexO-lJ , 
double [] output_data = rawData[i] ; 
if (p.isLogO) 

^ double[] antilog = new double[ p.getRangeO ]; 

antiloqro] = p.getMinimumO ; • r\ \. 

iSter?Slate( antilog. 0. p.aetRangeQ. p.get^^^^ ); 
doublet] fuzz = new doublet 1«FUZZY_BITS ], 

intetSilatei: fuzz. 0, 1«FUZZY_BITS , anti log[l] /anti logtO] ); 
shuffle( fuzz ); . 

^°Su^pSLditatji ='aStiiogrrnput_datatj] ] * fuzz[ j & 

(1«FUZZY_BITS-1) ]; 

else 

^ double g = p.getGainO; 

for (int j = 0; j.< n; ++]) . , 
output_datatj] = input_data[]] / g; 

} 

++progress; , j ^. 

SwingUtilities.invokeLater( progress_made ), 

} 

compensatedData = new doublefm] tn] ; 
for (int i = 0; i < m; 

^ doublet] output_data = compensatedDatati] ; 
for (int j = 0; j < m; ++]) 

^ doublet] input_data = rawDatatj]; 

Page 15 
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double coefficient = compensationMatnx[i] UJ ; 
for Cint k = 0; k < n; ++k) . 
output_data[k] += input_data[k] * coefficient; 

Swingutilities.invokeLaterC progressjnade ); 

} 

} 

newData = new int[m][n]; 

for (int 1=0; i < m; ++i) 

^ FCSParameter oldp = fcs.getParameter( sensorNames[i] ); 
FCSParameter newp = fcs_data.addParameterU ; , 
newp setAttribute? "$P", "N';. \["+sensorNamesCi ]+„] ); , 
newp.setAttributeC "$?", "s", oldp.getAttnbuteC $P . S ) ), 
newp.setBitsC oldp.getBitsC) ); 

"nSSJSSdLc'o?d?"e?Slc^L'sb. oldp.getMinin,u.O ); 
newp.setGainC oldp.getGainC) ); 

double[] input_data = compensatedData[i] ; 
int[] transformed_data = newData[i]; 
if Cnewp.isLogO) 

^ double[] lookup = logicle.transformC i, newp ); 
for Cint j = 0; j < n; ++d) 

^ int k = Arrays. binarySearchC lookup, input_data[j] ); 

if Ck < 0) ^ ^ ^.^ 
transformed_data[3] = -k-1; 

else , 
transformed_data[3J = k; 

} 

} 

else 

double g = newp.getGainC) ; 

^°?ran2for3ied_diita{j] ='Ci?it) Math.roundC input_data[j] * g ); 

} 

++progress; . , ^ 

Swingutilities.invokeLaterC progress_made ), 

} 

compensatedData = null; t ^. 

Swingutilities.invokeLaterC xfrm_coraplete 

catch CException ex) 

^ ex.printStackTraceC); 

} 

} 

}; 

Runnable progress_made = new Runnable C) 
^ public void run C) 

^ progressBar.setvalueC progress ); 
} 

}; 

Runnable xf rm_complete = new Runnable C) 
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parks_et_al_9.txt 

^ public void run O 

^ progressBar.setMinimumC 0 ); 

progressBar.setvalueC 0 ); ^^^T^+ain -4- 

progressBar.setMaximumC fcs_data.getTotal Q 
progressear.setEnabledC true ); „„^.K,,„on ^■ 

?tatusBar.setText( fcs.data.getFileO .getNameQ ), 

new ThreadC save_thread, "FCS writer" ).startO; 
fcs_timer. start O; 

} 

}; 

Runnable save_thread = new RunnableO 

^ public void run () 

^ Thread. currentThreadO.setPriorityC Thread.MlN_PRlORlTY ); 
fcs_exception = null; 
try 

^ fcs handler = fcs.data.getOutputiteratorQ ; 
for Cint i = 0; fcs_handler.hasMoreEventsO ; 

^ for Cint j = 0; j < originalData.length; ++j) 
fcs_hand1er.writevalue( onginaloataL:] Li J J. 
for (int j = 0; j < newData. length; ++3; 
fcs_hand1er.writevalue( newDataCj] [i J ), 

^cs_handler.close() ; 

catch (lOException ex) 

^ ex.printStackTraceO ; 
fcs_exception = ex; 

SwingUtilities.invokeLater( save_compl ete ); 

} 

}; 

Runnable save_compl ete = new Runnable () 
^ public void run () 

fcs_timer.stop() ; 

if (fcs_exception == null) 

{ 

> 

make_ready() ; 

} 

private 3 Label j Label 8 = new J Label (); 
private component componentl; 
private component component4; 

void function_action(ActionEvent e) 

^ if (spectrumTable.isEditingO) 

spectrumTabl e . getcel 1 Edi tor () . stopCel 1 Edi ti ng () , 

for (int i = 0; i < sensorNames. length; ++i) 

for (int 1=0; j < thi s . n_parameters ; ++]) 
Tor urn. J v, J 
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model parameters. setvalueAtC null, sensorNames. length + ], i 

} 

void help_actionCActionEvent e) 

^ helpDialog.setLocationRelativeToC this ); 
helpDialog.setvisibleC true ); 

} 

void help_dismiss(ActionEvent e) 

^ helpDialog.setvisibleC false ); 
} 

void about_action(ActionEvent e) 

^ DOptionPane.showMessageDialogC this, aboutMessage, "About Logicle" 
JOpti onPane . PLAIN_MESSAGE ) ; 
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package edu . Stanford . f acs . 1 ogl i ke ; 

mport java.awt.*; 
mport ] ava. awt. event. ; 
mport ^ava.io.BufferedReader, 
mport lava.io.File; 
mport ]ava.io.FileReader; 
mport java.io.IOException; 
mport java.util -Arrays; 
mport java.util .ArrayLi St; 
mport java.util .Iterator; 
mport ]ava.util .Map; 
mport java.util .Random; 
mport qava.util .StringTokemzer; 
mport lavax. swing.*; 
mport javax.swing.table.^; 

mport org. i sac.*; 
import javax. swing. event.*; 
i mpo r t com . bo r 1 and . i bcl . 1 ayout . * ; 
import javax. swing. border.*; 

* <p>Title: Log Like Data Transform</p> 

* <p>Description: </p> ^nn^^/n^ 

* <p>Copyright: Copyright. (c) 2002</p> 

* <p>Company: Stanford umversity</p> 

* ©author Wayne A. Moore 

* ©version 1.0 
*/ 

Dublic class converterFrame . , • ^ 

extends ]Frame implements TableModel Listener 

^ orivate J Panel contentPane; 
Snvatl JMenuBar jMenuBarl = new DMenuBarQ; 
orivate JMenu iMenuFile = new JMenuC;; ^^^^r\. 
lAllte SUenultem jMenuFileExit = new JMenultemQ , 
nrivate DMenu iMenuHelp = new DMenuu; 
pAlltl JMenuItem jMenuHelpAbout = new :Menuitem(); 
orivate JToolBar jToolBar = new DToolBarQ, 
private DButton openButton = new ^ButtonQ; 
private DButton saveButton = new DButtonC), 
private DButton jButtonS = new DButtonQ; 
private imagelcon imagel; 
private imagelcon imageZ; 

""Dn-OatfjMenuitei'openMenultem = new DMenuitemQ; 
lAvlte DMenuitem saveMenuitem = new jMenuitemQ; 
orivate DButton spectrumButton = new DButtonQ, 
lAllle DComboBox^unctionchoice = new DComboBoxQ; 
orivate D Label j Label 1 = new DLabelQ; 
pHvttl DDialog helpDialog = new DDialogO; 
orivate D Panel j Panel 1 = new DPanelQ; 
Sniate DButton helpDismiss = new DButtonQ; 
lAllll DScrollPane iScrollPaneZ = new.DScrollPaneO ; 
lAAle DEditorPane fielpPane = new DEditorPaneQ ; 
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nrivate JPanel aboutMessage = new jPanelQ; 
private 3Label jLabel2 = new ^LabelQ 
nr-ivate JLabel iLabelS = new DLabeKJ, . 
SrivSl GridLayout gridLayoutl = new GndLayoutQ; 
private 3Label jLabel4 = new DLabelQ; 
private component componentZ ; 
private DLabel ijLabelS = new ^LabelQ, 
private DLabel i Label 6 = new JLabelQ, 
private Border borderl; 
private Component components; . ^ 
private JLabel jLabel? = new JLabelQ; 
private D Panel 3 Panel 2 = new jPanelQ; 
nr-ivate DLabel statusBar = new DLabel U; „ 
lAATe DP?ogressBar progressBar = new DProgressBarQ ; 
oAvltl BordirLayout borderLayoutS = new BorderLayoutQ 
IrulTe ISfdlrLayout borderLayoutZ = new Border-Layout (); 
lAvltl DMenultem iMenuitemZ = new DMenultemO ; 
private DLabel jLabelS = new DLabel Q; 
private component component!; 
private Component componentA; 
nrivate D Panel errorMessage = new DPane U . 
grivltl GridLayout gridLayout2 = new GridLayoutQ ; 
private DLabel errorLabell = new DLabel 0, 
private DLabel errorLabel2 = new DLabel Q, 
private AxisDialog axisDialog; 

//construct the frame 
public converterFrameC) 

^ enabl eEvents (AWTEvent . window_event_mask) ; 
try 

Qnprt rumiabl e . setModel C model Parameters ); . i, • r>> 

ffi?a?Imetlrl!addTableModelListener( new TableModel Listener Q 

public void tableChanged(TableModel Event e) 

^ if (spectrumTable.getselectedColumnO >= 0 && e.getColumnQ == 

spect r umTabl e . getsel ectedcol umn () ) 
nomogram_actionO ; 

spectriiiiTable.getTableHeaderQ.setReorderingAllo^^^ false ); 
spect rumxabl e . setcol umnsel ecti onAl 1 owed ( t rue ) , 

Li stsel ecti onModel ^^ingle_selection^) ; ^^^^^ ^^^^^^^^ ^ . ^^^^^^^ 
new TableColumnModel Listener 0 

^ public void col umnsel ecti onchangedC Li stsel ecti onEvent Ise ) 
^ nomogram_action() ; 

public void columnMarginChangedC ChangeEvent ce ) 

Lblic void columnMovedC Tabl eCol umnModel Event tcme ) 
{ 

Jublic void columnRemovedC TableColumnModel Event tcme ) 
^ Page 2 
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Jublic void columnAddedC Tab! eCol umnModel Event tcme ) 



axUoialog = new AxisDialogC this. "Logicle scale", false ); 
reip?iSl°lefpaaP('converterFr^^^ "help.ht.V ) ); 

helpDialog.packO; 



catch (Exception e) 



^ e.printStackTraceO; 
} 

feeTo?d^5brnftir?h?^^^ Exception i 
W;c?n(edu"i?anford.facs.1oglikeXonverterFrame.class.get^^^^^ 

Lg;;c&u?I?anf ord . f acs . 1 ogl i ke XonverterFr^^^ 
illgS(edu"ILnf ord . facs J ogl i ke . converter Fram^ 

//setlconi™age(Toolkit.getDefaultToolkit0.createi.ageCConverterFra.e.class.getResou 

^"^;jreStSre'PB^anel) this,getContent ; 
comoonentZ = Box.createVerticalstrutes; , 
boX?l = BorderFactory.createEmptyBorder(10.15.10.25) . 
components = Box.createVerticalStrutCS) 
componentl = Box.createVert-icalStrutCS) . 
component4 = Box.createVerticalStrutCS). 
contentPane.setLayout(borderLayoutl); 
this setsizefnew Dimension(400. 300;;, ^.^^..n. 
thililtTitleC'Logicle Fluorescence compensation ). 
iMenuFi le . setTextC'Fi 1 e") ; 

^SPilfE^i'tilfdlSionLi^^ei^rCnew ActionListenerO { 
^ pSblic void actionPerformedCActionEvent e) { 
jMenuFileExit_actionPerformedCe; ; 

} 

jMenuHelp.setTextC'Help';); , 

public void actionPerformedCActionEvent e; i 
about_actionCe) ; 

> 

^ public void actionPerformedCActionEvent e) 

^ data_actionCe); 
} 

S^ptrSSSIS^So^lfS^Slf^ien PCS File"); 
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^ public void action Performed (Action Event e) 

^ save_action(e) ; 

} 

saveButton.setEnabled (false); c-ii*"^- 
saveButton.setToolTipTextC'Save Extended FCS File ), 

jButton3.setlcon(image3); . * ^w««-h A,-+-:oni i«:i-pnprn 

3 Buttons . addActi onLi stener (new ] ava. awt . event .Acti onu stener u 

public void actionPerformed(ActionEvent e) 

^ help_action(e) ; 
} 

j Buttons . setToolTi pText ("Hel p") ; 

spectrumMenuitem.seticon(new -:r>o"^^^■ 
imageicon (Conve rter Frame .class. getResource( spectrum . : pg ) ) ) , 

spectrumMenultem.setTextC'Spectrum . . . ); A^+iAnH^tPnerf') I 

spectrumMenuitem.addActionListener(new java.awt. event. Actionusteneru 1 
public void acti on Performed (Acti onEvent e) i 
spect rum_acti on (e) ; 

})1 

openMenultem. set Enabled (fal se) ; 
openMenuitem.setlcon(imagel) ; _ 

openMenultem. setTextC'Data .. . ); ow«.„i- Ari-inni ■i«;tf»nern 

openMenultem . addActi onu stener(new java . awt . event . Acti onLi stener u 

^ public void actionPerformed(ActionEvent e) 

data_action(e); 

} 

}) ; 

saveMenuitem.setEnabled (false) ; 
saveMenultem.setlcon(image2) ; 

saveMenuitem.setTextC'save AS ... ); A^t-i^ni -iQi-^nprn 

saveMenultem . addActi onLi stener(new java . awt . event .Acti onLi stener C) 

^ public void acti onPerformed (Acti onEvent e) 

^ save_action(e) ; 

} 

functionChoice.setFont(new java.awt.Font("Dialog" , 1- If)): 

f uncti onchoi ce . setMaxi mumSi ze(new Di mensi on (32767 ; , 

f uncti onchoi ce . setMi ni mumSi ze(new Di mensi on (138 . / 4 J ! 

f uncti onchoi ce.setPreferredSize (new Dimension(142, Z7;j,, 

functionChoice.setToolTipTextC'Select a Data Transformation ); 

f uncti onchoi ce.setModeKf uncti onModel); Ari-inm -i^rpnprfi 

f unctionchoi ce . addActi onLi stener(new java . awt . event . Acti onLi stener (.) 

^ public void actionPerformed(ActionEvent e) 

^ function_action(e) ; 
} 

spectrumButton.setMaximumSize(new Dimension(51, 27)); 
SDectrumButton.setMinimumSize(new DimensionQbi, 2/jj; 
spectrumButton.setToolTipTextC'Load spectrum Matrix ); 
spectrumButton.seticon(new 
imageicon (ConverterFrame.cl ass. getResource( spectrum. 3pg ))), 
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spect rumButton . addActi onLi stener (new j ava . awt . event . Acti onLi stene r u 

^ public void acti onPerformed (Acti onEvent e) 

spect rutiuacti on (e) ; 

} 

jLabell.setPreferredSize(new Dimension (17, 17)) ; 
spect rumrabl e . setRowSel ecti onAl 1 owed (fal se) ; 

hel pDi smi ss '. addlcti onLi stener (new j ava . awt . event . Acti onLi s tene r () 
^ public void acti onPerformed (Acti onEvent e) 

^ help_dismiss(e); 

} 

SeipDialog.setTitleC'Logicle Help"); ^ , 

helpDi alog.getContentPaneO .setLayout(borderLayout2) ; 

jMenultemZ.setTextC'Help . . ."); . * a^4.-;««i -ici-anorn 

|Menuitem2 . addActi onLi stener(new ] ava . awt . event . ActT onLi stener () 

public void actionPerformed(Acti onEvent e) 

help_action(e); 

} 

heipPane.setPreferredSize(new Dimension (400, 400)); 

l;l^^P^n^JSTei^("ih^^^Ar\^. <head>\r\n\r\n </head>\r\n <body>\r\n <p>\r\n 

Logic! e Data converter " + „ ^ , ^ /u^ i v \ "^ 
\r\n </p>\r\n </body>\r\n</html>\r\n ); 
he! pPane . setcontentType ("text/html ") ; . 
iScrollPaneZ.setPreferredSizeCnew Dimension (400, 400;;, 
lLabel2.setFont(new java.awt. Font ("Dialog" , 1, 14;;; 
1 Label 2 . setAl i qnmentY( (f 1 oat) 0.0); ^...-r.n^ 
1 Label 2 . setHori zontal Al i gnment (Swi ngconstants . CENTER) ; 
1 Label 2 . setHori zontalTextPosi ti on(Swi ngconstants . LEADING) ; 
lLabel2.setlconTextGap(0) ; , 
iLabel2.setText("Logicle Fluorescence compensation ); 
1 Label 3. setAlignmentx( (float) 0.5); 
3 Label 3 . setHori zontal Al i gnment (Swi ngconstants . CENTER) ; 
j Label 3 . setlconTextGap(O) ; 
] Label 3. setText ("Wayne A. Moore ); 
aboutMessage . setLayout (gri dLayoutl) ; 
gridLayoutl.setColumn5(l) ; 
gridLayoutl.setRows(ll) ; , 

1 Label 4 . setHori zontal Al i gnment (Swi ngconstants . CENTER) , 

iLabel4.setlconTextGap(0) ; 

] Label 4. setText ("October 2002"); . , „ ^ 

j Label 5. set Font (new java.awt . Font ( Dialog , 2, 10)), 

1 Label 5 . setHori zonta I Al i gnment (Swi ngconstants . CENTER) ; 

] Label 5 . setHori zontalTextPosi tion(Swi ngconstants . LEADING) ; 

jLabel5.setlconTextGap(0); . 

1 Label 5. setText ("copy right 2002, by The Board of Trustees ). 

iLabel6.setFont(new java.awt . Font ( Dialog , 2, 10)); 

1 Label 6 . setHori zontal Al i gnment (Swi ngconstants . CENTER) ; 

1 Label 6 . setHori zontal Text Posi ti on (Swi ngconstants . LEADING) ; 

lLabel6.setlconTextGap(0); 

] Label 6. setText ("of the Leland Stanford Jr. university ); 
aboutMessage. setBorder(b9rderl); , 
aboutMessage . setMi ni mumSi ze (new Di mensi on (60 , 40) ) , 
aboutMessage. setPreferredSize(new Dimension(300, 250;;; 
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i Label 7 . setHo ri zontal Al i gnment (Swi ngConstant s . CENTER) ; 
]Label7.setText("David R. Parks"); 
StatusBar. setTextC" "); 
i Panel 2 . setLayout(borderLayout3) ; _ 

oroaressBar.setString ("compensate ) ; ^^^,Tcrl^. 
?[aBel 8 setHori zontaUl i gnment (Swi ngConstants . CENTER) ; 
]Label8.setText("version 0.6"); 
errorMessage.setLayout(gridLayout2), 

gri dLayout2 . setcol umns (1) ; 
qridLayout2 .setRows(O) ; 
errorLabell.- — — " 
errorLabel2 . 

errorLabel2.- . ^- 

qMenul.setText("view ); _ 

ffiSSJSi: IdSIS^aU l«n^knew jav. , .wt . .v.nt . ActionLi stenerO 
public void actionPerformecl(ActionEvent e) 

^ seal cacti on (e); 
} 

iTool Bar . add (spect rumButton , null ) ; 
jToolBar.add(openButton) ; 
iToolBar.add(saveButton) ; 





iToo iBar.aaat^TuncxionChoice, nul I), 
iNienuFi 1 e . add(spectrumMenuitem) ; 
] MenuFi 1 e . add (openMenuitem) ; 
j Menu File. add (saveMenult em) ; 
iMenuFile.addSeparatorO; _ 
j MenuFi le. add (j MenuFi leExit); 
jMenuHelp.add(jMenuitem2) ; 
jMenuHelp.addseparatorO ; 
iMenuHelp.add(iMenuHelpAbout) ; 
]MenuBarl.add(qMenuFile) ; 
jMenuBarl.add(:jMenul) ; 
iMenuBarl . add (] MenuHel p) ; 

this.set3MenuBar(jMenuBarl); Kin^xu^ ■ 

contentPane. add (jTool Bar, BorderLayout .north) 
contentPane. add (j scroll panel. Border Layout. CENTER) , 
contentPane. add (] Panel 2, BorderLayout .SOUTH) . 
i panel 2 . add (statusBar , BorderLayout .CENTER) , 
1 Panel 2. add (prog ressBar, BorderLayout. EAST) ; 
n^rrollpanel aetviewport().add(spectrumTable, null); 
illpSUloS getconlXaJeO .add(5panell, BorderLayout. south) ; 

jScrollPaneZ.getviewportO-addChelpPane, nul Ij; 
aboutMessage . add ( j Label 2 , null); 
aboutMessage . add (component4 , nul U 5 
aboutMessage . add( j Label 8 , nul 1 ; ; 
aboutMessage . add(componentl , nul I J ; 
aboutMessage. addCjLabel 3 , null) ; 
aboutMessage. addqLabel7, null); 
aboutMessage . add (] Label 4 , null); 
aboutMessage . add (component2 .null); 
aboutMessage. addCqLabelS, null); 
aboutMessage. add (] Label 6, null) ; 
aboutMessage . add (components , null); 
errorMessage . add(errorLabel 1 , nul I ) ; 
errorMessage.add(errorLabel2, null); 
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jMenul.add(jMenulteml) ; 
//File I Exit action performed 

public void jMenuFileExit.actionPerformed(ActionEvent e) { 
system. exit (0) ; 

//Help I About action performed . ^ . . 
//overridden so we can exit when window is closed 
protected void processWindowEvent(windowEvent e) i 

super .processWindowEvent(e) ; 

if (e. get ID O == windowEvent.WiNDOW_CLGSiNG) { 
jMenuFi 1 eExi t^acti onPerf ormed(nul 1 ) ; 

} ^ 

private static final String[] stringArray = new String[0]; 
private static final int n^parameters = 2; 
private static final int fuzzy_bits = 8; 
private FCSFile fcs; 

private Strinq[] sensorNames = new StnngLOj; 
private int[][] original Data; 

[] newData; 
e[] [] rawData; 
e □ [] compensatedData ; 
e[] [] compensationMatrix; 
private int progress; . 
protected LogicleFunction logicle; . ^ ^- /^ 

protected LogicleFunction logarithm = new LogicleFunctionQ ; 

public class LogicleFunction 

^ double decades, minimum, maximum; 

public doubled transform ( 
int column, 

FCSParameter parameter ) 

^ doubled lookup = new double[ oarameter .getRangeQ ]; 
maximum = parameter .getMaximumQ ; 
minimum = parameter. getMi nimumQ ; 
decades = parameter .getDecadesO ; 
lookup [0] = minimum; . 
interpolateC lookup, 0, lookup. length, maximum ;; 

return lookup; 

} 

} 

public final LogicleFunction logLinearsplice = new LogicleFunction Q 

public doubled transform ( 
int column, 

FCSParameter parameter ) 
^ double linear_range = model Parameters. getParameterC sensorNames. length, 

column ); . 

if (Double. isNaNC linear_range )) 

return super. transformC column, parameter ;; 

maximum = parameter .getMaximumQ ; n 
double ln_range = Math.logC maximum / linear_range ) + 2, 
decades = ln_range / ln_10; 
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minimum = maximum / Math.expC ln_range ); 

doubleC] lookup = new doub1e[ parameter. getRangeQ ]; 
];;WiteriSikup: O. lookup.length. maximum ); 

MiEh'logC irn'el'JaSge^rSinimum ) * (double) parameter. getRangeQ / 
parameter .getDecadesCW^^nUO ^ ' p,,,„eter.getRange() - 1) 

^ double delta = (lookup [splice+1] - lookup [splice-l])/2; 
double base = lookup [splice] ; 
for (int i = 0; j < splice; ++:) 
lookup[jJ = base - (splice - ]) * delta, 

} 

return lookup; 

} 

}; 

public final LogicleFunction bi exponent! al = new LogicleFunction Q 

public double [] transform ( 
int column, 

FCSParameter parameter ) 
^ double width = modelParameters.getParameter( sensorNames. length, column ); 
if (Double. isNaN( width )) 

return super .transform( column, parameter 

decades = parameter.getDecadesO ; 
width /= decades; 

double extra = model Parameters. getParameter( sensorNames. length + 1, column 
if (Double. isNaN( extra ) II extra < 0) 

inf^Sro 2'(int) Math.round( extra * parameter .getRangeQ / (extra + decades) 
if (zero > 0) 

^ extra = zero * (decades + extra) / parameter .getRangeQ ; // so extra falls 

exactly on channel bound ry 

decades = extra * parameter. getRangeQ / zero, 

} 

maximum = parameter .getMaximumQ ; 
double positive_range = ln_10 * decades; 
minimum = maximum / Math.exp( positive_range ); 

doubleC] positive = new doubleC parameter. getRangeQ ]; 

?Sw?poli?i(%isitive, 0. positive. length, Math.exp( positive_range ) ); 

double negative_range = root( positive.range, width ); 
system. out. println( positive_range ); 
system. out. printlnC negative_range ); . 
doubleC] negative = new double [ positive. length ]; 

intfrpola?e(%e^ative, 0, negative. length. Math.exp( -negative_range ) ); 
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double s = Math.expC (positive_range + negative_range) * (width - extra / 

decades) ^)^.^^ i = 0; j < negative. length; 
negativetj] *= s; 

s = positive [zero] - negative[zero] ; 
for fint i = zero; j < positive. length; ++]; 
posi?ivelj] = tninimim * (positive[jl - negativeC]] - s); 

for (int 1=0; j < zero; 
positivetj] = -positive[2 * zero - ]]. 

return positive; 

} 

}; 

public final LogicleFunction hyperbolic = new LogicleFunction Q 

public double [] transform ( 
int column, 

FCSParameter parameter ) 

^ double linear.range = model Parameters. getParameterC sensorNames. length. 

column ) ; , , . 

if (Double. isNaNC linear_range )) 

return super .transform( column, parameter ), 

maximum = parameter. getMaxi mum() ; i/Math eV 

double scale_factor = linear.range / Cw^th.E - 1/Math.E; , 
double ln_range = Math.log( maximum / seal e_f actor J + ^, 
decades = ln_range / In 10; 
minimum = maximum / Math.exp( ln_range ), 

double[] lookup = new double[ parameter. getRangeO ]; 

iS;elipilate('looK';"o: lookup. length, maximum / scale_factor ); 

for fint i = 0: j < lookup. length; ++j) ^ ^ ^ 
lookSpEjj = (iookupCj] - l/looRup[j]) * scale_factor; 

return lookup; 

} 

}; 

public final LogicleFunction logPlus = new LogicleFunction Q 

^ public doubled transform ( 
int column, 

FCSParameter parameter ) 
^ double parameterl = model Parameters. getParameterC sensorNames. length, column 
^' double parameterZ = model Parameters. getParameterC sensorNames. length + 1. 

column ); . , ^ i ^^ 

if CDouble.isNaNC parameterl )) 

return super. transformC column, parameter J, 
if CDouble.isNaNC parameter2 )) 
parameter2 = 0; 

double[] lookup = super. transformC column, parameter ); 
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. 1 *^£th^loaHarameter.getMinimuir.O ) / ln_10; 

?nr fiJt i =0; j < lookup. length; ++3) 
lookSpLj] parameterl * j + parameter2; 

return lookup; 

} 

}; 

private class Model Parameters 
extends AbstractTableModel 

^ private double[][] parameters = new double[0][]; 
public int getRowCountO 
^ return parameters. length; 
} 

public int getColumnCountQ 
^ return sensorNames. length; 
} 

public String getcol umnName ( 
int columnindex ) 

^ return sensorNames[ columnindex ]; 
} 

public Class getcol umnCl ass ( 
int columnindex ) 

^ return super .getcol umnCl ass ( columnindex ); 

} 

public boolean i scel 1 Edi tabl e ( 
int rowindex, 
int columnindex ) 

{ 

return true; 

} 

public object getvalueAt ( 
int rowindex, 
int columnindex ) 

^ double value = parameters[ rowindex ][ columnindex ]; 
if (Double. isNaN( value )) 
return null; 

return string.valueof ( value ); 

} 

public void setvalueAt ( 
Object avalue, 
int rowindex, 
int columnindex) 

^ double dvalue; ^ . ^ 

if revalue instanceof Number; r\ 
dJilSe = (cSumber) aValue) .doubleValueQ ; 
else if (avalue instanceof Stnng)^^^ 
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{ 

try 



LI J 

^ dvalue = Double.valueOfC (String) avalue ) .doublevalueO ; 
catch (Number FormatExcepti on ex) 

{ . . _ 



dvalue = Double. NaN; 
} 

} 

else 



if''"(?i!;lll.??ilal(-d;Siue) M rewinds < sensorNan,es.le„gth) 
if ( rowindex == columnindex) 
J. .--I 1 . 



dvalue = 1; 
else 
dvalue = 0; 



naramPtersF Towindex ][ columnindex ] = dvalue; 
ffrl?lblelellSSa;id( rowindex. columnindex ); 

public void setParameters ( 
double[][] parameters ) 

this. parameters = parameters; 
this.fi reTableStructureChangedO, 



public double getParameter ( 
int rowindex, 
int columnindex ) 



^ return parameters [rowindex] [columnindex] ; 
} 

}; 

private Model Parameters model Parameters = new Model Parameters () ; 

^^s^^^n;?Pl^'"^o5's^=^^ - ^og x . b- >; 

functionNames ); i„in - wiai-h loaC 10 

rrrvftl SS%rr;;'g 'iS~iRROR r^Vel^taVr^ not compatible with the 

'PpHCaw'^inal int fcs_timer_tick = 200; 
private final Random random = new RandomO, 

private void interpolate ( 
double [] func, 



int i , 
int n, 
double X ) 



Snt m = n / 2; 

double y = Math.sqrt( func[i] * x ); 

func[ i + m ] = y; 
if (m > 1 ) 
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interpolateC func, i, m, y ); 
interpolateC func, i + m, m, x j, 

} 

} 

orivate void shuffle ( 
double [] sequence ) 

^ for (int i = 0; i < sequence. length; 

^ int j = i + random. nextintC sequence. length - i ); 
double t = sequence [i]; 
sequence [i] = sequence [3]; 
sequenceCj] = t; 

} 

} 

public void invert ( 
double[][] matrix ) 

^ int row = 0, col = 0, n = matrix. length; 

int pivot[] = new int[n]; 

int row_index[] = new int[n]; 

int col_index[] = new int[n]; 

for (int 1=0; i < n; ++i) 

^ double big = 0; . 

for (int ] = 0; 3 < n; ++^J 

^ if (pivot[j] != 1) , ,v 
for (int k = 0; k < n; ++k) 

^ if (pivot[k] « 0) 

^ double abs = Math.abs(matrix[j] [k]) ; 
if (abs >= big) 
{ 

big = abs; 
row = i; 
col = k; 

} 

hy^ril ^Slw°yiegalirgumentException( "matrix is singular" 

} 

> . r -.1 

++pivotLcoi J ; 
row_index[i] = row; 
col_index[i] = col ; 

if (row != col) . , - 

for (int k = 0; k < n; ++k) 

^ double t = matrix[row][k]; 
matrix[row][k] = matnx[col] [k] ; 
matrixCcol][k] = t; 

} 

'^Sr^jK^l'SilA^^uS^ntExcept^ "matrix is singular" ); 
double inverse = l/matnx[col] [colj^^ 



{ 
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matrix[coT][col] =1; 

for (int j = 0; D < n; ++3) 

matrix[col][j] *= inverse; 
for (int j = 0; ] < n; ++3) 

if (j != col) 

^ double t = matrix[j][col] ; 
matrix[j][col] =0: 
for (int k = 0; k < n; ++k) ^ 
matrix[j][k] -= matnx[col] [k] - t, 

} 

} 

for (int i = n-1; i >= 0; —i) 
if (row_index[i] != col_index[i]) 
for (int j = 0; j < n; ++]) 

^ double t = matrixm[row_indexm]; . . rill- 
matrixCj] [row_indexIi]] = mat nx []] [col .index [1]], 
matrix[j] [col_index[i]] = t; 

} 

} 

double root ( 
double b, 
double w ) 

double Xlo = 0; 

double Xhi = b; 

double Flo = Double. negative_infinity; 

double Phi = w * 2 * b; 

double d = (xlo + xhi)/2; 
double Dx = Math.abs( Xlo - Xhi ); 
double Dx_last = Dx; 

double Fb = -2 * Math.log(b) + w*b; 
double F = 2 * Math.log(d) + w*d + Fb, 
double DF = 2/d + w; 

if (w == 0) 
return b; 

for (int 1=0; i < 100; ++i) 

^ if (((d-Xhi)*DF-F)*((d-Xlo)*DF-F) >= 0 

II Math.abs( 2 * F ) > Math.abs( Dx_last * df )) 

^ system. out. println( "bisection " + Dx ); 
dx = (Xhi - xlo)/2; 
d = xlo + Dx; 
if (d == xlo) 
return d; 

} 

else 

system. out. println( "Newton " + Dx ); 
double t = d; 
d -= Dx; 
if (d == t) 
return d; ^_ 
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if (Math.absC Dx ) < IE- 12) 

return d; 
Dx_last = Dx; 

F = 2 * Math.log(d) + w*d + Fb; 
DF = 2/d + w; 
if (F < 0) 

xTo = d; 
Flo = f; 

} 

else 

^ Xhi = d; 
Fhi = f; 

} 

throw new illegalStateExceptionC "exceeded maximurmterations ). 

} 

void make_busy () 

^ spectrumButton . setEnabledC f al se ) ; 
spectrumMenuitem.setEnabledC raise 

openButton . setEnabl ed( f al se ) ; 
openNienuitem. setEnabledC false ), 

saveButton. setEnabledC false ); 
saveButton. setEnabledC raise 

progressBar.setEnabledC false ); 

} 

void make_ready C) 

^ spectrumButton. setEnabledC true ); 
spectrumMenuitem.setEnabledC true ), 

if CsensorNames != null) 

^ openButton. setEnabledC true ); 
openMenultem. setEnabledC true 

if Cfcs != null) 

^ saveButton. setEnabledC true ); 

saveButton. setEnabledC true ); ^ v. 
statusBar.setTextC fcs.getFileC) .getNameC) ). 

} 

e1 se 

statusBar. setTextC "" ); 



progressBar.setvalueC 0 ); 
progressBar.setEnabledC false ), 



} 
v< 
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^%a?rrSle^|4UlirE» .stopCellEditingO : 

make^usyO; 

''LfferedReader br = new Buff eredReaderC new FileReaderC 
fileChooser.getselectedFileO ) ): 

br.readLineO ; _ „ 

'sJnng^^^^^^^^^^ st = newStringTokenizerCbr.readLi \t ): 

ArrayList al = new ArravustQ , 
while (st.hasMoreTokensO) 

s.Slo?raler;TtHn|c?) ai.toArrayC strin^Array ); 
doubled [] parameters - new doubletsensorNames. length * 
2] [sensorNames. length] ; 

fSJ clr'" 0; r < sensorNames. length; ++r) 
^ - new stringTokenizerC br.readLineO, "\t" ); 
-S;Eeyrftc3V STrMt s^SUo^enO ); 

} 

far (• r < parameters. length; . .-.^ 

for'fint c = 0; c < sensorNames. length. ++0 
^°paraSeters[r][c] = rouble. NaN; 

model Parameters. setParametersC parameters j, 

'^iV'V- o"^^? < sensorNames. length; -f+r) 

if (fcs getParameterC sensorNames [r] ) == null) 

^ fcs = null ; 
break; 

} 

catch (Exception ex) 
^ ex.printStackxraceO; 
make_readyO ; 

} 

} 

FCSFile fcs data; 
FCSHandler fcs_handler; 
lOException fcs.exception; 

Timer fcs_timer = new Timer( fcs_timer_tick. new ActionListener () 

^ public void actionPerformed ( 
ActionEvent action ) 

^ progressBar.setvalueC fcs_handler.getEventO ); 
} ); page 15 
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void data_action(ActionEvent e) 

{ 

try 

^ File f = fileChooser.getselectedFileO; 

fcs data = new FCSFiTeC f ); i«^„*.u. 
for r?nt i = 0: i < sensorNames.length, 
if^(fcs_data getParameterC sensorNames[i] ) == null) 

^ joptionPane.showMessageDialogC thi|> SPECTRUM_ERROR. 
■""P^.. Error" , joptionPane.ERROR_MESSAGE ), 
return; 

} 

make_busyO; n ^^ 

progressBar-setMinimumC 0 ), 

gfSrrelfBfr'.lftrai^^y Us.data.getTotal () ): 

newThreadC data_thread. "FCS reader" ).startO: 
fcs_timer. start O; 

catch (Exception ex) 

^ ex.printstackTraceO; 
system. exit ( 1 ) ; 

} 

} 

} 

Runnable data_thread = new RunnableQ 

^ public void run Q 

^ Thread.currentThreadO.setPriorityC Thread.MlN_PRlORlTY ); 

fcs_exception = null; 
int n = fcs_data.getTotaiq, 

n, = fcs_data.getParaffletersp, 
illt[][] int_value = new int[m][n]; 
try 

^ fcs-handler = J"-data.|et^nP™ 
^°fo Cint i 2'o; fcSandler.hasMoreValuesO; ++3) 

^°int !!alueCj][i] = fcs_handler.readvalueO ; 
fcs_handler.closeO; 
original Data = int_value; 

catch (lOException ex) 

^ fcs_exception = ex; 
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swi nguti 1 i ti es . i nvokeLater ( data_compl ete ) ; 

} 

}; 

Runnable data_cotnpl ete = new Runnable O 
^ public void run O 



^ fcs.timer.stopO; 

if (fcs_exception != nuiu 

^ errorLabell.setTextC SirFifeC^re^PaS 

errorLabelZ.setTextC fcs^data.getFi leu -y ^ errorMessage, 
!Sp?ionPane.showMessageDialog( ConverterFrame.tms. ^ ^ . 



} 

else 

^ fcs = fcs_data; 
make_readyO; 

} 

}; 



void save_action(ActionEvent e) 
{ 

try 

if (f.existsQ) ^ . „,. 

' objectC] message = { f. "already exi«.. do you want to replace n. 

if (joptionPane.OlC-OPTION . , r^- messaqe, "overwrite", 

^gs?ro:ssr.Kr^s?;!is!SoSi?inra^e.3*4NiN<^.ss«E 



return; 

} 



int n = fcs.getTotalO; 
int m = sensorNames. length, 
rawData = null ; 
compensatedData = nu n , 
newData = nu 11 ; 

compen5ationMatrix.= new double[m] [m] ; 
for (int i = 0; i < m; 

'"JoiiSSsaj^onSUr^xHlbl" »odelParan,eters.getPara.eter( 1, j )i 

try 

^ invertC compensationMatrix ); 
catch (illegalArgumentException ex) 
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30ptionPane.showMessageSaloif^hii!°"sJectr^ Matrix is Singular", 
-spectrum Erro^.^.^^^^^^^^^^^^^^^^^^ ^. 
return; 

} 

make_busy() ; 

progress =0; . . , ^ \ 
progressBar.setMTmtnumC 0 * „ n. 

progressBar.setMaximutnC 2 * m + m m j, 
progressBar.setValueC 0 ); 
progressBar.setEnabledC true ); 
statusBar.setTextC "working ); 

fa data = new FCSFileC fileChooser.getSelectedFileQ ); 
fctdSt.geSTextslgmentO.readFromC fcs.getTextsegmentQ ). 

new Thread ( xf rm_thread, "Transform" ).startO; 

catch (Exception ex) 

^ ex.printstackTraceO ; 
} 

} 

} 

Runnable xfrm_thread = new RunnableQ 
^ public void run Q 
try 

^ int n = fcs.getTotalQ; 
int m = sensorNames. length; 

rawData = new double[m] [n] ; 
for (int 1=0; i < m; 

^ FCSParameter p = fcs.getParameter( sensorNames [i] ); 
int[] input_data = originalData[p.getindexO-l] ; 
double [] output_data = rawData[i]; 
if (p.isLogO) 

^ double [] antilog = new double [ p.getRangeO ]; 
r„UrS§taJe( Sn?nS^?'S™p'aetRangeO. p.getMaximumO ); 
double □ fuzz = new double [ 1«FUZZY_BITS ]; 
iSteSilate^ fuzz. 0. 1«fuzzy_bits . antilog[l]/antilog[0] ); 
shuffle( fuzz ); 

'°ou^p!!LditI[ji JaStiiog[1nput_data[j] ] * fuzz[ j & 

((1«FUZZY_BITS)-1) ]; 

else 

^ double g = p.getGainO; 

for (int j = 0; j.< n; ++D) ^ , ^. 
output_data[j] = input_data[3] / g, 
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sS^iliiies.invokeLaterC progress_made ); 

} 

compensatedData = new doubleM [n] ; 
for (int i = 0; i < m; ++^) 
^ doubleL] output_data = cotnpensatedData[i] ; 
for (int j = 0; j < m; ++3; 

'^SiiCLd^tacK ?nSit_ra?aCk] * coefficient; 
sl?SgStiliiies.invokeLater( progressjnade ); 

} 

} 

newData = new int[m][n]; 
for (int i = 0; i < m; 

^ FCSParameter oldp = fcs.getParameter( sensorNames[i] ); 
^FfsJlrlSfr newS =/cs_data.addPa™^^^ 

"nrp:irtftrr?Sre'( oid^!^gUtribute( "Sp". "s" ) ); 

newp . setBi t s ( ol dp . getBi t s O ) ; 

SlIftSLro^drg"^^^^^^^^ oldp.getMini.u.0 ): 
newp.setGain( oldp.getGainO ); 

doubl e [] i nput_data = compensatedData [i ] ; 
int[] transformed_data = newData[i], 
if (newp.isLogO) 

^ w^nMon lookuD = loqic1e.transfortn( i, newp.); 
Sewp.lPtDecades? logi?le. decades, logicle.mmmum ); 
for (int 3=0; j < n; ++]) 

^ k = Arrays.binarysearch( lookup. input_data[j] ); 

if (k < 0) ^ ^ ... ... 

transformed_dataLDJ = -^--'-i 

^ transformed_data[j] = k; 

} 

} 

else 

^ double g = newp.getcainO ; 

'^JrSforieLdUj] (i'nS' Math. round ( input_data[j] * g ); 

} 



} 



s;?Sgu?iliiies.invokeLater( progress_made ); 



SSrutili?its^nvSkrLker( xfrm_complete ); 
catch (Exception ex) 
^ ex . pri ntstackTraceO ; 
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} 

}; 

Runnable progressjnade = new Runnable O 

^ public void run Q 

^ progressBar.setvalueC progress ); 

Runnable xf rm_compl ete = new Runnable Q 

^ public void run Q 

^ progressBar.setMinimumC 0 ); 

S^°ol?fsllt?:rerMli??y ks_data.getTotalO ); 
mfr^U^feM^^^^ > = 

new ThreadC save_thread, "FCS writer" ).startO; 
fcs_timer.startO; 

} 

}; 

Runnable save_thread = new RunnableO 

^ public void run Q 

^ Thread. CurrentThreadO.setPriorityC Thread.MlN_PRlORiTY ); 
fcs_exception = null; 
try 

^ fcs handler = fcs_data.getoutputlteratorO ; 
for Cint i = 0; fcs.handler.hasMoreEventsQ , 

^ for Cint i = 0; j < originalData.length; ++j) 
^°fcs handler. wHteValueC oj-^ 9^ "^iJg? ^ U )• 
for cTnt i = 0; j < newData . 1 ength , ++3 J 
fcslhandler.w^iteValueC newDatatDl [i] ): 

fcs_handler.closeO ; 

catch (lOException ex) 

^ ex.printStackTraceO; 
fcs_exception = ex; 

swi nguti 1 i ti es . i nvokeLater ( save_compl ete ) ; 

} 

}; 

Runnable save.complete = new Runnable Q 
^ public void run O 

^ fcs_timer.stopp ; 
if (fcs_exception != null; 
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> 

make_readyO ; 

} 

private JMenu jMenul = new JMenuQ; 

private JMenultem jMenuiteml = new JMenuitemu; 

void function_action(ActionEvent e) 

^ if (spectrumTable.isEditingO) r^. 
spect rutnTabl e . getcel 1 Edi tor () . stopce 1 1 Edi ti ng U , 

for fint i = 0: i < sensorNames. length; 

for Cint i = 0; j < thi s . n_parameters ; ++]) • ^ ^. 

model Parameters. setValueAtC null, sensorNames. length + j, ^ ). 

} 

void nomogram_actionO 

^ int columnindex = spectrumTable.getselectedColumnQ ; 
if (columnindex < 0 1 1 f cs == null) 

return; 
FCSParameter p; 
try 

^ p = fcs.getParameterC sensorNames [columnindex] ); 

catch (lOException ex) 
{ 

return; 

} 

if (p.isLogO) 

^ LogicleFunction function = logicleFunctionsC 
functionChoice.getselectedindexC) ]; ^^i,.„„TnHov n ^• 

doublen right = function.transform( columnindex, p ), 
doublet] left = new double[ p.getRangeQ ]; 
leftrO] = function. mini mum; . . 

interoolateC left, 0, left. length, function. maximum ), ron . 

axisE?alog nomoReiative = new aouble[axisDialog.nomoLabel .length] [2] . 
for (int i = 0; i < axisDialog.nomoLabel .length; ++i) 

^ double X = Double. parseDouble( axisDialog.nomoLabel [i] ); 
int k = Arrays. binarysearch( left, x ); 

'^aiisSialog.nomoRelative[i][0] = (double) (-k-1) / (double) left.length; 
^^lxisDialog.nomoRelative[i][0] = (double) k / (double) left.length; 
k = Arrays. binarySearch( right, x ); 

'^axisSialog.nomoRelative[i][l] = (double) (-k-1) / (double) right. length; 
^^lxisDialog.nomoRelative[i][l] = (double) k / (double) right. length; 
axi sDi al og . repai nt () ; 

} 

else 

^ double g = p.getGainO; 
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} 

void help_action(ActionEvent e) 

^ helpDialog.setLocationRelativeToC this ); 
helpDialog.setvisibleC true ); 

} 

void help_dismiss(ActionEvent e) 

^ helpDialog.setvisibleC false ); 
} 

void about_action(ActionEvent e) 

^ jOptionPane.showMessageDialogC this, aboutwessage. "About Logicle", 
30pti onPane . PLAIN_MESSAGE ) ; 
public void tableChanged(TableModel Event e) 

^ /**@todo implement this j avax. swing. event. TableModel Li stener^^^ 
throw new java.lang.unsupportedoperationExceptionC Method tablechangeau not yet 

im^l emented."); 

void scale_action(ActionEvent e) 

^ axisDialog.setVisibleC true ); 
axisDialog.toFrontO ; 

} 
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borderLayoutl^.BorderLayout, false, false, , ,false.<default> 
imagel.imageicon, false, false, , ,false,<default> 
image2,imageicon, false, false, , .false,<default> 
imageS.imageicon, false, false, , .false,<default> 
iButtonl,D Button, false, false, . .false,<default> 
]Button2, 3 Button, false, false, , ,false,<default> 
] Buttons, JButton, false, false, , ,false,<default> 
iMenuBarl,JMenuBar, false, false, . ,false,<default> 
iMenuFi 1 e , JMenu , f al se , f al se , , , f al se , <dgaul t> 
TMenuFileExit.JMenuitem, false, false, , , f al se , <detau I t> 
iMenuHelp,DMenu, false, false , .false,<default> 
1 MenuHel pAbout , DMenuitem , f al se , f al se , , , f al se , <detau 1 1> 
iMenuiteml,DMenultem, false, false, , ,false,<detau it> 
:Menultem2,JMenultem,false,false, , ,false,<default> 
i scrol 1 Panel, JScroll Pane, false, false, , ,false,<default> 
iTool Bar, JTool Bar, false, false, , ,false,<default> 
menuspect rum, JNienultem, false, false, , .false,<detau it> 
spectrum_table,JTable, false, false, , .false,<default> 
SDectrumChooser,DFileChooser, false, false, . , f al se , <detau I t> 
statusBar, J Label, false, false, , ,false,<default> 
[iconNames] 
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package edu . Stanford . f acs . 1 ogl i ke ; 

import qava.awt.*; 
import java.awt. event.*; 
import java.io.BufferedReader; 
import :java.io. File; 
import ;java.io. FineReader; 
import :|ava.io.lOException; 
import lava.util -ArrayList; 
import ^ava.util .Iterator; 
import lava.util .Map; 
import qava.util .StringTokemzer; 
import qavax. swing.*; 
import ]avax. swing. table.*; 

import org. i sac.*; 



Parks_et_al_12.txt 



/ 



** 

* <p>Title: Log Like Data Transform</p> 

* <p>Description: </p> 

* <p>Copyright : copyright (c) 2002</p> 

* <p>Company: Stanford Uni versity</p> 



* ©author Wayne A, Moore 

* ©version 1.0 

V 

public class ConverterFrame extends JFrame { 
private J Panel content Pane; 
private JMenuBar jMenuBarl = new 3MenuBar(); 
private JMenu jMenuFile = new JMenuO; 
private JMenultem jMenuFileExit = new JMenultemO; 
private JMenu jMenuHelp = new JMenuQ; 
private jMenuitem jMenuHelpAbout = new JMenultemQ ; 
private JToolBar jToolBar = new JToolBarQ; 
private DButton jButtonl = new JButtonO; 
private JButton jButton2 = new JButton(); 
private DButton ]Button3 = new JButtonQ; 
private imageicon imagel; 
private imageicon image2; 
private imageicon imageS; 
private DLabel statusBar = new JLabelQ; 
private BorderLayout borderLayoutl = new BorderLayoutC) ; 
private JScrollPane jScrollPanel = new JScrollPaneQ ; 
private JTable spectrum_table = new JTableO; 

private JMenultem menuSpectrum = new JMenultemQ ; ^ ^ r 

private jFilechooser spectrumchooser = new JFileChooser( System. getPropertyC 
"user. home", "" ) ); 

//construct the frame 
public ConverterFrameO { 

enabl eEventsCAWTEvent . WIND0W_EVENT_MASK) ; 

try { 
jbinitO; 

catch(Exception e) { 
e.printstackTraceC) ; 

} 

} ...... 

//component imtianzation 

private void jblnitO throws Exception { 

imagelcon(edu?stanf ord . f acs . 1 ogl i ke . ConverterFrame . cl ass . getResource ("openFi 1 e . gi f ") 
^ ' Page 1 
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imagS (edu"f?anf ord . f acs . 1 ogl i ke . ConverterFr^^^ 

linagS (edu"s?anf ord . f acs J ogl i ke . ConverterFr^^^ 

//seticonltnageCToolkit.getDefaultToolkitQ.createimageCconverterFrame. class. getResou 

rce(" [Your icon]"))); .. ^„ 

contentPane = (3 Panel) this.getContentPaneQ , 
contentPane . setLayout (borderLayoutl) ; 
this.setSizeCnew Dimension (400, 300)); 
this.setTitle("Log Like Data Converter ); 
statusBar.setText(" "}',,^ 
iMenuFile.setText("File ); _ 

iMenuFileExit.setTextC Exit ); -.^.^^^^r^ t 

iMenuFileExit.addActionListener(new ActionListenerO i 
public void actionPerformed(ActionEvent e) i 
jMenuFi 1 eExi t_acti onPerf ormed(e) ; 

iMenuHel p . setText ("Hel p") ; 

]MenuHelpAbout.setText( About ); r 
]MenuHelpAbout.addActTonListener(new ActionListenerO i 
public void actionPerformed(ActionEvent e) 1 
jMenuHel pAbout_acti onPerf ormed(e) ; 

jButtonl.seticon(imagel) ; 
iButtonl.setToolTipText( open File ); 
iButton2.setlcon(image2); 
iButton2.setToolTipText ("Close File ); 
3 Buttons. setlcon (images); _ 
j Buttons. setToolTipText( Help ); 

menuSpectrum.setText("Spectrum ); A^i-inm i^tpnprr^ I 

menuspectrum . addActi onLi stener(new ] ava . awt. event . Acti onLi stener u 1. 
public void actionPerformed(ActionEvent e) i 
spect rum_acti on (e) ; 

} 

}) ; 

iMenuiteml! addict i onLi java. awt. event. ActionListenerO 
public void actionPerformed(ActionEvent e) 
data_action(e) ; 

} 

iMenultem2 . setText("Save") ; 
iToolBar.add(jButtonl) ; 
iToolBar.add(jButton2) ; 
iTool Bar . add(] Buttons) ; 
iMenuFile.add(menuSpectrum) ; 
]MenuFile.add(jMenulteml) ; 
iMenuFile.add(]Menuitem2) ; 
3MenuFile.add(iMenuFileExit) ; 
iMenuHelp.add(qMenuHelpAbout) ; 
]MenuBarl.add(]MenuFile) ; 
iMenuBarl.add(jMenuHelp) ; 

this.set3MenuBar(jMenuBarl); .,«„-ru^ . 

contentPane. add(jTool Bar, Border Layout. NORTH) , 
contentPane . add(statusBar , BorderLayout . SOUTH); 
contentPane.add(jScroll Panel, BorderLayout. CENTER); 
iscrollPanel.getviewportO .add(spectrum_table, null) , 
^ page 2 
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//File I Exit action performed . ^ r 

public void jMenuFileExit_actionPerformed(ActionEvent e; i 

System.exit(O); 

//Help I About action performed . ^ ^ r 

public void jMenuHelpAbout_actionPerformed(ActionEvent e) i 

//overridden so we can exit when window is closed 
protected void processWindowEvent(windowEvent e; { 

super. processWindowEventCe); 

if (e.getiDO == windowEvent.wiNDOW_CL0SiNG) { 
jMenuFi 1 eExi t_acti onPerf ormed(nul 1 ) ; 

} ^ 

private static final String[] stringArray = new String[0]; 
private JMenultem iMenulteml = new DMenultemQ; 
private JMenultem jMenultemZ = new DMenultemQ; 

void spectrum_action ( 
ActionEvent e ) 

^ int returnval = spectrumChooser.showOpenDialogC this ); 
ifCreturnval == jFileChooser.APPROVE_OPTiON) 
{ 

try 

^ BufferedReader br = new BufferedReaderC new FileReader( 
spectrumChooser.getSelectedFileQ ) ); 
br. readLineO ; 

br. readLineC) ; , - ^ u ^. ^. 

stringTokenizer st = new stnngTokemzerC br. readLineCJ , \t 
ArrayList al = new ArravListO; 
while (st.hasMoreTokensQ) 

al.addC st .nextTokenO ) ; _ ^ ^ . ^ . 

String [] columns = (Strinq[]) al.toArrayC stnngArray ). 
Double[][] rows = new Doubl e [columns. length] [J ; 
for (int r = 0; r < rows. length; ++r) 

^ Doubled row = new Doubl e [col umns. 1 ength] ; 
rows[r] = row; ^ - ^. 

st = new StnngTokemzerC br.readLineO , \t j, 
for (int c = 0; c < col umns. length; ++c) 

^ row[c] = new Double( Double. parseDouble( st. nextTokenO ) ); 
} 

spectrum_table.setModelC new DefaultTableModel ( rows, columns ) ); 
catch (Exception ex) 

^ ex.printStackTraceO ; 
} 

} 

} 

void data^acti on (ActionEvent e) 

^ int returnval = spectrumChooser.showOpenDialog( this ); 
if (returnval == DFileChooser.APPROVE_OPTiON) 
{ 

try 

Page 3 



parks_et_al_12.txt 



} 
} 

} 



^ Map. Entry entry = (Map.Entry) it.nextO; 
system. out. printC entr^.getKeyC) ); 
system.out.printC " = ); , ^. 
system. out. printlnC entry. getvalueO ), 

intnn event = new int[fcs.getTota1 ()] □ ; 
Listiidei^e^Stor Imi = fcs.getListM9deiteratorO ; 
for fint i = 0; Imi .hasNextO ; ++■>) 

eventCi] = (int[]) Imi .nextQ .cloneQ ; 
int total = fcs.getTotaiQ; 

catch (Exception ex) 

^ ex.printStackTraceO; 
} 
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package edu . Stanford . f acs . 1 ogl i ke ; 

mport java.awt.*; 
mport lava.awt.event 
mport ^ava.io.BufferedReader; 
mport lava.io. File; 
mport lava.io.FileReader; 
mport java.io.lOException; 
mport java.util .ArrayList; 
mport java.util .Iterator; 
mport java.util .Map; 
mport java.util .stringTokemzer; 

import lavax. swing.*; 
import javax. swing. table.*; 



import org. i sac.*; 



fick 

* <p>Title: Log Like Data Transform</p> 

* <p>Description: </p> ... , ^ 

* <p>copyriqht: Copyright (c) 2002</p> 

* <p>Company: Stanford university</p> 

* (aauthor Wayne A. Moore 

* ©version 1.0 
V 

public class converterFraine extends 3 Frame { 
private JPanel contentPane; 
private DNienuBar jMenuBarl = new DMenuBarC;, 
orivate DMenu jMenuFile = new 3MenuC;; 
pMvati JMenultem jMenuFileExit = new DMenuitemO; 
orivate :Menu iMenuHelp = new DMenuU; 
pMvatl DMenultetn jMenuHelpAbout = new DMenultemQ; 
private JToolsar jToolear = new jToolBarQ; 
private JButton iButtonl = new JButtonQ; 
private DButton ]Button2 = new DButtonQ; 
private JButton jButtonB = new DButtonU; 
private imageicon imagel; 
private imageicon imageZ; 
private imageicon imageS; u 
orivate JLabel statusBar = new JLabelCJ, . , 
Sr vate BorderLayout borderLayoutl = new BorderLayoutO ; 
pMvSl JScrollPane jscrollPanel = new DScrollPaneQ ; 
private JTable spectrum_table = new ^TableQ; 

"user. home", ) ); ^ T^-^^r^- 

nrivate JMenultem jMenulteml = new DMenultemU, 
private DMenuitem jMenultemZ = new DMenultemCJ ; 

//construct the frame 

public converterFrameO { ...^..^ 
enableEvents(AV\fTEvent.wiNDOW^EVENT.MASK) . 

try { 
jbinitO; 

catch(Exception e) { 
e.printStackTraceO ; 

} 

//component initialization 

private void jbinitO throws Exception { 

imagel = new 
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imagelcon (edu . Stanford . f acs . 1 ogl i ke . converterFrame . cl ass . getResource("openFi 1 e . gi f ") 

); . 

lmage??on (edu?I?anf ord . f acs J ogl ^ 

)); . 

imagelcon (edu"stanf 0 rd . f acs J ogl i ke . c^^^ 

//set iconimage (Tool ki t . getDef aul tTool ki t Q . createimage (ConverterFrame . cl ass . getResou 

rce("[Your icon]"))); 

contentPane = (D Panel) this.getContentPaneO ; 
contentPane . setLayout (borderLayoutl) ; 
this.setsize(new Dimension (400, 300)); 
this.setTitleC'Log Like Data Converter ); 
statusBar.setText(" "); 
iMenuFi 1 e . setText ("Fi 1 e") ; 
iMenuFileExit.setTextC'Exit"); 

] MenuFi 1 eExi t . addActi onLi stener (new Acti onLi stener () { 
public void actionPerformed(ActionEvent e) { 
jMenuFi 1 eExi t_acti onPerf ormed(e) ; 

iMenuHel p . setText ("Hel p") ; 

iMenuHel pAbout. setText ("About") ; . . r 

3MenuHel pAbout . addActi onLi stener(new Acti onLi stener () I 
public void actionPerformed(ActionEvent e) { 
jMenuHelpAbout_actionPerformed(e); 

iButtonl.setlcon(imagel) ; 
iButtonl.setToolTipTextC'Open File ); 
3Button2.setlcon(image2) ; 
jButton2.setToolTipText ("Close File ); 
iButton3.setlcon(image3) ; 
] Button3 . setToolTi pText ("Hel p") ; 

menuspectrum. setText ("spectrum") ; . . • ^ ^^ r 

menuspectrum. addActi onLi stener (new java.awt .event .ActionLi stener (J l 

public void acti onPerformed (Acti onEvent e) { 

^ spectrurTL.action(e) ; 

]Menuiteml.setText("Data") ; . ■ . • ^ 

jMenuiteml. addActi onLi stener (new ] ava.awt. event. Acti onLi stener c; 

public void acti onPerformed (Acti onEvent e) 

data_action(e) ; 

} 

}) ; 

iMenuitem2 . setText("Save") ; 

iTool Bar . add ( j Buttonl) ; 

iToolBar.addqButton2) ; 

]ToolBar.add(]Button3) ; 

iMenuFi 1 e . add(menuSpect rum) ; 

iMenuFi 1 e . add(jMenulteml) ; 

iMenuFile.add(qMenultem2) ; 

IMenuFi 1 e . add(qMenuFi 1 eExi t) ; 

IMenuHel p. addQMenuHelpAbout) ; 

]MenuBarl.addqMenuFile) ; 

jMenuBarl.addOMenuHelp) ; 

thi s . set3MenuBar(jMenuBarl) ; 

contentPane . add ( jTool Bar , BorderLayout . NORTH) ; 

contentPane . add(statusBar , BorderLayout . SOUTH) ; 
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contentPane . add ( j Scrol 1 Panel , BorderLayout. center) ; 
j scrol 1 Panel . getvi ewport () . add(spect rum^tabl e , null); 

//File I Exit action performed . ^ n r 

public void jMenuFileExit_actionPerformed(ActionEvent e) i 

System.exit(O); 

//Help I About action performed ^ ^ s 

public void jMenuHelpAbout_actionPerformed(ActionEvent e) { 

//Overridden so we can exit when window is closed 
protected void processWindowEvent(WindowEvent e) { 

super . processwi ndowEvent (e) ; ^ ^ 

if (e.getiDQ == wi ndowEvent .window_closing) { 
jMenuFi 1 eExi t^acti onperf ormedCnul 1 ) ; 

} ^ 

private static final String [] stringArray = new String[0]; 
private String [] sensor^names; 

void spectrum_action ( 
ActionEvent e ) 

^ int returnval = spectrumchooser.showOpenDialogC this ); 
ifCreturnVal == jFileChooser .APPR0VE_0PTI0N) 
{ 

try 

^ BufferedReader br = new BufferedReaderC new FileReaderC 
spectrumChooser.getselectedFileQ ) ); 
br. readLineC) ; 

br. readLineO ; . . , . j. ■ "\^" ^. 

StringTokenizer st = new stringTokemzer( br . readuneQ , \t ), 
ArrayList al = new ArrayListQ; 
while (st .hasMoreTokensC)) 

al.addC st .nextTokenO ); ^ ^ ^ ■ * ^ 

String[] columns = sensor^names = (String[]) al.toArrayC stnngArray ); 
Doublet] [] rows = new Double[columns.length] [] ; 
for (int r = 0; r < rows. length; ++r) 

^ Doubled row = new Double[columns.length] ; 
rowsir] = row; , . „^ „ . 

st = new stringTokenizer( br. readLineO , \t ); 
for (int c = 0; c < columns. length; ++c) 

^ row[c] = new Double( Double. parseDouble( st .nextTokenO ) ); 
} 

spectrum_table.setModel( new DefaultTableModel ( rows, columns ) ); 
catch (Exception ex) 
ex.printStackTraceO ; 

} 

} 

} 

private void interpolate ( 
doubled func, 
int i , 
int n, 
double X ) 
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int m = n / 2; 

double y = Math.sqrtC func[i] * x ); 

func[ i + m ] = y; 

if (m > 1 ) 

interpolateC func, i, m, y ); 
Interpol ateC func, i + m, m, x ); 

} 



void data_action(ActionEvent e) 

int returnval = spectrumChooser. showOpenDialogC this ); 

if(returnval == JFileChooser.APPROVE_OPTION) 

{ 

try 
{ 

FCSFile fcs = new FCSFileC spectrumChooser. getselectedFileQ ); 

for (Iterator it = fcs.getKeywordsO -entrySetQ .iteratorQ ; it.hasNextQ; 

^ Map. Entry entry = (Map. Entry) it.nextO; 
System. out. print( entry .getKeyO ); 
system. out .print( " = ' ); 
System.out.println( entry .getvalueO ); 

int n = fcs.getTotal () ; 
int m = fcs.getParametersO ; 
int[][] int_value = new int[m][n]; 
ListModelterator Imi = fcs.getListModelterator() ; 
for (int i = 0; Imi . hasMoreEventsO ; ++i) 
for (int j = 0; Imi .hasMorevaluesQ ; 
int_value[j] [i] = Imi .nextValueO ; 
for (int j = 1; j <= m; ++j) 

FCSParameter p = fcs.getParameter( ] ); 
if (p.isLogO) 

doubled antilog = new double[ p.getRangeO ]; 

antilog[0] = p.getMinimumO ; 

interpolate( antilog, 0, p.getRangeO, p.getMaximum() ); 
for (int i = 0; i < n; ++i) 

int int^in = int_value[j] [i] ; 
double value = antilogC int_in ] ; 

int int_out = java.util -Arrays, binarySearch( antilog, value ); 
System. out. print ( int_in ); 
System.out.print( " -> " ); 
System. out. print( value ); 
system. out. print ( " -> " ) ; 
System. out. println( int^out ); 
if (int_out < 0) 
int_out = -int_out; 

} 

} 

^ } 

catch (Exception ex) 
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. pri ntStackTraceO ; 
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package edu . stanf o rd . f acs . 1 ogl i ke ; 

mport qava.awt.*; 
mport lava.awt. event.*; 
mport java.io.BufferedReader; 
mport lava.io.File; 
mport lava.io.FileReader; 
mport ^ava.io.lOException; 
mport lava.util .Arrays; 
mport java.util .Arrayust; 
mport java.util .Iterator; 
mport lava.util .Map; 
mport Java.util .StringTokemzer; 

import lavax. swing.*; 
import javax. swing. table.*; 

import org. i sac.*; 



<p>Title: Log Like Data Transform</p> 

* <p>Description: </p> 

* <p>copyright: Copyright (c) 2002</p> 

* <p>Company: Stanford Uni versity</p> 

* ©author Wayne A. Moore 

* ©version 1.0 



V 



public class converterFrame extends 3Frame { 
private 3 Panel contentPane; 
private DMenuBar jMenuBarl = new JMenuBarU; 
nrivate DMenu iMenuFile = new 3MenuC;; 
brivate DMenultem jMenuFileExit = new DMenultemQ, 
orivate JMenu iMenuHelp = new JMenuC); 
lAvltl DMenultem jMenuHelpAbout = new DMenultemQ; 
private DToolBar jToolBar = new DToolBarQj; 
private DButton jButtonl = new DButtonQ; 
private DButton iButtonZ = new DButtonQ; 
private DButton DButton3 = new DButtonQ; 
private imageicon imagel; 
private imageicon imageZ; 
private imaqeicon image3; u 
nr-ivatp iLabel statusBar = new DLabelQi . 

Erivate DScrolTPane j scroll Panel = new JScrollPaneQ, 
orivate DTable spectrum_tab1e = new jTabTeQ; 

5!:{;|£fp?5^JpSsr1?lft?r^hoSl^r'^^S«!;,SkhooserC syste^.getPropertyC 

"pnWatrDMenuitem'ijMenulteml = new ^Menultem^ 
private DMenultem DMenuItemZ = new DMenultemQ, 

//construct the frame 
public ConverterFrameO { 

enableEvents(AWTEvent.wiNDOW_EVENT^MASK) ; 

try { 
jblnitQ; 

catch (Exception e) { 
e.printStackTraceQ ; 

} 

//component initialization . 
private void jblnitQ throws Exception { 
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); 

imagS Cedu" sLnf ord . f acs J ogl i ke . C^^^^ 

)); . 

imagel??S Cedu" Stanford . facs J ogl i ke . con^ 

/ /seticonimage (Tool ki t . getDef aul tTool ki t () . createimage (ConverterFrame . cl ass . getResou 

rceC'CYour icon]"))); ^ , . r^ 
contentPane = (JPanel) this.getContentPaneO ; 
contentPane.setLayout(borderLayoutl) ; 
this.setSizeCnew Dimension(400, 300)); 
this.setTitle("Log Like Data converter ); 
statusBar.setTextC' ") ; 
jMenuFi 1 e . setText ("Fi 1 e") ; 
iMenuFileExit.setTextC'Exit") ; . . 

iMenuFi 1 eExi t . addActi onLi stener(new Acti onLi stener () { 
public void actionPerformed(ActionEvent e) { 
jMenuFi 1 eExi t«acti onPerf ormed(e) ; 

iMenuHelp.setTextC'Help") ; 
iMenuHelpAbout. setText ("About"); . . 

iMenuHel pAbout . addActi onLi stener(new Acti onLi stener C) 1 
public void actionPerformed(ActionEvent e) { 
jMenuHel pAbout_actionPerformed(e) ; 

jButtonl.setlcon(imagel) ; 
]Buttonl.5etToolTipText("open File ); 
jButton2.setlcon(image2) ; 
iButtonZ.setToolTipText ("Close File ); 
1 Buttons . seticon (i mageS) ; 
3 Buttons. setToolTi pText("Help") ; 

menuspectrum. setText ("spectrum"); , ^ , ^- ..-.^^^^.r^ / 

menuspect rum. addActi onLi stener (new java . awt . event .Acti onLi stenerC) i 
public void actionPerformed(ActionEvent e) { 
spectrum_action(e) ; 

})! 

iMenulteml.setText("Data"); , ^- . • ^ 

j Menulteml . addActi onLi stener (new ] ava . awt . event . Act i onLi s tene r O 

public void acti onPerformed (Acti onEvent e) 
data_action(e) ; 

} 

}) ; 

iMenuitem2 . setText ("Save") ; 

]ToolBar .add(jButtonl) ; 

iToolBar.add(jButton2) ; 

iTool Bar. add(] Buttons) ; 

iMenuFi 1 e . add(menuspect rum) ; 

iMenuFile.add(qMenulteml) ; 

iMenuFile.add(^Menultem2) ; 

iMenuFi 1 e . add(iMenuFi 1 eExi t) ; 

iMenuHel p . add(:jMenuHel pAbout) ; 

iMenuSarl.addqMenuFile) ; 

]MenuBarl.add(3MenuHelp) ; 

thi s . set JMenuBar (jMenuBarl) ; 

contentPane . add( jTool Bar , BorderLayout . NORTH) ; 
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rnntpntpane addCstatusBar, Bo rderLayout. SOUTH) ; 
coS^IS? PaSI : Idd ( j sc rol 1 Panel , sorderLayout center) ; 
j scrol 1 Panei . getvi ewport () . add (spectruin_tab1 e . null ) ; 

iofd'lMlS:i??erx?t!ac??onPerfor.ed(ActionEvent e) i 
system.exit(O) ; 

b]? i„t§°55e;fH;?Wo~onPerfo™edCAct^o„Eve„t O { 

//overridden so we can exit when window is closed 
protected void processwi ndowEvent (Wi ndowEvent e; i 
super, processwi ndowEvent (e) ; ^. ^^-r...^^ / 

if (e get ID O == wi ndowEvent. WIND0W_CL0SING) { 
jMenuFi 1 eExi t_acti onPerf ormed(nul 1 ) ; 

} * 

private static final string [] stringArray = new String[0]; 
private FCSFile fcs; 



private String 
private int[][. 
private double 
private double 
private double 



] sensor__names; 
int^data; 
raw^data ; 
compensated^data ; 
compensati on^mat ri x ; 



void compensate O 
throws lOException 

^ int n = fcs.getTotalO ; 
int m = sensor^names. length; 

raw.data = new double [ 15 It " , 
compensated__data = new doublet m JL n j; 
for (int i = 0; i < m; 



FCSParameter p = fcs.getParameter( sensor names [i] ), 
int[] input_data = int_data[p.getindexO-l] ; 
doubleC] output_data = raw«data[ij; 
if (p.isLogO) 

^ doublet] antilog = new doublet p.getRangeQ ]; 

antiloqtOl = p.getMinimumO; . n . 

iSterpSlateC antilog, 0, p.getRangeQ, p.getMaximumQ ); 

for (int j = 0; j < t?; ++1) ^ . ^ r-:i i- 
output_data[j] = antilogt input_datat]] J, 

} 

else 

double g = p.getGainQ; 
for (int j = 0; j < n; 
output_data[j] = g * input_data[3] ; 

} 

} 

for (int i = 0; i < m; ++i) 

^ doublet] mult = new doublefm]; 
multti] = 1; 

doublet] output_data = compensated_datati] ; 
for (int j = 0; j < n; ++j) 
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^ double X = 0; 
for (int k = 0; k < m; ++!<).,, 
X += raw_data[k][j] * mult[k]; 
output_data[j] = x; 

} 

} 

for (int i = 0; i < m; 

^ FCSParatneter p = fcs.getParameter( sensor_names[i] ); 
doubled input_data = compensated_dataLi J , 
int[] transfortned_data = new intLnj; 
if (p.isLogO) 

^ doubleC] lookup = new double[ p.getRangeO ]; 

lookup[Ol = p.getMinimumQ ; ^■ 
interpolateC lookup. 0. p.getRangeQ. p.getMaximumQ ), 

for (int j = 0; j < n; 

^ int k = Arrays. binarySearchC lookup, input_data[j] ); 

if (k < 0) . ^.^ 

transformed_dataLD] = -K! 
el se 

transforined_data[j] = k; 

} 

} 

else 

^ double g = p.getGainQ; 
^°[raisforLd_dltatj] 2' (iJt)' Math. floor ( input_data[j] / g ); 

} 

} 

} 

void spectrum_action ( 
ActionEvent e ) 

^ int returnval = spectrumChooser.showOpenDialogC this ); 
if(re?urnval == DFileChooser.APPROVE_OPTlON) 
{ 

try 

^ BufferedReader br = new BufferedReader( new FileReaderC 
spectrumChooser.getSelectedFileQ ) ); 
br.readLineO; 

sJringSR^ st = new stringTokenizerC br.readLineO. "\t" ); 
ArrayList al = new ArravListQ; 
while (st.hasMoreTokensQ) 

stringtf ?oiS»~?sor:n,™es = (stringn) al.toArrayC strin,Array ); 
Doubl e h [1 rows = new Doubl e [col utnns . 1 ength] LJ , 
for (int r = 0; r < rows. length; ++r) 

^ Double □ row = new Doubl e [col umns.l ength] ; 

rowsir] = rovy; . ^„„r^ -m-" 

St = new stringTokemzer( br.readLineO , \t j. 
for (int c = 0; c < col umns.l ength; ++c; 

^ row[c] = new Double( Double. par seDouble( st.nextxokenO ) ); 
^ Page 4 
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spectrum.tab1e.setModel( new oefaultTableModeU rows, columns ) ); 
catch (Exception ex) 

^ ex.printStackTraceO; 
} 

Lpensation_matrix = new doubleC sensor names length ][ sensor_names. length ]; 
fS7(int i=0; i.<sensor_names. length. 
compensation_matnxLT] [1] = i-u, 

} 

private void interpolate ( 
double [] func, 
int i , 
int n, 
double X ) 

double ~y"='^Math.sqrt( func[i] * x ); 
func[ i + m ] = y; 
if (m > 1 ) 

^ interpolateC func, i, m, y ); 
interpolateC func. i + m. m, x 



{ 



} 



} 

void data_action(ActionEvent e) 

^ int returnval = spectrumChooser.showOpenDialogC this ); 
ifCreturnval == DFileChooser.APPROVE_oPTiON) 

{ 



try 



) 



^ Map. Entry entry = (Map.Entry) it.nextO; 
system. out. printC entr^j.getKeyQ 
Svstem.out.printC " = ); , % 
system. out. printlnC entry. getValueQ ); 

int n = f cs . getTotal O ; 

int m = fcs.getParametersQ; 

intrin int value = new intLmJLnj; ^ /-n . 

LisiiodeliSitor imi = fcs.getListModeiteratorQ ; 
for Cint i = 0; Imi.hasMoreEventsQ; 
for (int i = 0; Imi .hasMorevaluesQ ; ++d) 
int_value[j][i] = Imi .nextvalueO ; 
for (int j = 1; j <= ++3-' 

* int[] ia = int_value[ j-1 ]; , ■ ^ . 
FCSParameter p = fcs.getParameter( 3 :) , 

if (p.isLogO) 

^ doublet] antilog = new double[ p.getRangeO ]; 
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?SrJ°p'offieC ?n?nS;?'rt^U«an,eO, p.getMaxlmu^O ); 

for (int i = 0; i < n; 

^ int int_in = ia[i]; . . ^ . -,. 

frt'^fntTolf .°jlJa:;??lU"r'rirsi^n"y«archC antilo,, value ); 

system. out. printC int_in ); 
system. out. pnntC -> )\ 
system. out. printC value__); 
system, out. pnntC -> j; 
system. out. pnntlnC int_out 
if (int_out < 0) 
int_out = -int_out; 

} 

} 

this.fcs = fcs; 

catch (Exception ex) 

^ ex.printStackTraceO; 
} 
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package edu . Stanford . f acs . 1 ogl i ke ; 

mport java.awt.*; 
mport ^ava.awt.event.*, 
mport Dava.io.*; 
mport java.uti I 

mport javax. swing.*; 
mport ^avax.swT-'^ "^^ 
mport ]avax.swi 



import iavax. SWT ng. event. , 
import ]avax. swing. table.*; 



ng. Timer; 
ng. border.*; 



import org. 1 sac. 

^**<p>Title: Log Like Data Transform</p> 

* <p>Description: </p> onn7^/o> 

* <p>Copyright: copyright, (c) ?002</p> 

* <p>company: Stanford umversity</p> 

* ©author Wayne A. Moore 

* ©version 1.0 

v 



pHvate JButton saveButton = new DButtonQ, 
private DButton jButtonS = new JButtonQ, 
private imageicon imagel; 
private imageicon imageZ; 

Si! ir«%s^i 

nrivate JButton spectrumButton = new JButtonQ, 

lAlTe DcSoBox^functionChoice = new DComboBoxQ; 

nrivate D Label jLabeU = new DLabelQ; 

I'AllTe 3Dialog\elpDialog = new 3Dia ogO; 

nrivate DPanel jpanell = new DPanelU, 

lAATe DButton-'helpDismiss = new ^ButtonQ; 

private DScrollPane jScrollPaneZ = "^w jscroll PaneU . 

brivate DEditorPane helpPane = new DEditorPaneu , 

private JPanel aboutMessage = new DPanelQ. 

private 3 Label j Label 2 = new ^LabelQ, 

nrivate 3 Label i Label 3 = new 3Label(); 

IrYAll GridLayJut gridLayoutl = new GrndLaycutO ; 

private 3 Label j Label 4 = new 3LabelO, 

private Component coyponentZ; 

private 3Label jLabelS = new 3LabelO, ^ 
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private D Label i Label 6 = new JLabelQ; 
private Border borderl; 
private Component components; . 
private 3 Label i Label 7 = new JLabelQ; 
private J Panel ] Panel 2 = new JPanelC); 
private DLabel statusBar = new JLabelQ; 
private JProgressBar progressBar = new jProgressBarU , 
orivate BorderLayout borderLayoutS = new BorderLayoutQ ; 
private BorderLayout borderLayout2 = new BorderLayoutQ; 
private JMenultem iMenuItem2 = new JMenuItemQ; 
private J Label j Label 8 = new JLabelQ; 
private Component componentl; 
private Component component4; 
private JPanel errorMessage = new JPanelQ; 
private GridLayout gridLayout2 = new GndLayoutQ; 
private J Label errorLabell = new JLabelQ; 
private JLabel errorLabel2 = new JLabelQ; 
private JMenu menuscale = new JMenuQ; 
private JMenuItem iMenuiteml = new JMenuitemQ; 
private JMenuItem jMenultem3 = new JMenuitemQ; 
private JMenu jMenuB = new JMenuQ; 
private JMenuItem jMenuitem4 = new JMenuitemQ; 
private JMenuItem jMenuitemS = new JMenuitemQ; 
private AxisDialog axisDialog; 

//Construct the frame 
public ConverterFrameQ 

^ enabl eEvents(AWTEvent . wiNDOW_EVENT.MASK) ; 
try 

soectrumTable-setModel ( model Parameters ); ^ n ■ ^ 

model Parameters. addTableModel Li stenerC new TableModel Listener Q 

public void tablechanged (TableModel Event e) 

^ if (spectrumxable.getSelectedColumnQ >= 0 && e.getColumnQ == 

spect rumTabl e . get Sel ectedcol umn Q ) 
nomogram_actionO ; 

spectrumTable.getTableHeaderQ.setReorderingAllowedC false ); 
spect rumTabl e . setcol umnsel ecti onAl 1 owed( true ) ; 
SDectrumTable.setRowselectionAllowedC false ) ;. ^ _ ^. 

spect rumTabl e . getcol umnModel Q , getsel ecti onModel Q . setSel ecti onModeC 
ListSelectionModel .SINGLE^SELECTION ); ^ J -,. . ^ r 

spect rumTabl e . getCol umnModel Q . addCol umnModel Li stener C 
new TableCol umnModel Li stener Q 

^ public void col umnSel ecti onChangedC ListSelectionEvent Ise ) 
^ nomogram_actionQ ; 

public void columnMarginChangedC changeEvent ce ) 

public void columnMoved( Tablecol umnModel Event tcme ) 
{ 

public void columnRemovedC TableCol umnModel Event tcme ) 
^ Page 2 



public void columnAddedC labl eColumnModel Event tcme ) 
LUoialog = new AxisDialogC this, "Logicle scale", false ); 

helpDialog.packO; 

catch (Exception e) 

^ e.printstackTraceO; 
} 

teMd^Si'ntUf'^^^^^^^^ -ception ( 
xIgSU?s?anford.facs.logli.e 

LgSU"sLnford.facs.log1i.e.onverterPr^^^ 

images = new iogiike.ConverterFrame.class.getResource("help.gif ")) ; 

rceC" tYour icon] ))) ; ^ . g getcontentPaneO ; 

iMenuFi 1 e . setrext ( Fi l e ) ; 

jNienuFi 1 eExi t_acti onPerf o rmed Ce) . 

} 

iMinuHelp.setTextC"Help';) ; . 

>enuHelpAbout.setText( About ..^ ActionListenerQ { 

'TuSfn'ord'aftfoSrM e) ( 

about_action(e); 

} 

^ public void actionPerformedCActionEvent e) 

^ data_action(e); 
} 

lS«rn:ld5«t?afs??n4r(ne«java.a«.event.Act1on^ 
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public void actionPerformedCActionEvent e) 

^ save_action(e) ; 
ifBrt°o"n:lft^S^^^^^^^^^^ -tended PCS .He"); 

^ public void actionPerformedCActionEvent e) 

^ help_action(e); 

} 

i Bitton3 . setTooiTi pTextC'Hel p") ; 
ipectrumMenuitem.setlconCnew ^ 

spectrum_action(e) ; 

} 



openMenultem.settnduicu^.---.. > 

openMenuitem.setlconOmagel), . . ^ 

^ public void actionPerformedCActionEvent e) 

^ data_actionCe); 
} 

saiewenultem . setEnabl ed Cf a1 se) ; 
saveMenuitem.seticonCimagez;, , 
saveWenuitem.setTextC save AS . 
saveMenultem . addActi onLi stene r ^new j * v 
^ public void actionPerformedCActionEvent e) 

^ save_actionCe) ; 

f uncti onchoi ce . !!^5l"]:J"';|;ii z2 (new Di mensi on (142 , 27)) ; 
functionChoice.setPreferredSize^ne ^^^^ Transformation ), 

^ public void actionPerformedCActionEvent e) 

If 

function_action(e) ; 

fprc?SK5;lS™nir«u"L<>ad spectrum Matr,x ); 
spectrumButton.setlconCnew _^„--o„pceC"spectrum. jpg"))) ' . ^ 
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^ public void actionPerformedCActionEvent e) 

^ spectrum_action(e); 
} 

cptPreferredSizeCnew Dimension (17, 17)); 
tlllroi-islirdSlS^nfis?!^^^ java.awt.event.ActionustenerO 
^ public void actionPerformedCActionEvent e) 



^ help_dismiss(e) ; 
} 



public void actionPerformedCActionEvent e) 

^ help_actionCe); 
} 



hlioPane setPreferredSizeCnew Dimension (400, 400)); 
^I^be•^e^«hffiirrl^ <head>\r\n\r\n </headAr\n <bodyAr\n <P>\r\n 
Logicle Data^converter^;^^. ,/b,dyAr\n</html>\r\n") ; 

iboS^iesslge . setLayout Cgri dLayoutl) ; 
gri dLayoutl . setCol umns CD ; 
?[iSef4°feiHrriSntilil^gnmentCswingCon 
1 Label 4 . setlconTextGapCO) ; 

Tiabel4 setTextC January 2003 ), .,..,„„ , -in^')- 

UabelS.seticonTextGapCO); . gg^rd of Trustees"); 

iLabelS.setTextC'Copynght 2002, zuua, oy , 

ttSlll:l»CvKfj4und Stanford :r. universit/^ 
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jLabel7.setText("David R. Parks"); 

fpfnelfsriSJoftCborS^rLayoutB) ; . 

gri dLayoutZ . setcol umns W , 

erroPLabelZ.setTextC iLabellO ). 
menuscal e . setEnabl ed (Tal se) , 
menuscal e . setText ( Seal e ) , 

* public void actionPerformedCActionEvent e) 



^ view_action(e); 



} 

isSiifdiss^^^^^^ 

^ public void actionPerformedCActionEvent e) 
^ new_scale_action(e); 



} 

}); 



|rlriS5:rd5i^^^^^ 

^ public void actionPerformedCActionEvent e) 
^ automati c_acti on_l_neg Ce) ; 



} 

}); 



jM^nuS . setTextC-AutomatiO ; . 
fMrnSrm'sllfdlS^n^fstfnlKn^ew ^ 
^ public void actionPerformedCActionEvent e) 
^ automati c_acti on_5_neg Ce) ; 



} 

}); 



ffiri?e:i:ifdi^^^^^ 

^ public void actionPerformedCActionEvent e) 
^ automati c_act i on_0_neg Ce) ; 



} 

}); 



iTool Bar . add Cspect rumButton , nul 1 ) , 
iTool Bar . addCopenButton) ; 
iToolBar.addCsaveButton) ; 
iToolBar.addCjButtonS); 
ijool Bar . add Q Label 1 . nu] 1 ) : , . . 
iiool Bar . addCf uncti onChoi ce , nul 1 ) . 
]ienuFi 1 e . add Cspect rumMenultem) , 
^MenuFile.addCopenwenultem), 
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iMenuFile.addCsaveMenuitem) ; 
i MenuFi 1 e . addseparator Q ; _ 
jMenuFi 1 e . addOMenuFi 1 eExi t) , 
iMenuHelp.add(3Menultem2), 

iMenuBarl.add(3MenuFi1e); 
iMenuBarl.add(menuScaTe), 
]MenuBarl.add(jMenuHelp); 

thi s . setDMenuBar (DMenuBarl) • ^ ^ . NORTH) ; ^ 

contentPane.addpToolBar, Borae ^^^^^ 
contentPane . addOscrol 1 panel .Boraer^ ^ y^^^^ 

contentPane -addO Panel 2 . BorderLay ^^^^^^^ . 

iPaSeU.addChelpDismiss, nullK ,1 2 BorderLayout. CENTER); 

aboutwessage . add ( j Label 2 , null), 
lboXssage.addCcomponent4 nu^ 

aboutwessage . add(] Label 8 , nul U , 
aboutMessage.addCcomponentl, nul U . 
aboutMessage . add Q Label 3 , nul 1 ) , 
aboutwessage . add 0 Label 7 , nul 1 ) , 
aboutwessage . add (] Label 4 .nul 1 ) . _ 
aboutwessage. add(component2nu^^ 
aboutwessage . add (3 Label 5 . nul I ) , 
aboutwessage . add p Label 6 , nu 1 1; , 
aboutwessage. addCcomponentS, null) , _ 
errorwessage . add(er rorLabel 1, nu 1 U . 
Irrorwessage . add (er rorLabel 2 , nul 1 ) , 
menuScale.addaMenuiteml); 
tnenuscale.addgwenultem3) , 
tnenuscal e . add(3Wenu3) , 
iwenu3 . addCjwenultemb) ; 
]Wenu3.add(iwenuitem4); 
]Wenu3 .add(]wenultem5) ; 

ki^c Uy5«S:i??eKfflonPerfo™edCAct1o.Event e) I 

system.exitCO) ; 

h"'^' lHS™''L'Se'?an''lx1?™hen windm. 1s closed 
jMw5FileExit-actionPerfonnedCnull), 

i iiii B s'-s-i;?:;-' ■ ■" 

private FCSFile fcs; 
private int f csTotal ; 

private int[]t] onginalData; 
private intG^p newOata;^^ 



private doub 
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r.r-i v/ai-P doubl e ri [] compensatedData ; ^ 

nrivate int progress; r-icx- 

'pfrvSI t2|icte?oS iSliri'thn = new Logicl.FunctionO ; 

public class LogicleFunction 

^ double decades, minimum, maximum; 

public doubl e[] transform ( 
int column, 

FCSParameter parameter } 
^ doubleG lookup = new doubl e[ parameter. getRangeO ]; 
SSmum = parameter .getMaxirnum(5 ; 
minimum = Parameter.getMinimumQ , 
decades = parameter. getoecadest; , 

ISeligffiterlSoK; 0. lookup.Urgth, maxmum ); 
return lookup; 

} 

public void automati cScal e ( 
int column, 

FCSParameter parameter , 
double scale ) 



public final Logiclefunction logLinearsplic. = new LogicleFunction O 



public doubleG transform ( 
int column, 

FCSParameter parameter ) 
' double linear_range = model Parameter s.getParameterC sensorNames. length. 

"^"•""iV (Double isNaNC linear_range ) II linear range <= 0) 
re?Cm super. transformC column, parameter ), 

roJirinTanr=^M^^^ / ^^""^^ ' ' 

I :»/MKip( ln_range ); 
double[] lookup = new double[ parameter. getRangeO ]; 
rnfe"r?i?it;(1SS: 0, lookup. length, maximum ); 

int solice = (int) Math.roundC 2 * parameter getRangeO / ln_range ), 

if\sp?ice > 0^&& splice < parameter. getRangeO - D 

^ double delta = (lookup [splice.l] - lookup[splice-l])/2; 

double base = lookup [splice]; 

""^^lii I gise^MSVife'- delta; 
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return lookup; 

} 

public final LogicleFunction bi exponents al = new LogicleFunction 0 
i -. I, 1^ h« suitable for immunofluorescence on lymphocytes 

55^Lr-/?°o"le l%>mlf Ar types of data 
double useroecades = 4.5; //i-J- 

public double □ transform C 
int column, 

FCSParameter parameter ; . , 

t double »1dtK - H^elPara^eters.getPara^aterC sensorKa.es.U.gth. column ), 

^*r^rrn'%ip"e?"^™sSri^C column, parameter 

double extra - ^odel.ara^eters.getParaneterC sensorNa^s. length . 1, column 
if CDouble.i5NaN( extra ) II extra < 0) 

rF'2i'^h^'U'.ro^( extra ' parameter.getRangeO / Cextra . decades) 

^' if (zero > 0) - 

^ extra - zero * (decades . extra) / parameter .getRangeQ ; // so extra falls 

exactly 02^^5f3"^lexra'*'parameter. getRangeQ / zero; 

} 

// width /= decades; //l-O 
width /= 2 * decades; //l-l 

maximum = Parameter .getMaxi mump; 

S?S?:S/^^:Ji?S=S7Satb^25jc po^f^fvll^ange ) ; 

doubled positive - new doublet parameter.getRangeO ]; 

?;re?ii?Sl(°piinlve, 0, positive.length, Math.expC positive.range ) y. 

?Sfr;o?i?U%ektive, 0, negative. length, Math.expC -negative.range ) ) ; 
g^V ;^^alJh?Sy('^pS^i?«;r/SSi?ve.range) • Cwidth . extra / 

decadesy;.^^ . j < negative. length; 

negativetj] *= s; 
s = positiveCzero] - nega^^^^^^ -j) 
^^posiJivetj] "minimlm *'(positive[j! - negative[:] - s) . 
// reflect around zero 
//or (int j=0; 3 < zero; --3)^^^^ ^ 
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positiveCj] = -positive[2 * zero - j]; 

return positive; 

} 

public void automati cScal e ( 
int column, 

FCSParameter parameter, 
double scale ) 

II ^ double minimum = parameter. getMinimumQ; //l.O * i„ m ^. 

double minimum = parameter .getMaximumQ / Math.expC userDecades * ln_10 

//l-l , . . s 

if (scale <= 0 I I scale <= mi m mum) 

^ model Parameters. setParameterC 0, sensorNames. length, column ); 
model Parameters. setParameterC Double. NaN, sensorNames. length + 1, column 

else 

// ^ double decades = Math.logC scale / parameter. getMinimumQ ) / ln_10; 

^^^'^ double decades = Math.logC scale / minimum ) / ln_10; //l-l 

model Parameters. setParameterC decades, sensorNames. length, column j; 
// model Parameters. setParameterC decades, sensorNames. length + 1, column j; 

^^^'^ model Parameters. setParameterC Double. NaN, sensorNames. length + 1, column ); 
/A.l J 

} 

}; 

public final LogicleFunction hyperbolic = new LogicleFunction Q 

public doubled transform ( 
int column, 

FCSParameter parameter ) 

^ double linear^range = model Parameters. getParameter( sensorNames. length, 

column ); 

if (Double.isNaN( linear_range )) 

return super .transform( column, parameter ); 

maximum = parameter. getMaximumO; . , , ^ , i. n 
double seal e_f actor = linear_range / (Math.E - 1/Math.Ej; 
double ln_range = Math.logC maximum / seal e_f actor ) + 2; 
decades = ln_range / ln_10; 
minimum = maximum / Math.expC ln_range ); 

doubled lookup = new double[ parameter. getRangeC) ]; 
lookupCO] = 1 / Math.E; ^ , , , x ^ ^ 

interpolateC lookup, 0, lookup -length, maximum / scale_factor j; 

for Cint i = 0; i < lookup. 1 ength ; 
lookup [j] = Clookup[j] - l/lookupCj]) * seal e_f actor; 

return lookup; 

} 

}; 

public final LogicleFunction logplus = new LogicleFunction C) 

Page 10 
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^ public doubleG transform C 



int coiumii, ^ 

KSParameter parameter ) „„„,3^terC sensorNames.length, column 

double para™eter2 - ^delPara^ters.getPara^eterC 

"'"^''rrr5^%irer"tr55;ro?fc;^co|L parameter ) ; 
if (MSble.isNaN(^parameter2 » 

:s.per.ransW ^ ^ ^^^^^^^^^^^^ , / 



} 



parameterl.*= h?' -v-^-ioSkup. length; 
*°?o^Kui "X-'pari^"?! .'i . Parametera, 

return lookup; 



public void automaticscale ( 
int column, 

FCSParameter oarameter, 



..al to.cle.„c.o. logPlu— = .o.cle.„c.o„ 0 



^ public doubleG transform ( 
int col umn , 



int column, V 

FCSParameter parameter j <;pnsorNames. length, column ); 

"iferSSHr^Sc^^rpUamet^^ 

for Cint i = O; j < lookup. length; 
^ lookupCj] -= constant; 

return lookup; 



} 

Dublic void automaticscale ( 
int column, 

FCSParameter oarameter, 
Hoiihle scale ) 



FCSPardiiicLci K 



} 
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private class Model Parameters 
extends AbstractTableModel 

^ private double[][] parameters = new double[0][]; 
public int getRowCountO 
^ return parameters. length; 
} 

public int getColumnCountO 
^ return sensorNames. length; 
} 

public String getcol umnName ( 
int columnindex ) 

^ return sensorNames [ columnindex ]; 
} 

public Class getcol umnCl ass ( 
int columnindex ; 

^ return super. getcol umnCl ass ( columnindex ); 

} 

public boolean i scel 1 Edi tabl e ( 
int rowindex, 
int columnindex ) 

^ return true; 
} 

public Object getvalueAt ( 
int rowindex, 
int columnindex ; 

^ double value = parameters[ rowindex ][ columnindex ]; 
if (Double. isNaN( value ); 
return null; 

return string.valueof ( value ); 

} 

public void setvalueAt ( 
object avalue, 
int rowindex, 
int columnindex) 

^ double dvalue; „..„i,or^ 

else if (avalue instanceof String) 
{ 

try 

^ dvalue = Double.valueOf( (string) aValue ) .doublevalueO ; 

catch (Number FormatExcepti on ex) 

^ dvalue = Double. NaN; 

J page 12 
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r&il.uSalc-raiue) rowlndex < sersorNames. length) 

if ( rowlndex == columnindex) 

dvalue = 1; 
else 

dvalue = 0; 

} 

Dublic void setParameters ( 
double[][] parameters ) 

^ thi s . parameters = parameters; 
this.fi reTableStructureChangedO ; 

} 

public double getParameter ( 
int rowlndex, 
int columnindex ) 

^ return parameters [rowlndex] [columnindex] ; 
} 

public void setParameter ( 
double value, 
int rowlndex, 
int columnindex J 

^ if (Double.isNaNC value )) ^^^„„n^n,\e^yc •)■ 
setvalueAtC null, rowlndex, columnindex J, 

^^KtvalueAtC new DoubleC value ). rowlndex. columnindex ); 

} 

. -. mnrtol Parameters = new Model Parameters Q ; 

nrivate Model Parameters . mode iParameiier > 

^^.fe'xpS^?y^'«i!OTc^=^o, imear splice", ■'log (x . cV, "x . a log x 
* priUte LogicleFurctionn logicleFunctiM , ,„s„nstant, logPlus 1; 

lrvts°rfi;MirxSi;^itesir^new'D^^^^^ 

'"ijlSra^srr^iiJ 1"p^i?«uKo-« ^Xil no. compatible wU. tKe 

"r^iii'ksii isSdoS"7i;3i^?MSi.o; 

private void interpolate ( 
doubled func, 
int i , 
int n, 
double X ) 

^ // exp((x+y)/2) = sqrt(exp(x)*exp(y)) 
iS'uble V=Valh.sqrt( func[i] * x ^)^^^ 
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funcC i + m ] = y; 
if (m > 1 ) 

^ Interpol ate( func, i , m , y ); 
interpolateC func, i + m, m, x j. 

> 

} 

nrivate void shuffle ( 
double □ sequence ) 

Sor (int i = 0; i < sequence. length; ++i) 

^ j = i + random. nextintC sequence. length - i ); 
double t = sequence [t]; 
sequenceCi] = sequence[3]. 
sequencer j] = t; 

} 

} 

public void invert ( 
double[][] matrix ) 

^ int row = 0, col =.0. n = matrix. length; 

int pivot[] = new int[n]; 

int row_index[] = new intCn] 

int col_index[] = new int[n], 

for (int i = 0; i < n; ++i) 

^ double big =0; 
for (int ] = 0; ] < n; ++]) 

^ if (pivotCj] != 1) 

for (int k = 0; k < n, ++k; 

^ if (pivotCk] == 0) 

^ double abs = Math.abs(matrix[j] [k]) ; 
if (abs >= big) 

^ big = abs; 
row = j ; 
col = k; 

} 



} 

++pivot[col 
row_index[i 
col_index[i 



hlr^ 'Sii°^RUu'rgumentEx«ptionC "n.atrix is singular" ): 



= row; 
= col ; 



if Crow != col) . 
fir (int k = 0; k < n; ++k) 

^ double t = matrix[row][k]; 
matrixCrow][k] = matnx[col] [k] , 
matrix[col][k] = t; ^^^^ 
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} 

double inverse = l/matnx[col] [colj , 

matrix [col] [col] =.1; 

for (int D = 0; D < n; ++3) 

matrix[col][j] *= inverse; 
for (int j = 0; D < n; ++DJ 

if (j != col) 

^ double t = matrix[j][col]; 
matrix[j][col] = 0; 
for (int k = 0; k < n; ++k) ^ 
raatrix[j][k] -= matnx[col] [k] t, 

} 

} 

for (int i = n-1; i >= O- 

if (row index [i] != col_index[i]) 
for (Tnt j = 0; 3 < 

matrix[j][col_index[i]] = t, 

} 

} 

double root ( 
doubl e b , 
double w ) 

^ double Xlo = 0; 

double Xhi = b; 

double Flo = Double. negative_infinity; 

double Fhi = w * 2 * b; 

double d = (xlo + xhi)/2; 

double Dx = Math.abs( Xlo - xm }, 

double Dx_last = Dx; 

double Fb = -2 * Math.log(b) + w*b; 
dSubll F = 2 * Math.log(d) + w*d + Fb; 
double DF = 2/d + w; 

if (w == 0) 
return b; 

for (int 1=0; i < 100; ++i) 

// ^ system. out. println( "bisection " + Dx ); 
ox = (xhi - xlo)/2; 
d = xlo + Dx; 
if (d == Xlo) 
return d; 

else 

= page 15 
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// system. out. printlnC "Newton " + Dx ); 

double t = d; 
d Dx; 
if (d == t) 
return d; 

if (Math.absC Dx ) < lE-12) 

return d; 
Dx^last = Dx; 

F = 2 * Math.log(d) + w*d + Fb; 
DF = 2/d + w; 
if (F < 0) 

xlo = d; 
Flo = F; 

} 

else 

xhi = d; 
Fhi = F; 

} 

} 

throw new illegalstateException( "exceeded maximum iterations" ); 

} 

void make_busy () 

^ spectrumButton.setEnabledC false ); 
spectrumMenultem.setEnabledC raise j; 

openButton.setEnabledC false ); 
openMenultem.setEnabledC false ); 

saveButton.setEnabledC false ); 
saveButton.setEnabledC false ); 

progressBar.setEnabledC false ); 

} 

void make_ready () 

^ spectrumButton.setEnabledC true ); 
spectrumMenultem.setEnabledC true 

if CsensorNames != null) 

^ openButton.setEnabledC true ); 
openMenultem.setEnabledC true ); 

} 

if Cfcs != null) 

^ saveButton.setEnabledC true ); 
saveMenultem.setEnabledC true ); 
menuScale.setEnabledC true ); 
statusBar.setTextC fcs.getFileC) -getNameC) J, 

} 

else 

^ saveButton.setEnabledC false ); 
saveMenultem.setEnabledC false 

Page 16 
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menuScale.setEnabledC false ); 
statusBar . setTextC ) ; 

} 

progressBar.setvalueC 0 ); 
progresssar.setEnabledC false ), 

} 

void spectrum_action ( 
ActionEvent e ) 



i f ( re?J rnJal == DFi leChooser . APPR0VE_0PTI0N) 

^ make_busyO; 
try 

^ BufferedReader br = new Buff eredReaderC new FileReaderC 
fi leChooser. getSelectedFTleO ) 
br.readLineO; 

^Jring?oi2nPz4r st = new stringTokenizerC br.readLineO. "\t" ); 
ArrayList al = new ArravustQ; 
while (st.hasMoreTokensQ) 

doubleClC] parameters = new double[sensorNames. length + 
2] [sensorNames . 1 ength] ; 

fSJ (r'= 0; r < sensorNames. length; ++r) 

^ St = new StringTokenizerC br.readLineO, "\t" ); 

br.closeO; 

for C; r < parameters. length; . ,_>, 

for fint c = 0; c < sensorNames. length, ++0 
narametersFr] [c] = Double. NaN; 
mode??a?ame?ers.ietParameters( parameters ); 

iff^J"(r': < sensorNames. length; ++r) 

if (fcs. get parameter ( sensorNames [r] ) == null) 

^ fcs = null; 
break; 

} 

catch (Exception ex) 
^ ex.printStackTraceO; 
make_readyO; 

} 

J page 17 
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boolean get_compensation() 

^ if (spectruniTable.isEditingO) nr^-it-innr^ 
spectrumxabl e . getcel 1 Edi tor () . stopcel 1 Edi ti ng () 

int m = senscrNames. length; 
compensationMatnx = new doubl e [mj Lmj , 
for (int i = 0; i < m; ++i) 



^°^o^p2SsajionSitr^x[j]h]'=^modelParameters.getParameter( i. j ) 



try 

^ invertC compensationMatrix ); 



catch (illegalArgumentException ex) 

^ joptionPane.showMessageDialogC this, "spectrum Matrix is Singular 
Er ror " 

DOptionPane.ERROR^MESSAGE ); 
return false; 

} 

return true; 

} 

FCSFile fcs_data; 
FCSHandl er f cs^andl er ; 
lOExcepti on f cs.excepti on ; 

Timer fcs.timer = new TimerC FCS_TIMER_TICK, new ActionListener Q 

^ public void actionPerformed ( 
ActionEvent action ) 

^ ''p^S?isfBfr'setvalu'e?fcs_handler.getEventO ); 

} 

} ); 

void data_action(ActionEvent e) 

^ if (spectruniTable.isEditingO) r■^^^c^^^.i„nr^^ 
spect r umTabl e . getCel 1 Edi tor () . stopCel 1 Edi ti ng () , 

fileChooser.setDialogTitle( "Load FCS. Data File" ); 
int returnval = fileChooser.showOpenDialog( this ), 
if(returnval == jFileChooser.APPROVE_OPTiON) 
{ 

try 

^ File f = filechooser.getselectedFileO; 

fcs data = new FCSFile( f ); -, 
for~Cint i = 0; i < sensorNames.length; 

if (fcs_data.getParameter( sensorNames[i] ) == null) 

^ jOptionPane.showMessageDialog( this, SPECTRUM_ERROR, 
"Error", DOptionPane,ERROR_MESSAGE ); 
return; 

} 

inake_busy() ; ^ 
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progressBar.setMinimumC 0 ); 
progressBar.setValueC 0 ); 

progressBar.setMaximumC fcs_data.getTotal Q ); 
progressBar.setEnabledC true ); 
statusBar.setTextC fcs_data.getFileO .getNameQ ); 

new Thread ( data_thread, "FCS reader" ).startO; 
fcs_timer. start O ; 

catch (Exception ex) 

ex.printstackTraceO ; 
^ System. exit( 1 ); 

} ^ 

Runnable data_thread = new RunnableC) 
public void run () 

Thread. currentThreadQ .setPriorityC Thread.MlN_PRiORlTY ); 

fcs_exception = null; 

try 

int n = fcs_data.getTotal () ; 
int m = sensorNames. length; 

int[][] int_value = new int[fcs_data.getParametersO] [n] ; 
double[][] double_value = new double[m] [n] ; 

fcs_handler = fcs.data.getlnputlteratorQ ; 
for (int i = 0; fcs_handler.hasMoreEvents() ; ++i) 
for (int j = 0; fcs_handler.hasMoreValuesO ; 
int_value[j] [i] = fcs_handler . readValueO ; 
fcs_handler-close() ; 

for (int i = 0; i < m; ++i) 

FCSParameter p = fcs_data.getParameter( sensorNames [i] ); 
int[] input_data = int_value[p.getindex()-l] ; 
double [] output_data = double_value[i] ; 
if (p.isLogO) 

doubled antilog = new double[ p.getRangeO ]; 
antilog[0] = p.getMinimumO ; 

interpolate( antilog, 0, p.getRangeO, p.getMaximumO ); 

double[] fuzz = new double[ 1«FUZZY_BITS ]; 
fuzz[0] = 1; 

interpolate( fuzz, 0, 1«FUZZY_BITS, antilog[l]/antilog[0] ); 
shuffle( fuzz ); 

for (int j = 0; j < n; 

output_data[j] = antilogC input_data[j] ] * fuzz[ j & 
((l«FUZZy_BlTS)-l) ]; 
} 

else 

double g = p.getGainO; 

for (int j = 0; j < n; 

output_data[j] = input_data[j] / g; 

} 

} 
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original Data = int_value; 
rawData = double_value; 

catch (lOException ex) 

^ fcs_exception = ex; 



} 

}; 



LingUtilities.invokeLaterC data_compl ete ); 



Runnable data_complete = new Runnable O 
^ public void run O 

^ fcs_timer.stopP ; 
if (fcs_exception != null; 



^ ■ coi-Tpytf fcs exception.getMessageO ); 

fr^^rtlbl 2:reSfx't^( fctdata^g^^^^^^^ 

Spti onpane . showMessageDi al og C ^° I^^^j^pti onPane . error_MESSAGE ) ; 



} 

else 

^ fcs = fcs_data; 

fcsTotal = fcs.getTotalO, 

mal<e_readyO; 

} 

}; 

void save_action(ActionEvent e) 

if (!get_compensation O) 
return; 



fileChooser.setDialogXitleC "save New FCS Data File" ); 



{ 

try 



if (f.existsQ) 

* objectCl ™essag= - { f, "already exists, oo you want to replace -.t. >, 

if (DOptionPane.OieOPTlON , ^. messaae, "Overwrite", 



return; 
} 



make«busyO ; 

int m = sensorNames. length; 

progress =0; . . ^ n ^. 
progressBar.setMinimumC 0 ), 
progressBar.setMaximumC m + m ^ m^^ 
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progressBar.setvaTueC 0 )) 
progressBar.setEnabledC true ), 
statusBar.setTextC "working ), 

IrsIdlS.ieS^ex^lgmlS^c'fU.getTextSegmentC) ) ; 
new ThreadC xfrm_thread. "Transform" ).startO; 
catch (Exception ex) 

^ ex . pri ntStackTraceO ; 
} 

} 

} 

void compensate_dataO 

Snt m = sensorNames. length; 

compensatedData = new double[m] [fcsTotai J , 
for (int i = 0; i < m; 

^ double[] output_data = compensatedData[i] ; 
for (int j = 0; j < m; 

'ZJ^J^& ?nSrjlS[k]^t^Lefficient; 
sl?!;gul[iliiies.invokeLater( progress_made ); 

} 

} 

} 

Runnable xfrm_thread = new RunnableO 
^ public void run () 
try 

compensate_data() ; 

int m = sensorNames. length; 
newData = new int[m] [fcsTotal J ; 
for (int 1=0; i < m; 

^ FCSParameter oldp = f"-H^jrSa?aL?rrO ^' 
prqparamGter newp = fcs_clata.aaclParanieteri,j , 

newp . setBi t s ( ol dp . getBi ts () ) ; 

SS;g;U«s?'^?drg«SS%k). oldp.getMini«.0 ); 
newp.setGain( oldp.getGainO ;; 

doubleP] input_data = compensatedpata[i] ; 
int[] transformed_data = newData[i]; 
if (newp.isLogO) 

^ ,i«..Mon lookuD = loqicle.transform( i, nevyp ) ; 
Sp feWcades? logi?le. decades, loaicle. minimum ); 
for^int j = 0; J< f«Tot^al^; ^ 



Parks_et_al_15 . txt 
t int k = Arrays. binarySearchC lookup, input_data[j] ); 
^^transformed_clata[j] = -k-1; 



} 

} 

else 



6l S6 

transfortned_data[j] = k; 



} 
} 

}; 



^ double g = newP-9^^^^i"^iL,i. 

^^Jr^ifforJeLSUj] ="(TSf Mat^riundC input_dataCj] 

} 

sl?!igU^ilities.invokeLater( progress_made ); 

} 

S!lt??n?fs^nvSkrLlter( xfroucomplete ); 

catch (Exception ex) 
^ ex.printStackTraceO; 



Runnable progress_made = new Runnable O 
^ public void run Q 
^ progressBar.setvalueC progress ); 

};^ 

Runnable xf rnucomplete = new Runnable 0 
^ public void run O 

^ progressBar.setMinimumC 0 ); 

rro^g'etslt^-.lftSiSL? 'f^s_data.getTotalO ); 

?S?S^se?feMVa?a%UpileO .getNa.eO ) ; 

new ThreadC save_thread, "FCS writer" ).startO; 
fcs_timer.startO; 

} 

}; 

Runnable save_thread = new RunnableQ 

^ public void run () 

^ Thread.currentThreadO.setPriorityC Thread. MIN_priority ); 
fcs_exception = null; 

^ for (int j = 0; j < original Data length; 
^ page 22 
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^°fcslhand^er.Sritivalue( newDatablL^l )• 
fcs_handler.closeO; 
catch (lOException ex) 

^ ex.printStackTraceO; 

fcs_exception = ex; 
LingUtilities.invokeLaterC save_co.plete ); 

} 

lunnable save.co.pl ete = new Runnable 0 
^ public void run O 

make_readyO; 

} 

void function_action(ActionEvent e) 

'^^^Wfbf^^^ 

n- 1 < sensorNames. length; 

} 

void nomogranuactionO 

t .„.„.nxnde. - sP=ct«Tab.e ,«S.Uc«acolu.nO , 

if (columnindex < 0 1 1 res 

return; 
FCSParameter p; 

'^'p = fcs.getParameterC sensorNames [columnindex] ); 

catch (lOException ex) 
{ 

return; 

} 

if (p.isLogO) 

^ LogicleFunction function = logicleFunctions[ 

^--^sska-^rfgSM^ p 
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doubleG left = new doubleC p.getRangeQ ]; 

left[0]= function. mini mum; function.maximum ) ; ^ut r^i • 

'^aSsSialg.no^opeUtweCiltO] = (double) (-k-U / Cdouble) left.length; 
''"isDialog.nonH.Relative[1][0] - (double) k / (double) left.length; 
k = Arrays. binarySearchC right, x ); 

''aSs5i2^o,.no»oRelative[i][l] - (double) (-k-1) / (double) right.lengtb, 
''KisDialog.no™oRelative[1][l] = (double) k / (double) right.length; 
axi sDi al og . repai nt Q ; 

} 

else 

^ double g = p.getGainQ; 
} 

} 

void help_acti on (Action Event e) 

^ helpDialog.setLocationRelativeToC this ); 
helpDialog.setVTSibleC true ), 

} 

void help_dismiss(ActionEvent e) 

^ helpDialog.setvisibleC false ); 
} 

void about_action(ActionEvent e) 

^ ,OptionPane.showMessageDialog( this. aboutMessage. "About Logicle". 
DOptionPane.PLAiN_MESSAGE ); 
public void tablechangedCTableModelEvent e) 

implemented.") ; 

void view_action(ActionEvent e) 

'^^^l?I^^S^^g4Ul^^Ei^??r^0.stopCellEditing(); 

axisDialog.setvisibleC true ); 
axi sDialog.toFrontQ ; 

} 

void scale_action(ActionEvent e) 

' ^^^fc?ffle^.9etS?i» . stopCelUditi.9() ; 
if (fcs « null) p^g^ 24 
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return; 

fiUcKooser.setPialogxitleC Scale" ); 



f-ileChooser.setDialogTixie*. aavc u^y 

String fn = fcs.getFileO .getNameO ; 
irCfnAoUpperCaseO.endsWntK ^^^^ ). 
f^'^lclolle^^^^^^^^^ fcs.getPi^eO.getParentMleO. fn 



{ 

try 



if (f.existsO) . ,,. 

« Obj.ct[] message = ( f. ^ " 

if (DOptionPane.OieOPTION , .j... messaae, "Overwrite", 

return; 

} 

int c = spectruraTable.getselectedColumnO; 
FCSpHm4ter p = fcs.getParameterC sensorNames[c] ): 

pw.printlnC scale[i] ); 
pw.closeO ; 

catch (Exception ex) 

^ ex.printStackTraceO ; 
} 

} 

} 

void new_scale_action(ActionEvent e) 

' ^^^KKfe'.geiS^ 
if (fcs == null) 
return; 

fileChooser.setDialogTitleC "save Logicle Scales" ); 

String fn = fcs.getFileO ,getNatneO; 

if (fn.toUppercaseQ.endsWithC -FCS IJ „ -^^^„^^^Q ); 

f!4hoSehSa^lk2i4?4rr^^leC fS%;l?mo"gUparent.ileO. ^ 

{ 

try 



parks^et^al. 15.txt 

if (f.existsQ) 



* Objeott] .es.a,e = { f. "alrMdy exists. 00 you want to replace ,t. 
if (30ptionPane.OK_OPTlON , , t^. message, "overwrite". 



return; 



^ OK^nn^-wr-itPrf new FileOutputstreamC f ) ); 

^ FCSParameter p - fcs.getParameterC sensor»ames[i] ); 

p„.println( spectru»,Na.eV r«ensor«a™est1]V] <".sensorNa™esCi]. 

I^Eg-s^aW^^^^^^ 

pw.printlnC scaleL;]] )' 
pw.prnntlnC p.getMaximumQ ), 



pw. close C) ; 



} 
} 

} 



catch (Exception ex) 
^ ex.printStackTraceO ; 



private double autoscaleParameter; 
void autoscale_action O 
if (!get_compensation O) 

logHcle - logicleFunctionsC functionchoice.getselectedlndexO ], 

try 

tnake_busyO; 

int m = sensorNames. length; 

progress =0; . . „^ ^ ^^ 
progressBar.setMimmumC 0 * „ 
brogressBar.setMaximumC m + m * m ), 
progressBar.setvalueC 0 ); 
progressBar . setEnabl ed( true ) , 
statusBar.setTextC "working ), 

newThreadC autoscale_thread. "Autoscale" ).startO; 
catch (Exception ex) 
^ ex.printStackTraceO; 
} 

} 

Runnable autoscale_thread = new RunnableQ 
^ public void run () p^^^ 26 



)]; 
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{ 

try 

compensate_data() ; 

int m = sensorNames. length; 
for (int i = 0; i < m; 

^ double[] 5ample_data = cotnpensatedData[i] ; 
Arrays. sortC sample_data 

III clSlt-SS'. O; sa™ple_data[scale.tatal] < 0 &S scale.total < 
''v^!5l'Hf»plLSL[Cint) Math.roundC autoscalePara.eter • scale_total 
l„gicl..auto»aticscale( i, fcs.getMrameterC sensorNamesCi] ), value ): 
S^SJiltiies.invokeLaterC progressjade ); 

} 

SSg!l""tfe" ^vSkeUterC autoscale^conplete ); 

catch (Exception ex) 

^ ex.printStackTraceO ; 
} 

} 

}; 

Runnable autoscale_complete = new Runnable O 
^ public void run () 

^ nomogram_action() ; 
make_ready() ; 

} 

pi-ivate DMenuitem jMenultem6 = new DMenuitemQ; 
void automatic_action_L.neg(ActionEvent e) 

^ autoscaleParameter = .01; 
autoscale_action() ; 

} 

void automatic_action_5_neg(ActionEvent e) 

^ autoscaleParameter = .05; 
autoscale_action(); 

} 

void automatic_action_0_neg(ActionEvent e) 

^ autoscaleParameter = 0; 
autoscale_action() ; 

} 

} 
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package edu . Stanford . f acs . 1 ogl i ke ; 

import ^avax. swing.*; 
import java.awt. ; 

''**<p>Tit1e: Log Like Data Transform</p> 

: X?op5H?hi:"^o;(?ight (c) 2002</p> 

* igiSpany: Stanford universe ty</p> 

* ©author Wayne A. Moore 

* ©version 1.0 

*/ 

nublic class DataConverter { . 
private boolean packFrame = false; 

//construct the application 

"IJi^ereJ^are^frS i ne» converterFra.e0 1 

HlJL:r?Kt^KKe^-fS^"-f-^^^ sue info, e.,. fro. their layout 

if CpackFrame) { 
frame. packO ; 

} 

else { ^ . . ^ . 
frame. vail date O ; 

framesize. height = screenSize. height, 

f.a™e.setuocationCCscree„s1ze.width - fra^Slze.widtW / 2, (screens, ze. he, 

framesize. height)/ 2); 
frame . setvi si bl eCtrue) , 

Sffli"c S'tfc void n«incstring[] ar«0 I 
^'•JliLager.setuookAndFeeUuiManager.getsysten.LookAndFeeKlassNameOJ; 

catch (Exception e) { 
e . pri ntstackTraceO • 

new DataConverter O; 



} 
} 
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package edu.stanford.facs.loglike; 

import javax. swing.*; 
import ]ava.awt.*; 

/r<p>Title: Log Like Data Transform</p> 

* <p>Description: </P>.. ?002</d> 

* ©author Wayne A. Moore 

* ©version 1.0 

//construct the application 

public DataConverterO i converterFrameQ; 

if (packFrame) { 
frame. pack U ; 

frame . val i date () ; 

}/center the window ^j^^t.getDefaultToolkitQ .getScreenSizeQ ; 

f ramesize.wiuLii ^ . , .n / t rcrrppn<;ize.hei 

la.e.set.ocat1o„C(scree„S.ze.»1dth - f r»eSlze.w,dth) / 2, (scr«.s,ze 

framesize.height) / 2) , 
frame . setvi si bl eCtrueJ , 

kli"c s?a5ic void main(string[] args) { 

catch (Exception e) { 
e.printStackTraceO, 

new DataConverterO; 

} 

} 
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<html> 
<head> 
</heacl> 
<body> 

<center><hl>Logicle Fluorescence 
Compensati on</hl></center> 

'hrllll'.T. I'Z^^n.WT"^^^ range raw data 

that are typically. 

plotted on a loganthimc scale. 

l5lrfa;^^C?%^(ron^;f Sw jluSle values or "fluorescence 

na?S?alirp"rocluces small and negative values that are 
not handled well by <b>log</b> scaling. 

:JiL?<b?S<%Wcode> data files to include various 
<i>li^e</i> representations of the data, which are intended to 
bl«er properties for small and negative values. 
</p> 

and one or two 

additional parameters. ^^«^*..,Mm 
For example for a simple two color spectrum 
<code><b>Fluor</b></code> and 
<code><b>PhyEry</b></code> 

<center> „ ... „-„ 

<table border=;;r cellpadding= 4 
cellspacing="2"> „ n „ 

<col width="2*" rules="cols > 
<col width="l*" 

<t?><th></th><th><code><b>Fluor</b></code><^^ 

PhyEry</b></code></th></tr> 

<t rxth al i gn=" ri ght " styl e=" col or : 

hlijp"><code><b>Fluor</b></code> 

CoStrol</th^<td>1.0</td><td>0.1983</td></tr> 

^[l2"colo?rbi;i"^<code><b>PhyEry</b></code> 
cSt?ol </th>;td>0 . 0088568</td><td>l. 0</td></t r> 

2yl2"cSlo?rred^>'parameter l</th><td></td><td></td></tr> 

aligS-righf style="color: red">Pararaeter 
2</th><td></td><td></td></tr> 

</tabl e> 
</center> 

2ode>!;Kn</b></code> rows are the ordinary (normalized) spillover 
coefficients for each of 
the dyes used. 

S™«.JsVari«'derendin, on the function choosen. 

YOU can load a 

new spectrum at any time. 

</P> 
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once you've loaded a spectum you 

can open an <code><b>FCS</b></code> data file 

trappj55Ha?e!''??e:"'it contains data for all the fluorochromes 

finally you can save an extended <code><b>FCS</b></code> file 
that contains the new compensated data 
as well as any previous 

gS!pSrE?5!vbivcSa:^S?^^^?Srt?Sjsfon™ed data will be 

<b>Fl Sw?o<font si ze="-2"><sup>TM</sup></f ont></b> compensated 
parameters <code><b>&l t ; Fl uor&gt ; </b></code> 

<Jode><b><PhyEry></b></code> will be preserved unchanged if 

present. 

</P> 

<Elockquote> ^ „. j ^ 

<font size="+l" style=" indent: 

JSS^'can alter the function parameters to adjust them to your 

data and you can alter . - ^ 

the spectral overlap values if you wish to 

tinker with the 

compensation. 

</font> 

</blockquote> 

</p> 

<h2>Functions 
Avail abl e</h2> 

jRe first two functions are parameterized by a 
single value that characterizes 
thp desired range (in the units r- ^l. 

encoded in the ?code><b>FCS</b></code> file) of the 

<b>bi polar ^ , . 

linear</b> region of the scale. 

^?^eK^:rsSi?S</L;^:/K^^ of the transformed data 

Segative linear boundry is the minimum value plotted, the zeroes will 
at one factor of <b><code>e</code></b> above the minimum and the 

linear boundry „j„^^/u^ 
is at <b><code>e<sup>2</sup></code></b>. 

</p> 
<h3>log 

linear splice</h3> 

tK data are transformed linearly up to the 
parameter value and logarithm! cly beyond it. 

tflX'il i§n?i5SoSs in the first derivative although not in the 
5^^°"^- page 2 
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The scale values at the boundry and above are identical to 
the normal values. 
</p> 

<h3>hyperbol i c</h3> 

jRe data are scaled ^ , , o-. ^ / ^ 

to bring the linear region to <code>|x| < l</code> 

<code><b>sinh</b></code> is taken and scaled in turn to produce 
the ordinary 

logarithm for large values. 

The <code><b>sinh</b></code> 

function is odd and linear for small values and 

exponential at large . ,.1.1^% 

values (more than a decade above the linear threshold), 

while the 

derivatives are as well behaved as they get. 
</P> 

<h3>x + a log x + 
b</h3> 

A'^normal antilog table is computed for the data then the 

values in the table 

are offset by a linear function. 

With suitable , . ^ _ . . 

parameters, the small values will be approximately linear and 
large 

values will be offset by insignificant amounts. 

The first parameter ^ , .-.n ^ j ^ 

gives the approximate linear value that will be plotted at 

one decade 

above the minimum. 

The second (optional) parameter can be used to 
displace the linearized region 
to include negative 
numbers. 

</p> 
<hr> 

<center><i><font size="-l"> 

Copyright © 2002 

by The Board of Trustees<br> 

of the Lei and Stanford 3r. 

University 

</f on tx/i ></cente r> 

</body> 
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<html> 
<head> 
</head> 
<body> 

<center><hl>Logicle Fluorescence Compensation</hl></center> 

tKis program computes compensated fluorescence data from high dynamic range raw data 
that are typically plotted on a logarithmic scale. 

The <b>linear</b> process of spectral overlap correction of raw fluorescence values 
or "fluorescence compensation" 

naturally produces small and negative values that are not handled well by <b>log</b> 
scaling. 

This program extends <code><b>FCS</b></code> data files to include various 
<b>log</b> <i>like</i> representations of the data, which are intended to 
have better properties for small and negative values. 
</p> 

T^e transformations require as input a <b>spectrum</b> matrix exported from 
<b>FlowJo<font size="-2"><sup>TM</sup></font></b> 
and one or two additional parameters. 

For example for a simple two color spectrum <code><b>Fluor</b></code> and 

<code><b>PhyEry</b></code> 

<center> 

<table border="l" cellpadding="4" cellspacing="2"> 
<col width="2*" rules="cols"> 
<col width="l*" span="2"> 

<tr><th></th><th><code><b>Fluor</b></code></th><th><code><b>PhyEry</b></code></th></ 
tr> 

<tr><th align=" right" style="color : blue"><code><b>Fluor</b></code> 
Control </th><td>l . 0</td><td>0 . 1983</td></t r> 

<tr><th align=" right" style="color ; blue"><code><b>PhyEry</b></code> 
control </th><td>0 . 0088568</td><td>l . 0</td></t r> 

<tr><th align=" right" style="color: red">Parameter l</th><td></td><td></td></tr> 
<tr><th align=" right" style="color : red">Parameter 2</th><td></td><td></td></tr> 
</table> 
</center> 

The first <code><b>n</b></code> rows are the ordinary (normalized) spillover 
coefficients for each of 
the dyes used. 

The meaning of the two parameters varies depending on the function choosen. 

YOU can load a new spectrum at any time. 

</p> 

<P> 

once you've loaded a spectum you can open an <code><b>FCS</b></code> data file 
for which the spectrum is appropriate, i.e., it contains data for all the 
fluorochromes listed. 

Finally you can save an extended <code><b>FCS</b></code> file that contains the new 
compensated data 

as well as any previous contents. 

For example, if you initially had a file with raw data for <code><b>Fluor</b></code> 
and <code><b>PhyEry</b></code>, 

then new parameters <code><b>[Fluor]</b></code> and <code><b>[PhyEry]</b></code> 
with the transformed data will be added. 

Note that the compensated values are computed independantly of <b>Flow3o<font 
size="-2"><sup>TM</sup></font></b> 

and that <b>FlowJo<font size="-2"><sup>TM</sup></font></b> compensated parameters 
<code><b>&l t ; Fl uo r&g t ; </b></code> 

and <code><b><PhyEry></b></code> will be preserved unchanged if present. 

</p> 

<Clockquote> 

<font size="4-l" style="indent : 10px"> 

You can alter the function parameters to adjust them to your data and you can alter 
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the spectral overlap values if you wish to tinker with the compensation. 

</font> 

</blockquote> 

</p> 

<P> 

If you first select one of the data columns and then choose <b>view</b> from the 
<b>Scale</b> menu a side by side comparison of the default <b>FlowJo<font 
size="-2"><sup>TM</sup></font></b> display scale and 
the <code><b>Loqicle</b></code> scale is displayed. 

If you choose <D>Save</b> from the <b>Scale</b> menu the scale for that parameter 
will be 

saved in a file that <b>Flow3o<font size="-2"><sup>TM</sup></font></b> can use to 
display the scale correctly. 

Finally, if you choose <b>Automatic</b> then <code><b>Logicle</b></code> will 
attempt to 

set the scale appropriately for the data according to the chosen criterea. 
</p> 

<h2>Functions Avail able</h2> 

<h3>Bi exponenti al </h3> 
<p> 

The bi exponential transformation uses a function of the form 
<center> f(x) = a e<sup>bx</sup> - c e<sup>-dx</sup></center> 
to transform the data. 

The first parameter is the number of decades (onto the original scale) where 

the transition from the negative to the positive exponential occurs. 

By default the same value is used to set the range of negative values shown. 

The second parameter is optional and if present, indicates the number of extra 

decades worth of space added to negative range. 

</p> 

<h3>log (x+c) </h3> 

<P> 

Adds a constant to the value and then takes the ordinary logarithm 

</p> 
<p> 

The next two functions are parameterized by a single value that characterizes 

the desired range (in the units encoded in the <code><b>FCS</b></code> file) of the 

<b>bi polar linear</b> region of the scale. 

For both functions a <b>FlowJo<font size="-2"><sup>TM</sup></font></b> plot of the 
transformed data 

the negative linear boundry is the minimum value plotted, the zeroes will fall 
at one factor of <b><code>e</code></b> above the minimum and the linear boundry 
is at <b><code>e<sup>2</sup></code></b>. 
</p> 

<h3>log linear splice</h3> 

T^e data are transformed linearly up to the parameter value and logarithmicly beyond 
it. 

The splice is made so that it is continuous in the first derivative although not in 
the second. 

The scale values at the boundry and above are identical to the normal values. 
</p> 

<h3>hyperbol i c</h3> 
<p> 

The data are scaled to bring the linear region to <code>|x| < l</code> then 
<code><b>sinh</b></code> is taken and scaled in turn to produce the ordinary 
logarithm for large values. 

The <code><b>sinh</b></code> function is odd and linear for small values and 
exponential at large values (more than a decade above the linear threshold), 
while the derivatives are as well behaved as they get. 
</p> 
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<h3>x + a log X + b</h3> 

<P> 

A normal antilog table is computed for the data then the values in the table 
are offset by a linear function. 

with suitable parameters, the small values will be approximately linear and 
large values will be offset by insignificant amounts. 

The first parameter gives the approximate linear value that will be plotted at 
one decade above the minimum. 

The second (optional) parameter can be used to displace the linearized region 

to include negative numbers. 

</p> 

<hr> 

<center><i><font size="-l"> 

copyright © 2002 by The Board of Trustees<br> 
of the Lei and Stanford Jr. university 
</font></i></center> 
</body> 



</body> 
</html> 
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